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ABSTRACT 


Two  separate  computer  codes  were  coupled  with  a con- 
strained function  minimization  code  to  produce  automated 
marine  condenser  design  and  optimization  programs  of  vastly 
different  complexity.  The  first  program,  OPCODEl,  was 
developed  from  the  Heat  Exchange  Institute's  Standards  for 
Steam  Surface  Condensers  (HEI).  The  second  program,  0PC0DE2, 
was  developed  from  the  sophisticated  ORCONl , a computer  code 
produced  by  the  Oak  Ridge  National  Laboratory.  CONMIN,  the 
optimization  program,  was  developed  at  the  Ames  Research 
Center. 

OPCODEl  was  well  verified  using  .main  condenser  input  data 
of  an  aircraft  carrier  and  a destroyer  escort.  Verification 
of  0PC0DE2,  using  main  condenser  data  of  an  aircraft  carrier, 
was  less  satisfactory  due  to  the  conservative  nature  of 
flooding  effects  on  the  outside  film  heat  transfer  coeffi- 
cient used  in  ORCONl. 

OPCODEl  is  an  excellent  design  tool  for  the  conceptual 
design  of  a marine  condenser.  Optimized  test  cases  run  with 
OPCODEl  show  that  a condenser  designed  by  the  HEI  method  is 
nearly  optimum  with  respect  to  volume. 

Test  cases  with  0PC0DE2  show  that  enhancing  the  heat 
transfer  on  the  shell-side  by  80  percent  yields  a condenser 
with  ten  percent  less  volume  than  the  unenhanced  case. 
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NOMENCLATURE 


English  Letter  Svmbols 


G.  (X) 


tube  sheet  area/number  of  tubes,  ft'^/tube 

2 

heat  transfer  area,  ft 

internal  tube  area  per  linear  foot,  ft^/ft 

2 

external  tube  area  per  linear  foot,  ft  /ft 

flow  area  for  a single  tube,  ft" 

coefficient  for  calculation  of  overall  heat 
transfer  coefficient 

tube  outside  diameter,  ft 

tube  inside  diameter,  ft 

internal  heat  transfer  enhancement  factor 

external  heat  transfer  enhancement  factor 

objective  function 

average  tube  flooding  factor 

tube  flooding  factor  for  the  n-th  tube 

fouling  factor 

material  correction  factor 

temperature  correction  factor 

tube-side  friction  factor 

cooling  water  flow,  gpm 

inequality  constraint 

tube  flow  factor 


acceleration  due  to  gravity,  ft/hr 
acceleration  due  to  gravity,  ft/sec*” 
equality  constraint 
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h 


enthalpy,  BTU/lb 


h.  - internal  film  heat  transfer  coefficient, 

^ BTU/ (hr) (ft’) (“F) 

h - external  film  heat  transfer  coefficient, 

° BTU/ (hr) (ft-)  (°F) 

h*  - external  film  heat  transfer  coefficient 

° corrected  for  inundation,  BTU/ (hr) (ft^) (°F) 

h£g  - latent  heat  of  condensation,  BTU/lb 

ICALC  - control  flag  for  COPES 

k - tube  constant,  k = s/g 

k^  - fluid  conductivity,  BTU/ (hr) (ft) (°F) 

k^  - wall  conductivity,  BTU/ (hr) (ft) (°F) 

k^  - vapor  conductivity,  BTU/ (hr) (ft) (®F) 

L - tube  length,  ft 

LMTD  - log  mean  temperature  difference,  ®F 

m - mass  flow  rate,  Ib/hr 

N ' number  of  tubes 

NAC  - number  of  active  constraints 

NCON  - number  of  constraints 

NDV  - number  of  design  variables 

n*  - number  of  tubes  in  a vertical  row  above 

the  n-th  tube 

P - absolute  pressure,  psia 

PP  - pumping  power,  ft>lbf/sec 

AP  ^ - tube  exit  loss,  ft  w.c. 

ext 

AP  ^ - tube  entrance  loss,  ft  w.c. 

ent  ’ 

AP  - sum  of  all  tube-side  pressure  losses,  ft  w.c. 

AP^  - internal  tube  friction  loss,  ft  w.c. 
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f 


1 


Q 

q 

R. 


s 

s 

T 

AT 


LM 


U 

U. 


VLB, 


VUB: 


w.  c. 


w 

X 


heat  transferred,  BTU/hr 
iteration  number 

fouling  resistance,  (f t ) (hr) (°F) /BTU 
tube  inside  radius,  ft 
tube  outside  radius,  ft 
search  direction 

outside  area  of  tube  per  linear  foot,  ft“/ft 
temperature,  “R 

log  mean  temperature  difference,  ®F 
cooling  water  inlet  temperature,  °F 
cooling  water  outlet  temperature,  ®F 
saturation  temperature,  '’F 
vapor  temperature,  °F 
wall  temperature,  °F 

uncorracted  overall  heat  transfer  coefficient, 
BTU/(hr) (ft2) (°F) 

corrected  overall  heat  transfer  coefficient, 
BTU/(hr) (ft2)  (»FJ 

overall  heat  transfer  coefficient  based  on 
inside  tube  area,  BTU/ (hr) (ft^) (®F) 

cooling  water  velocity,  ft/sec 

lower  side  constraint  on  i-th  design  variable 

upper  side  constraint  on  i-th  design  variable 

steam  flow,  Ib/hr 

water  column 

weight  of  cooling  water,  Ib/gal 
vector  of  design  variables 
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Dimensionless  Groups 
Pr  - Prandtl  number 

Re  - Reynolds  number 


Greek  Letter  Symbols 

a*  - move  parameter  in  optimization  problem 

3 - ratio  of  Ap/A^ 

parameter  in  method  of  feasible  directions 
absolute  roughness,  ft 

push  off  factor  in  method  of  feasible  directions 
fluid  absolute  viscosity,  lb  sec/ft*" 
condensate  density,  Ib/ft^ 

3 


3 

“f 


sw 


sea  water  density,  Ib/ff 
vapor  density,  Ib/ft^ 
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I.  INTRODUCTION 


A.  BACKGROUND 

In  the  last  ten  years  a revolution  has  swept  the  marine 
power  plant  industry  that  could  result  in  the  obsolescence 
of  the  marine  steam  power  plant.  The  revolution  was  caused 
primarily  by  the  use  of  the  gas  turbine  as  an  alternative 
to  marine,  and  more  recently,  to  naval  propulsion.  Gas 
turbines  brought  about  power  plants  which  were  more  compact, 
lighter,  but  less  fuel  efficient  than  the  steam  power  plants 
which  they  displaced. 

When  compared  with  the  compact  gas  turbine,  the  massive 
size  and  weight  of  the  marine  steam  power  plant,  which 
evolved  in  stride  with  the  behemoth  of  the  power  industry, 
the  stationary  steam  power  plant,  made  this  means  of  propul- 
sion less  desirable  for  naval  vessels.  Therefore,  it  has 
become  imperative  for  the  naval  engineering  community  to 
develop  a more  efficient,  more  compact,  lighter  weight  steam 
power  plant. 

As  steam  power  plants  are  durable  and  they  can  burn  a 
variety  of  fuels  - both  essential  qualities  when  considering 
the  Navy's  combat  readiness  — they  must  not  be  allowed  to 
be  overdesigned  out  of  existence.  To  make  steam  propulsion 
competitive  with  marine  gas  turbines,  advanced  concepts  must 
be  explored  in  all  areas  of  steam  propulsion.  Such  concepts 
as  pressurized  boilers,  super  critical  cycles,  enhanced 
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condenser  tubes,  and  dropwise  condensation  must  be  developed 
further.  Above  all,  overdesign  by  the  use  of  unnecessary 
safety  factors  must  be  curtailed,  and  the  minimum  safe  design 
must  be  developed  and  identified. 

B.  METHODOLOGY 

In  the  United  States  the  most  prevalent  criterion  for  the 
design  and  specification  of  surface  condensers  is  based  on 
the  "square  root  of  V"  relationship  as  developed  by  the  Heat 
Exchange  Institute  (HEI)  [1].  Using  this  method,  the  overall 
heat  transfer  coefficient  is  calculated  as  a function  of  the 
square  root  of  the  cooling  water  velocity  multiplied  by 
correction  factors  for  inlet  cooling  water  temperature,  tube 
wall  thickness  and  material,  and  fouling. 

The  HEI  method  was  adopted  by  the  Department  of  the  .Navy, 
Bureau  of  Ships  (now  Naval  Sea  Systems  Command)  for  the 
specification  of  U.  S.  Navy  condensers  by  issuing  Design 
Data  Sheet  4601-1  (DDS)  [2]  in  1953.  Henceforth  this  thesis 
will  designate  the  preceding  methodology  as  the  HEI/DDS 
method. 

With  the  advent  of  the  high  speed  digital  computer, 
numerical  methods  of  solving  complex  engineering  problems 
are  now  possible.  A computer  code  has  been  developed  to 
calculate  the  local  heat  transfer  and  thermodynamic  properties 
of  a large  surface  condenser  on  a row  by  row  basis.  Known  as 
ORCONl,  this  code  was  developed  by  Oak  Ridge  National 
Laboratory  (ORNL)  under  contract  to  the  Office  of  Saline 
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Water  during  the  period  from  1968  to  1970  [3],  The  program 
was  based,  in  part,  on  the  work  performed  by  Eissenberg  [4]. 
Eissenberg's  experimental  results  led  to  correction  factors 
on  the  basic  Nusselt  equation  to  account  for  inundation 
effects  of  tubes  within  a condenser  tube  bundle.  Addi- 
tionally, logic  was  developed  to  account  for  the  pressure 
loss  caused  by  the  steam's  passage  through  the  tube  bank 
with  the  accompanying  reduction  in  saturation  temperature; 
the  heat  resistance  due  to  the  presence  of  a noncondensable 
gas  film;  heat  transfer  enhancement  factors  on  both  sides  of 
the  tubes;  and  other  important  factors  to  yield  a program 
which  could  calculate  heat  flux,  overall  heat  transfer 
coefficient,  noncondensable  gas  concentration,  and  fifteen 
additional  parameters  on  a local,  row  by  row  basis. 

Search  [5]  has  used  ORCONl  to  perform  parametric  studies 
of  an  actual  naval  condenser.  Tube  enhancement,  the  high 
velocity  flow  allowed  with  titanium  tubes,  and  dropwise 
condensation  were  investigated.  The  penalties  paid  in 
increased  pumping  power  and  increased  cost  were  weighed 
against  the  gains  realized  with  a more  compact  and  with  a 
lighter  condenser. 

In  the  open  literature  there  is  but  one  reference  [6] 
to  coupling  a condenser  analysis  and  design  program  with 
an  optimization  procedure  that  is  capable  of  improving  a 
given  design. 


I 
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I The  case  for  utilizing  an  optimizing  design  scheme  for 

the  condenser  portion  of  a steam  propulsion  plant  can  easily 
be  made  by  re-emphasizing  the  fact  that  in  order  for  steam 
propulsion  to  remain  a viable  contender  when  the  naval 
vessels  of  the  late  1980's  are  designed,  it  must  compete 
and  succeed  in  an  area  that  is  rapidly  being  dominated  by 
gas  turbines.  All  components  of  the  naval  ship's  steam 
propulsion  plant  will  have  to  be  critically  designed  to 
ensure  that  the  minimum  design  will  still  perform  as  and 
when  required. 

C.  OBJECTIVES 

There  were  two  primary  objectives  of  this  thesis. 

The  first  objective  was  to  develop  a computer  code  based 
on  the  HEI/DDS  method  of  condenser  design,  as  this  method 
is  considered  the  industry  standard.  Coupling  the  HEI/DDS 
code  with  a numerical  optimization  program  yields  a complete 
design  package.  The  design  package  can  be  used  for  trade- 
off studies,  first  cut  analysis,  and  conceptual  design. 

The  second  objective  was  to  couple  ORCONl,  with  its 
capability  of  varying  both  internal  and  external  tube 
enhancement  factors,  with  a numerical  optimization  program 
to  produce  a more  detailed  design  program. 

These  design  tools  provide  the  Naval  architect  and  the 
Naval  engineer  with  the  means  to  optimize  size,  weight,  ? 

' . I 

j design,  and  cost  of  the  marine  steam  propulsion  plant  for  i 

I ships  of  the  1980's;  the  enhanced  design  provides  the  optimum  f 
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streamlined  design  of  a steam  plant  resulting  in  its 
reinstatement  and  continuance  as  a viable  alternative  to 
gas  turbine  propulsion. 


II.  NUMERICAL  OPTIMIZATION 


A.  BACKGROUND 

Nearly  all  design  problems  require  either  the  minimiza- 
tion or  maximization  of  a parameter  or  function.  This 
parameter  will  be  called  the  problem's  objective  function 
or  design  objective  [7].  For  the  design  to  be  acceptable, 
it  must  satisfy  a set  of  design  constraints.  For  example, 
if  an  engineer  was  designing  a piping  system  to  achieve  the 
minimum  in  required  pumping  power,  the  minimum  allowable 
flow  delivered  would  be  a meaningful  constraint.  Likewise 
a constraint  that  required  the  inside  diameter  of  the  pipe 
to  be  less  than  the  outside  diameter  of  the  pipe  would  be 
a necessity. 

If  the  problem  could  be  formulated  analytically  with  a 
great  deal  of  simplicity,  the  minima  or  maxima  could  be 
found  by  using  the  methods  of  differential  calculus  or  the 
calculus  of  variations.  However,  these  methods  would  fail 
for  all  but  the  very  simplest  of  problems. 

A numerical  method  that  would  be  satisfactory  for  rela- 
tively small  scale  problems  would  be  an  iterative  solution 
technique.  A simple  computer  program  could  be  written 
containing  a series  of  nested  iteration  loops  that  would 
vary  the  design  parameters  and  solve  the  problem  for  a 
variety  of  values  for  each  of  the  parameters.  For  small, 
easily  formulated  problems,  the  cost  in  central  processor  (CPU) 
time  would  not  be  excessive,  and  this  method  would  be 
satisfactory . 
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However,  for  a serious  engineering  problem,  this  method 
quickly  becomes  too  costly  to  pursue.  For  example,  if  the 
engineer  had  ten  design  parameters  for  which  he  wanted  to  try 
ten  separate  values,  he  would  need  to  make  ten  billion 
calculations.  If  each  calculation  took  ten  CPU  seconds, 
the  solution  would  be  available  in  approximately  3200  years! 
Thus  a rational  approach  to  design  automation  and  optimiza- 
tion is  obviously  needed. 

There  are  many  optimization  schemes  available  to  the 
engineer.  The  various  methods  fall  into  three  broad  cate- 
gories based  on  the  type  of  problem  to  be  solved: 
unconstrained  minimization,  solution  of  constrained  problems 
by  unconstrained  minimization,  and  direct  methods  for 
solution  of  constrained  problems  [8].  An  optimization 
program  based  on  the  last  method  was  chosen  for  this 
research  work. 

B.  CONSTRAI.MED  FUNCTION  MINIMIZATION  (CONMIN) 

Vanderplaats  [9]  developed  an  optimization  program, 
CONMLN,  capable  of  optimizing  a very  wide  class  of 
engineering  problems.  CONMIN  is  a FORTRAN  program,  in 
subprogram  form,  that  optimizes  a multi -variable  function 
subject  to  a set  of  inequality  constraints. 

Three  basic  definitions  are  required  [10]: 

Design  Variables  - Those  parameters  which  the 
optimization  program  is  permitted  to  change  in 
order  to  improve  the  design.  Design  variables 
appear  only  on  the  right  side  of  an  equation, 
are  continuous  and  have  continuous  first  derivatives. 

I 
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Design  Constraints  - Any  parameter  which  must 
not  exceed  specified  bounds  for  the  design  to  be 
acceptable.  Design  constraints  may  be  linear 
or  non  linear,  implicit  or  explicit,  but  they 
must  be  functions  of  the  design  variables. 

Design  constraints  appear  only  on  the  left  side 
of  equations. 

Objective  Function  - The  parameter  which  is  going 
to  be  minimized  or  maximized  during  the  optimization 
process.  The  objective  function  may  also  be  linear 
or  non  linear,  implicit  or  explicit,  and  must  be  a 
function  of  the  design  variables.  The  objective 
function  usually  appears  on  the  left  side  of  an 
equation.  The  only  exception  is  if  the  objective 
function  is  also  a design  variable. 

As  can  readily  be  seen  by  the  definitions  above,  design 
constraints  and  objective  functions  are  usually 
interchangeable . 

The  number  of  design  variables  being  utilized  for  an 
optimization  is  equivalent  to  the  dimension  of  the  design 
space  in  which  the  design  is  being  calculated.  Thus,  if 


an  optimization  problem  has  four  design  variables  specified, 
the  design  will  be  carried  out  in  four-space. 

Assuming  that  the  optimization  process  requires  the 
minimization  of  a particular  objective  function,  the  general 
optimization  problem  can  be  stated  as; 

Find  the  vector  of  design  variables,  X,  to 
minimize  FCX)  subject  to  the  constraints 
Gj  (I)  < 0.0  , j « 1 ,\CON 


VLB.  1 X < VUB.  , i =■  1,NDV 


(1) 

(2) 


In  the  general  problem  statement,  FCX)  is  the  objective 
function,  there  are  NDV  design  variables,  and  NCON  constraints. 
VLB^  and  VUB^  are  the  lower  bounds  and  upper  bounds  respec- 
tively on  the  i-th  design  variable.  If  the  inequality 
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U condition  of  equation  (1)  is  violated,  (G^ (X)  > 0) , for  any 

constraint,  that  constraint  is  said  to  be  violated.  If  the 
equality  condition  is  met,  (G^ (X)  = 0),  the  constraint  is 
active.  If  the  inequality  condition  is  met,  (G^.  (X)  < 0), 
the  constraint  is  inactive.  Because  of  the  numerical 


problems  involved  in  representing  exact  zero  on  a computer 
with  a finite  number  of  significant  figures,  the  equality 


condition  is  represented  by  a band  around  the  value  Gj  C^)  = O' 
The  band  is  equal  to  twice  the  constraint  thickness  (CT)  and 


is  shown  in  Figure  1.  Figure  1 is  a representation  of  a 


two-variable  design  space  with  the  values  of  G^ (X)  = 0 
plotted.  In  this  instance. 


For  n-space,  G^ (X)  would  appear  as  a hypersurface 
whereas  for  n = 2,  the  hypersurface  degenerates  to  a single 
curve  that  can  be  easily  represented. 

Any  design  which  satisfies  the  inequalities  of 
equations  (1)  and  (2)  is  referred  to  as  a feasible  design. 
If  a design  violates  any  of  these  constraints,  it  is  an 
infeasible  design.  The  minimum  feasible  design  is  said  to 
be  optimal.  Note  that  if  it  is  desired  to  maximize  an 
objective  function,  the  process  reduces  to  minimizing  the 
negative  of  the  objective  function. 

CONMIN  requires  an  initial  vector  of  design  variables, 
7,  which  may  or  may  not  yield  a feasible  design.  If  the 
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initial  design  is  infeasible,  CONMIN  moves  towards  a 
feasible  solution  with  a minimal  increase  in  the  objective 
function  [9].  The  optimization  process  then  proceeds  in 
an  iterative  fashion  with  the  following  recursive 
relationship : 


+ a*S^ 


where  q is  the  iteration  number  and  a*  is  the  move  parameter, 
a scalar,  which  defines  the  distance  of  travel  in  the 
direction  of  search,  S. 

The  optimization  process  is  divided  into  two  parts. 

The  first  is  the  determination  of  J which  will  reduce  the 
objective  function  without  violating  any  constraints.  The 
second  is  the  determination  of  the  scalar  ct*  so  that  F(X) 
is  minimized  in  this  direction,  a new  constraint  is 
encountered,  or  a currently  active  constraint  is  encountered 
again . 

For  the  sake  of  discussion,  consider  a condenser  design 
problem  with  two  design  variables,  and  X,,  where 

X^  = condenser  tube  outside  diameter,  and 


^2  = tube  pitch  to  diameter  ratio. 


The  objective  function  is  condenser  volume,  VOL(X).  .‘\ssume 
that  the  tube  bundle  diameter  must  be  greater  than  a 
specified  value,  and  that  the  cooling  water  pumping 

power  must  be  less  than  a specified  value,  Figure  2 
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is  a graphical  representation  of  the  problem.  Note  the 
lines  of  constant  objective  function,  with  VOL^ (X)  > VOL2 (X) , 
and  the  initial  design  at  point  (a) . 

It  must  be  reiterated  that,  while  this  example  assumes 
a feasible  initial  design,  this  is  not  a requirement  and 
CONMIN  will  also  optimize  when  given  an  infeasible  initial 
design. 

The  optimization  begins  by  calculating  the  gradient  of 
the  objective  function  by  using  the  finite  difference  tech- 
nique. A perturbation  of  0.01  is  applied  to  each  of  the 
design  variables  in  a single  forward  step.  The  gradient 
vector  is  therefore 


’a  cvoif 

1 

“a  VOL 

3Xi 

AXi 

7F(X)  = VVOL  = 

9 CVOL) 

A VOL 

ax. 

» J 

Because  no  constraints  are  active  or  violated,  the  greatest 
improvement  in  the  objective  function  will  be  realized  by 
moving  in  the  direction  of  steepest  descent  so  that 

S * - WOL 

as  shown  in  Figure  3. 

With  the  value  of  S’  now  determined,  it  remains  to  find 
the  move  parameter,  a*,  that  will  allow  the  greatest 
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improvement  in  the  objective  function.  A one-dimensional 
search  is  carried  out  in  the  S’  direction  until  the  value  for 


a*  is  found.  This  is  point  @ on  Figure  3.  The  location 
of  point  (b)  terminates  the  first  design  iteration. 

The  second  design  iteration  is  begun  by,  once  again, 
perturbing  X to  find  WOL.  Instead  of  moving  in  the 
direction  of  steepest  descent,  however,  a new  S is  found 
by  the  method  of  conjugate  directions,  developed  by 
Fletcher  and  Reeves  [11].  With  this  method,  S is  calculated 
by  the  following  relationship; 


This  is  shown  in  Figure  4. 

The  advantage  of  using  the  Fletcher  - Reeves  method 
instead  of  the  steepest  descent  is  that  convergence  to  an 
optimum  is  much  faster.  With  the  new  search  direction,  5^, 
a search  is  performed  in  this  direction  until  a constraint 
is  encountered.  This  occurs  at  point  on  Figure  4 at 
the  pumping  power  constraint,  thus  terminating  the  second 
design  iteration. 

The  third  design  iteration  begins  with  the 
constraint  active  at  point  @ on  Figure  5.  As  7V0L  is 
found,  the  gradient  of  the  active  constraint  is  found  using 
the  information  from  the  same  forward  finite  difference  step. 

Not  only  is  a new  required  that  will  reduce  the 
objective  function,  but  this  new  S must  not  violate  the 
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active  constraint.  This  problem  may  be  solved  by  the  method 
of  feasible  directions  developed  by  Zoutendijk  [12]  and 
implemented  by  Vanderplaats  and  Moses  [13]. 

The  finding  of  a new  S has  now  become  a sub-problem 
which  is  a linear  programming  problem  with  a single 
quadratic  constraint.  This  sub-problem  [14]  is  stated  as: 

Find  a vector  S to  maximize  B subject  to  the  constraints 


7 F(X)  . S’  + B < 0 (3) 

V Gj  (X]  • S’  + 9^3  < 0 j = l.NAG  (4) 

S’  • S’  < 1 (5) 

where,  for  the  case  at  hand,  7 F(X)  =7V0L,^Gj(X)  = -7hP, 
and  Gj  (X)  = - HP.  NAG,  the  number  of  active  constraints, 

is  one  in  this  instance. 


If  equation  (3)  is  satisfied  and  3 is  positive,  the 
resulting  search  direction  will  reduce  the  objective  function 
and  is  defined  as  a usable  direction.  If  equation  (4)  is 
satisfied  and  3 is  positive,  S is  a feasible  direction 
because,  for  a small  move  in  this  direction,  no  constraints 
will  be  violated.  9^  is  defined  as  the  push-off  factor  for 
Gj (T)  and  has  the  effect  of  pushing  the  design  away  from 
the  active  constraint.  9^  must  be  zero  or  positive  to 
maintain  a feasible  design.  For  9^  * 0,  ? would  be  tangent 
to  the  active  constraint.  For  9j  >>  1.0,  S'  would  be 
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pushed  away  from  the  active  constraint  and  very  nearly 
tangent  to  a line  of  constant  objective  function.  For 
9j  * 1.0,  the  angle  between  constant  objective  function  and 
the  active  constraint  would  be  approximately  bisected.  If 
the  maximum  value  of  6 obtainable  from  equations  (3) , (4) 
and  (5)  is  zero,  then  no  direction  exists  which  will  both 
reduce  the  objective  function  and  satisfy  the  constraint, 
and  the  current  design  is  optimal  or  is  at  least  a local 
minimum.  In  this  example,  a usable  and  feasible  direction 
exists  and  a one  dimensional  search  leads  to  point  in 

Figure  5 where  the  minimum  bundle  diameter,  BD  ■ , constraint 

is  encountered.  This  ends  the  third  iteration  in  the 
optimization  process. 

The  fourth  iteration  begins  as  before,  with  the  calcula- 
tion of  the  gradient  of  the  objective  function  and  the 
active  constraint.  The  sub-problem  of  equations  (3)  through 
C5)  is  again  solved  for  a new  S. 

It  should  be  noted  that  the  minimum  bundle  diameter 

constraint  is  assumed  to  be  linear;  therefore  in  equation  (4) 

the  push-off  factor,  0^,  is  set  to  zero  allowing  to 

follow  the  constraint  as  shown  in  Figure  6.  A one 

dimensional  search  is  again  carried  out  in  the  new  S’  direction 

until  a new  constraint  is  encountered  or  an  active  constraint 

is  re-encountered.  Thus,  the  activated  BD  . constraint  is 

min 

"ridden"  until  no  further  design  improvement  is  realized. 


This  occurs  at  point  (e)  on  Figure  6 and  the  fourth  design 
iteration  is  terminated. 


For  the  fifth  iteration,  the  same  procedure  is  followed 
and  the  solution  to  the  sub-problem  characterized  by 
equations  (3)  through  (5)  yields  a value  of  zero  for  3. 

Thus,  there  is  no  direction  that  will  both  reduce  the 
objective  function  and  satisfy  the  constraint  and  the 
current  design  at  point  (e)  in  Figure  7 is  optimal. 

Figure  8 illustrates  the  value  of  using  optimization 
techniques  to  solve  the  design  problem.  Assume  an  initial 
design  at  point  @ such  that  the  minimum  tube  outside 
diameter  is  active.  A reduction  in  the  tube  pitch  to 
diameter  ratio  will  yield  a design  improvement  until  the 
minimum  bundle  diameter  constraint  is  encountered  at  point 
(D  . At  this  point,  neither  the  pitch  to  diameter  ratio 
nor  the  outside  tube  diameter  can  be  reduced  independently 
as  required  to  reduce  the  objective  function,  without 
violating  the  active  constraints.  Only  by  changing  the 
two  design  variables  in  a particular  manner  can  the  minimum 
value  of  the  objective  function  at  point  © be  achieved. 

This  discussion  of  the  methodology  involved  with  CONMIN 
would  not  be  complete  without  citing  the  program's  limita- 
tions. The  number  of  design  variables  (NDV)  directly 
affects  the  computational  time  required  to  reach  an  optimum. 

Since  the  calculation  of  the  gradient  information  required 

for  each  design  variable  at  the  beginning  of  each  design 

iteration  is  found  by  using  a single  forward  finite  difference  4 

step,  requiring  a complete  pass  through  the  analysis  portion  | 

of  the  program,  there  is  a subsequent  increase  in  CPU  time 
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as  NDV  increases.  Also,  problems  with  many  design  variables 
tend  to  converge  more  slowly  due  to  the  interaction  between 
the  design  variables  and  because  of  the  numerical  inaccuracy 
(machine  related)  generated  during  the  optimization  process. 
The  number  of  design  variables  should  therefore  be  kept 
small  in  order  to  expedite  the  optimization  process. 
Vanderplaats  [7]  recommends  a practical  limit  of  twenty 
design  variables. 

The  number  of  design  constraints  used  is  not  limited  in 
the  same  manner.  Recall  that  the  only  time  gradient 
information  is  stored  for  a constraint  is  when  that  constraint 
is  active.  Therefore,  the  sheer  number  of  constraints  will 
not  adversely  affect  the  optimization. 

As  can  be  seen  in  Figure  9,  it  is  quite  possible  that 
the  optimal  design  found  is  actually  a relative  optimum  and 
not  a global  optimum.  This  problem  can  be  overcome  by 
starting  the  design  with  several  different  initial  vectors, 

X,  until  the  same  optimal  design  is  repeated. 

Equality  constraints  of  the  type 

= 0 

are  very  difficult  to  provide  for  in  the  general  optimization 
scheme  [7].  By  defining  an  objective  function  in  which  a 
weighting  factor  multiplies  the  parameter  to  be  held  constant, 

Y,  and  this  product  is  summed  with  the  parameter  to  be 
optimized,  X, 

objective  function  ■ X + weighting  factor  x Y 
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the  parameter  to  be  held  constant  can  be  forced  to  the 
appropriate  bound.  The  weighting  factor  is  generally  one 
that  will  keep  both  terms  in  the  same  order  of  magnitude 
but  trial  and  error  is  sometimes  required.  Parameters  can 
be  held  "constant"  within  tO.5  percent. 

C.  CONTROL  PROGRAM  FOR  ENGINEERING  SYNTHESIS  (COPES) 

Recall  that  the  optimization  program,  CONMIN,  was 
written  in  subroutine  form.  Vanderplaats  [15]  has  developed 
a main  program  to  simplify  the  use  of  CO.NMIN  and  to  further 
aid  in  the  design  optimization  process. 

The  user  must  supply  an  analysis  subroutine  with  the 
name  A.NALIZ.  .ANALIZ,  in  keeping  with  general  good  program- 
ming practice,  must  have  three  segments:  input,  analysis  and 
output.  Based  on  the  value  of  a flag  from  COPES 
(ICALC  * 1,2  or  3),  .ANALIZ  performs  the  proper  function. 

Finally,  CONMI.N  and  A.NALIZ  do  not  communicate  directly  with 
each  other  as  COPES  is  the  main  program. 

COPES  is  constantly  being  revised  by  Vanderplaats  to 
better  meet  the  needs  of  the  engineer. 

The  COPES  program  currently  provides  four  specific 
capabilities  : 

1.  Single  analysis  - just  as  if  COPES/CONMIN  were 
not  used. 

2.  Optimization  - minimization  or  maximization  of  a 

multivariable  function  with  limits  imposed  on  other  I 

I 

functions. 
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3.  Sensitivity  analysis  - used  to  investigate  the 
effect  of  changing  one  or  more  design  variables 
on  one  or  more  calculated  functions. 

4.  Two-variable  function  space  - provides  analyses 
of  all  specified  combinations  of  two  design 
variables . 

This  work  is  concerned  only  with  the  application  of 
COPES/CONMIN  to  the  optimized  design  of  a marine  condenser, 
therefore,  items  3 and  4 are  included  only  for  completeness. 

If  there  is  no  relative  or  absolute  change  in  the  value 
of  the  objective  function  for  three  design  iterations, 
the  optimum  is  found  and  COPES  prints  the  appropriate 
message  and  terminates  the  program.  If  no  feasible  design 
can  be  found  after  ten  design  iterations,  COPES  prints  the 
appropriate  message  and  terminates  the  program. 

Experience  has  shown  that  most  design  problems  can  be 
optimized  within  20  design  iterations  and  the  maximum  number 
of  design  iterations  permitted  is  defaulted  to  this  value. 
Thus,  if  the  optimum  value  has  not  occurred  in  20  design 
iterations,  COPES  will  print  the  appropriate  message  and 
terminate  the  procedure. 

COPES  has  simplified  the  procedures  involved  in  using  a 
sophisticated  program  such  as  CONMIN.  Thus  the  engineer  is 
freed  from  the  unwanted  role  of  systems  analyst  and  may 
devote  his  talents  to  engineering. 


III.  OPTIMIZED  CONDENSER  DESIGN,  VERSION  1 (OPCODEl) 


A.  DEVELOPMENT 

The  Bureau  of  Ships  adopted  the  Heat  Exchange  Institute's 

Standards  for  Steam  Surface  Condensers  [1]  by  issuing  Design 

Data  Sheet  DPS  4601-1  [2]  in  October,  1953.  The  DDS  is  still 

being  used  to  specify  naval  condensers. 

The  technique  involved  is  based  on  calculating  a value 

of  the  overall  heat  transfer  coefficient,  LF  , based  on  the 

c 

cooling  water  velocity  through  the  condenser  tubes,  the 
condenser  tube  wall  thickness  and  material,  the  tube  fouling 
factor,  and  the  cooling  water  inlet  (injection)  temperature. 
Knowledge  of  these  parameters  and  their  associated  correction 
factors  leads  to  the  simple  formulation  of 


U. 


F^FjFj  C /y 


The  correction  factors  F, , F,,  and  F,  are  tabulated  in 

12  j 

references  [1]  and  [2].  Reference  [2]  specifies  a value  of 
C * 270  for  0.625  and  0.750  inch  outside  diameter  (o.d.), 

13  Birmingham  Wire  Gauge  (BWG)  tubes.  A value  of  C = 263  is 
used  for  0.875  and  1.00  inch  o.d.,  18  BWG  tubes. 

Because  of  the  simplicity  of  the  HEI/DDS  method  of 
condenser  design,  and  since  the  DDS  is  still  the  specifying 
document  for  naval  condensers,  the  author  chose  to  implement 
the  HEI/DDS  method  for  automated  design.  The  combination  of 
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COPES/CONMIN  with  the  HEI/DDS  method  was  named  OPCODEl 
(optimized  COndenser  Design,  Version  1_)  . The  development  of 
OPCODEl  is  presented  in  Appendices  A and  B. 

To  add  greater  versatility  to  the  HEI/DDS  method, 
algorithms  were  developed  for  OPCODEl  to  calculate  a 
circular  bundle  geometry  and  the  cooling  water  pumping 
requirement.  The  circular  bundle  geometry  was  designed  from 
the  long  axis  out.  A central  12  inch  diameter  void  to  serve 
as  the  collection  header  for  noncondensable  gasses  was 
provided  along  the  bundle's  longitudinal  axis.  Once  the 
number  of  tubes  was  calculated,  the  tubes  were  placed  in 
rows  concentric  to  the  central  void  and  spaced  with  a 60 
degree  triangular  pattern.  To  simplify  the  algorithm,  and  to 
provide  continuous  functions  to  COPES/CONMI.V , partial  tubes 
were  permitted  in  a row  and  the  outermost  row  was  allowed 
to  be  partially  filled. 

The  calculation  of  the  mass  flow  rate  of  cooling  water 
required  began  the  calculation  of  the  sea  water  (S-W)  pressure 
drop  and  pumping  power  requirement. 

The  S-W  pressure  drop  was  divided  into  a component  based 
on  the  sudden  expansion  and  contraction  losses  caused  by  the 
flow  entering  and  exiting  the  waterboxes;  a component  based 
on  the  sudden  contraction  and  expansion  losses  caused  by 
flow  from  the  inlet  waterbox  into  the  inlet  tube  sheet  and 
from  the  outlet  tube  sheet  into  the  outlet  waterbox;  and  a 
component  based  on  the  normal  frictional  losses  associated 
with  internal  tube  flow. 
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Reference  1 presents,  in  graphical  form,  the  pressure 
drops  for  all  the  components.  The  graphical  data  from 
Figure  SF-9  of  reference  1 for  the  calculation  of  waterbox 
losses  was  implemented  using  a systems  library  interpolating 
subroutine.  The  S-W  velocity  in  the  waterboxes  was  taken 
as  75  percent  of  the  velocity  in  the  condenser  tubes  as 
specified  in  reference  16. 

To  simplify  the  tube-side  pressure  loss  algorithm,  the 
friction  factor,  f,  was  first  calculated  by  solving  the 
transcendental  Colebrook  equation  [17].  Details  of  the 
solution  are  provided  in  Appendix  B.  An  absolute  roughness 
of  5 x 10  ^ feet  was  assumed  as  representative  of  drawn 
tubing.  With  the  value  of  f calculated,  the  tube-side 
pressure  drop,  was  calculated  using  the  familiar 

Darcy-Weisbach  formula  [17]: 

L 

"Pt  = ^ . 

The  tube  sheet  entrance  and  exit  losses  were  calculated 
by  utilizing  the  area  ratio  technique  developed  in 
reference  18.  With  this  method,  the  flow  area  associated 
with  each  tube  was  calculated 

A » tube  sheet  area 

"F  number  of  tubes 

followed  by  the  calculation  of  the  internal  flow  area  for 
an  individual  tube,  A£  and  the  ratio  of  the  two  areas: 
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With  the  value  of  S calculated,  the  entrance  loss  was 
calculated  with  the  relation 


= 0 . 5 ( 1 . 0 - 3 ^)  ( [feet  of  water  column] 

and  the  exit  loss  was  calculated  with  the  relation 


^2  2 

AP  = (1.0-  3“]  [feet  of  water  column]. 

*0 


Although  the  heat  transfer  calculations  performed  in 
OPCODEl  were  for  a single  pass  condenser,  expansion  of  the 
program  to  multiple  pass  condensers  is  feasible.  To 
simplify  this  future  expansion,  the  waterbox  pressure  loss 
was  multiplied  by  the  number  of  tube  passes,  read  as  an 
input  variable. 

The  pumping  power  was  now  easily  calculated 


PP  = N • V • A.  • p • P 
f sw 


rft-  Ibf, 

'■  sec  ' 


The  input  variables  required  for  the  use  of  OPCODEl 
are  listed  in  the  User's  Manual  provided  as  Appendix  C. 


B.  OPCODEl  VERIFICATION 

It  was  desirable  to  verify  the  performance  of  OPCODEl 
as  a predictor  of  condenser  performance  by  comparison  with 
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actual  experimental  data.  Complete  and  accurate  data  was 
not  available,  therefore,  another  method  of  model  verification 
was  sought. 

Using  the  design  data  from  the  condenser  technical 
manuals  for  two  classes  of  ships,  the  CVA  67  and  the  DE  1040, 
an  attempt  was  made  to  repeat  the  design  of  the  condensers. 

The  two  condensers  are  vastly  different  in  size;  the  CVA  67 
has  a heat  transfer  area  of  16,011  square  feet  while  the 
DE  1040  class  condenser  has  a heat  transfer  area  of 
6600  square  feet. 

The  results  of  the  design  of  the  CVA  6"'  by  OPCODEl 
are  tabulated  in  Table  I with  the  data  from  the  condenser 
technical  manual  [19]  included  for  comparison.  Very  close 
agreement  was  achieved  for  all  parameters  except  tube-side 
pressure  drop.  This  difference  was  attributed  to  the  tube 
sheet  layout  by  OPCODEl  and  the  fact  that  no  steam  lanes 
were  accounted  for.  The  addition  of  steam  lanes  would  tend 
to  increase  the  tube-side  pressure  drop  since  more  tube  sheet 
area  would  be  unused  as  tube  sites.  Thus  the  entrance  and 
exit  losses  associated  with  the  tube  sheets  would  be  greater 
than  predicted  by  OPCODEl. 

The  results  of  the  design  of  the  DE  1040  class  condenser 
by  OPCODEl  are  tabulated  in  Table  II  with  the  data  from  the 
condenser  technical  manual  [20]  included  for  comparison. 
Excellent  correlation  was  achieved  with  all  parameters  within 
three  percent  of  the  specifications  except  for  bundle 
diameter  and  tube-side  pressure  drop.  Since  no  provision  was 
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made  for  the  calculation  of  circumferential  steam  lanes 
in  OPCODEl,  the  bundle  diameter  calculated  was  less  than 
the  prototype. 

Since  the  DE  1040  condenser  has  a basically  circular 
bundle  geometry,  close  agreement  was  expected  between  the 
calculated  tube-side  pressure  drop  and  the  value  specified 
in  the  technical  manual.  This  was  not  the  case.  One 
plausible  explanation  could  be  that  the  designers  specified 
a large  factor  of  safety  for  this  parameter. 

These  results  confirm  that  the  designers  of  the  CVA  b"’ 
and  thn  DE  1040  main  condensers  did,  in  fact,  use  the 
HEI/DDS  method  for  the  calculation  of  the  heat  transfer 
parameters  and,  with  the  notable  exceptions  of  tube-side 
pressure  drop  and  bundle  diameter,  OPCODEl  accurately 
predicts  these  parameters  at  the  full  power  design  point. 

If  the  limitations  imposed  on  tube-side  pressure  drop 
and  bundle  diameter  by  the  exclusion  of  steam  lanes  from 
OPCODEl  are  acknowledged,  the  results  received  from  OPCODEl 
can  be  viewed  as  a "first  cut”  analysis.  This  was  the 
original  purpose  of  developing  this  program. 

C.  LIMITATIONS  OF  OPCODEl 

OPCODEl  is  insensitive  to  shell-side  conditions.  The 
saturation  steam  pressure  and  the  saturation  steam  temperature 
are  assumed  to  remain  constant  as  the  steam  passes  through 
the  bundle,  whereas  the  steam  flow  actually  experiences  a 
pressure  drop  with  a resulting  decrease  in  saturation 
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temperature  as  the  steam  passes  through  each  row  of  tubes. 
.As  a result,  OPCODEl  makes  no  attempt  to  specify  steam 
lanes,  either  circumferential  or  radial,  since  the  proper 
design  of  steam  lanes  is  a strong  function  of  local  steam 
pressure . 

OPCODEl  has  no  provision  for  the  application  of  tube 
enhancement  factors  due  to  the  simple  manner  in  which  the 
corrected  overall  heat  transfer  coefficient  is  computed. 
The  usual  method  for  the  calculation  of  [21]  is: 


U. 
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.A. 
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In  (r^/b^) 
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R ^ ^ 
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where  the  internal  film  heat  transfer  coefficient  is  given 
by  the  Colburn  form  of  the  convective  film  heat  transfer 
correlation  [21]: 


h. 

1 


0.023  (Ref  • ^ 


(7] 


The  outside  film  heat  transfer  coefficient  is  given  by 
Nusselt's  equation  [21]  corrected  for  the  inundation,  or 
condensate  rain,  effects  of  upper  tubes  on  lower  tubes  in 
the  bundle  by  including  the  factor  n*  in  the  denominator, 
where  n*  is  the  number  of  tubes  in  a vertical  row  above  the 
i-th  tube 
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With  and  in  the  forms  given  by  equations  (7)  and 
(8),  enhancement  factors  can  be  applied  as  multipliers  of 
these  equations  and  equation  (6)  will  yield  a value  for 
the  enhanced  overall  heat  transfer  coefficient, 

None  of  the  foregoing  is  possible  with  the  HEI/DDS 
based  OPCODEl. 

Finally,  since  OPCODEl  is  insensitive  to  steam  conditions, 
the  effect  of  noncondensable  gasses  on  the  heat  transfer 
characteristics  of  the  condenser  are  unknown. 
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IV.  OPTIMIZED  CONDENSER  DESIGN,  VERSION  2 (OPCODE2) 


A.  BACKGROUND 

In  the  late  1960's,  engineers  at  the  Oak  Ridge  National 
Laboratory  developed  a sophisticated  computer  code  under 
contract  to  the  Office  of  Saline  Water.  This  code,  called 
ORCO.Nl  [3],  was  generated  to  aid  in  the  analysis  and 
parametric  study  of  large,  generally  circular  steam 
condensers  for  use  in  large  scale,  multistage  distillation 
plants  for  the  production  of  potable  water  from  sea  water 
by  the  flash  evaporation  process. 

.Much  of  ORCONl  was  dependent  on  Eissenberg's  research 
work  [4]  on  the  effects  of  condensate  rain  on  the  shell- 
side  convective  heat  transfer  coefficient. 

The  program  analyzes  a single  pass,  circular  or  semi- 
circular condenser,  with  steam  flowing  on  the  shell-side 
of  the  tubes  and  variable  salinity  water  flowing  on  the 
tube-side.  .An  optional,  rectangular  air  cooler  bundle  is 
provided  for,  as  well  as  elementary,  shell-side  baffles. 

The  bundle  is  divided  into  30  degree  sectors  and  symmetry 
about  the  central  axis  may  be  employed  to  reduce  computa- 
tional effort.  The  tubes  are  placed  on  a 60  degree  equilat- 
eral triangular  pattern  of  concentric  rows  with  the  rows 
added  from  the  outermost  row  to  an  inner  void  provided 
along  the  bundle's  longitudinal  axis.  This  serves  as  a 
collection  header  for  noncondensable  gasses  prior  to  passage 
through  the  air  cooler,  if  specified.  The  steam  is  assumed 
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outside  of  the  bundle  to  the 


to  flow  radially  from  the 
central  void.  Figure  10  shows  the  ORCONl  model  of  a steam 
condenser  with  optional  baffles  and  cooler. 

ORCONl  proceeds  with  sector  by  sector,  row  by  row 
calculations  of  the  following  quantities: 

a)  Saturation  temperature  of  the  steam  entering  the 
row. 

b)  Pressure  of  the  steam  plus  noncondensable  gas 
entering  the  row. 

c)  Steam  flow  entering  the  row. 

d)  Steam  and  noncondensable  gas  velocity  at  the  minimum 
cross  section  in  the  row. 

e)  The  fraction  of  noncondensable  gas  by  weight. 

f)  The  overall  heat  transfer  coefficient  for  the 
average  tube  in  the  row. 

g)  The  steam-side  condensing  coefficient. 

h)  The  tube-side  heat  transfer  coefficient. 

i)  The  shell-side  film  heat  transfer  coefficient 
composed  of  the  noncondensable  gas  filh.  . lus  the 
condensate  film. 

j)  The  shell-side  Reynolds  number  based  on  the  mass 
flow  at  the  minimum  cross  sectional  area  in  the  row. 

k)  The  heat  flu.x  per  square  foot  of  condenser  tube. 

l)  The  shell-side  friction  factor. 

m)  The  mass  flow  rate  of  steam  plus  noncondensable  gas 
at  the  minimum  cross  section  in  the  row. 


n) 
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the 
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q)  The  number  of  tubes  per  row. 


r)  The  cumulative  shell-side  pressure  drop  from  row  1 
to  row  n. 

In  addition  to  the  above  parameters,  the  area  weighted 
overall  heat  transfer  coefficients  for  the  condenser  section, 
the  cooler  section,  and  the  combined  condenser  are  used  to 
calculate  the  "back  calculated"  log  mean  temperature 
difference  (LMTD) . This  value  of  the  LMTD , when  compared 
with  the  LMTD  calculated  by  standard  means  using  the 
saturation  steam  temperature  and  the  average  cooling  water 
inlet  and  exit  temperatures,  represents  the  loss  in  average 
thermal  driving  force  due  to  pressure  drops  within  the 
bundle . 

The  tube  internal  film  heat  transfer  coefficient  is 
calculated  by  the  Colburn  equation  multiplied  by  the 
internal  enhancement  factor 

h.  = 0.023  (Re)°-®  (Pr)^/^  (^ . ) . 

The  uncorrected  external  film  heat  transfer  coefficient 
is  calculated  using  the  basic  Nusselt  equation  [21]: 

1/4 

(E^)  . (9) 


0.725 


V 3 2 , 

0^  h g, 


The  value  of  h^  calculated  in  equation  (9)  is  for  a 
single  tube.  Eissenberg  [4]  corrects  for  inundation  effects 
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by  first  calculating  a tube  flooding  factor  using  the 
following  relation: 

F = 0.6F,  + (1  - 0. 5647F  ,)n*^- . 

n d d 

Fj  is  an  input  parameter  and  is  a function  of  both  tube 
spacing  to  diameter  ratio  and  tube  orientation  [3,22]. 

The  condensate  film  coefficient  for  the  typical  tube  in 
the  n-th  row  is  then  calculated  by  correcting  the  value  of 
h^  from  equation  (9) ; 

h^*(n)  = [nF„  - (n  - 1)F„  , ]h^ 

0 n n-io 

To  model  the  main  condenser  of  the  CVA  67,  one  quarter 
of  the  full  circular  bundle  was  specified.  With  two  bundle 
quarters  placed  back  to  back  as  shown  in  Figure  11,  the 
steam  lanes  and  tube  arrangement  closely  approximated  the 
CVA  67  tube  sheet  as  shown  in  Figure  12  from  reference  19. 

The  combination  of  the  COPES/CO.NMIN  optimization  package 
with  a suitably  modified  ORCONl  produced  0PC0DE2  (OPt imi zed 
condenser  DEsign , Version  _2 ) . 

B.  MODIFICATIONS  TO  ORCONl 

The  original  ORCONl  had  neither  tube-side  pressure  drop 
calculations  nor  volumetric  calculations.  The  subroutines 
developed  for  OPCODEl  to  calculate  the  tube-side  pressure 
drop  and  pumping  power  were  installed  in  ORCONl  when 
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converting  it  to  0PC0DE2.  An  equation  for  the  calculation 
of  condenser  volume  based  strictly  on  the  model  as  developed 
for  the  CVA  67  main  condenser  and  shown  in  Figure  10  was 
added  to  ORCONl. 

The  flowchart  given  in  Figure  13  illustrates  the 
original  program  logic  for  ORCONl.  Once  all  calculations 
were  completed,  a test  for  the  exit  steam  fraction  was  made 
against  an  input  target  value.  If  the  target  value  was  not 
met,  and  depending  on  the  value  of  an  input  flag,  either 
the  length  of  the  condenser  tubes  or  the  quantity  of  inlet 
steam  was  varied,  and  the  calculations  were  repeated  until 
the  target  value  and  the  calculated  values  agreed.  The 
adjustment  of  tube  length  or  inlet  steam  was  performed  in 
subroutine  .ADJUST. 

It  was  felt  that  the  quantity  of  inlet  steam  was  generally 
the  value  which  should  drive  the  entire  design  of  the 
condenser  and  should  therefore  remain  a constant. 

Similarly,  the  tube  length  was  a good  candidate  for  inclusion 
in  the  optimization  program  as  a design  variable.  For  these 
reasons,  subroutine  ADJUST  was  removed  from  ORCONl. 

The  original  ORCONl  provided  logic  for  multiple  data 
runs.  This  capability  was  removed  since  COPES/CONMIN 
requires  a "once  through”  analysis  program.  Figure  14  shows 
the  sections  of  ORCONl  removed  during  the  conversion  to 
OPCODE2. 

Search  [5]  had  previously  modified  ORCONl  so  as  to 
calculate  the  shell-side  heat  transfer  coefficient  either 
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with  the  original  equation  from  the  work  of  Eissenberg  [4] 
or  by  using  the  relation  specified  by  the  manufacturer  of 
Korodense  condenser  tubes  [23].  This  facility  was  left  in 
0PC0DE2 . 

ORCONl  was  not  originally  a machine  independent  program 
as  it  had  been  tailored  for  operation  on  IBM  (International 
Business  Machines)  equipment.  This  was  not  a desirable 
characteristic  and  programming  changes  were  made  to  ensure 
that  ORCONl  was  machine  independent  and  that  0PC0DE2  would 
remain  machine  independent. 

ORCONl  uses  iterative  techniques  to  solve  for  quantities 
such  as  condensation  rate,  steam  mass  flow  rate,  and  pressure 
drop  balance  between  sectors.  The  description  of  subroutine 
ADJUST  is  an  e.xcellent  e.xample  of  such  an  application. 

CONMIN  uses  perturbation  techniques  to  calculate  the 
gradient  information  required  for  each  design  variable  and 
for  each  active  design  constraint  during  an  optimization 
iteration. 

Since  ORCONl  has  the  capability  to  make  design  decisions, 
a perturbation  by  CONMIN  would  cause  an  adjustment  by  ORCONl. 
The  two  programs  therefore  worked  at  cross  purposes. 

It  was  not  possible  to  remove  all  of  these  design 
decision  points  from  ORCONl  because  of  the  strong  coupling 
between  the  subroutines.  The  decision  was  made  to  remove 
only  ADJUST  and  to  leave  the  other  design  decision  points 
intact . 
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This  problem  affected  the  choice  of  parameters  that 
could  be  used  as  design  variables  and  the  range  of  values 
that  the  chosen  design  variable  was  permitted  to  assume. 

Another  problem  area  developed  because  many  of  the  thermo- 
dynamic properties  were  calculated  in  ORCONl  by  subroutines 
that  use  logarithmic  functions  to  approximate  the  thermo- 
dynamic curves.  If  the  arguments  for  these  functions  take 
on  values  less  than  or  equal  to  :ero,  the  function  is 
undefined.  While  constraints  could  be  added  as  part  of  the 
optimization  process,  the  computation  would  stop  during  the 
pass  through  the  analysis. 

Tests  were  put  before  all  calculations  that  involved 
logarithmic  evaluation  and  before  all  function  evaluations 
where  the  denominator  could  take  on  values  very  close  to 
zero.  These  tests,  when  activated,  would  make  a "design 
decision"  by  setting  the  function  to  an  approximate  value 
such  as : 

if  X £ 0.0001  , then  y = -10.0 
otherwise,  y = In  x . 

These  approximate  values  would  cascade  due  to  the  large 
number  of  times  the  function  was  evaluated.  Other  mathematical 
instabilities  also  occurred. 

The  design  variables  or  the  values  of  the  design  variables 
that  would  trigger  the  original  instability  were  found.  The 
tests  were  removed  and  the  troublesome  design  variables  and/or 
the  particular  values  that  would  cause  the  undesirable 
response  were  avoided. 


47 


c.  opcode:  verification 

Verification  of  OPCODE2  as  a predictor  of  condenser 
performance  was  attempted  by  inputing  the  design  values 
from  the  CVA  67  technical  manual  [19]  and  comparing  the 
condenser  designed  by  0PC0DE2  with  the  values  given  by  the 
condenser  technical  manual  [19].  Both  the  Eissenberg  [3] 
and  the  Korodense  [23]  relations  for  tube  inundation 
effects  were  used. 

Because  only  one  half  of  the  condenser  is  designed  by 
ORCO.Nl , the  quantity  of  inlet  steam  and  the  number  of  tubes 
specified  by  reference  19  were  halved. 

The  program  would  not  accept  an  inlet  steam  rate  of 
greater  than  200,000  pounds  per  hour  (400,000  pounds  per  hour 
for  the  entire  bundle). 

Holman  [21]  states  that  up  to  a 20  percent  increase  in 
heat  transfer  rate  may  be  realized  by  the  ripples  set  up  in 
the  condensate  film  by  steam  passing  over  the  film.  The 
Korodense  literature  [23]  indicates  an  enhancement  due  to 
the  same  phenomena  as  less  than  the  20  percent  reported  by 
Holman  but  also  shows  the  enhancement  to  be  a weak  function 
of  tube  location  in  the  bundle. 

An  enhancement,  due  to  film  ripples,  of  14  percent  was 
assumed  for  the  current  work. 

The  results  of  the  verifications  using  both  the  Eissenberg 
and  the  Korodense  relations  are  tabulated  in  Table  III  with 
the  data  from  the  CVA  67  technical  manual  [19]  included  for 
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comparison.  Both  relations  yielded  unsatisfactory  design 
verification  and  in  both  cases  a very  conservative  condenser 
was  designed. 

All  parameters,  as  calculated  by  0PC0DE2,  were  ten  to 
twenty  percent  less  than  those  from  the  actual  condenser. 

With  the  removal  of  subroutine  .ADJUST,  the  condenser  designed 
by  0PC0DE2  vented  in  e.xcess  of  20  percent  of  the  inlet  steam 
without  condensing  it. 

It  is  believed  that  the  effects  of  tube  inundation  give 
a shell-side  film  heat  transfer  coefficient  which  is  too 
conservative  as  reflected  in  the  reduced  value  of  the  overall 
heat  transfer  coefficient. 

D.  LIMITATIONS  OF  0PC0DE2 

Because  ORCONl  has  the  capability  to  make  design  decisions 
and  COPES/CONMIN  requires  a "once  through"  analysis,  the 
coupling  of  these  two  programs  in  OPCODE2  created  a situation 
where  the  programs  were  working  against  each  other.  This 
placed  a limitation  on  which  parameters  could  be  used  as 
design  variables  and  on  what  range  of  values  these  design 
variables  could  use. 

The  condenser  designed  by  0PC0DE2  was  very  conservative, 
a condition  caused  primarily  by  the  conservative  value  of 
the  overall  heat  transfer  coefficient  calculated. 

Accepting  these  limitations,  it  was  felt  that  OPCODE2 
should  be  further  developed  and  that  case  studies  should  be 
performed , 
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V.  RESULTS 


A.  EXPLANATION  OF  THE  CASE  STUDIES 

The  case  studies  were  devised  to  best  exercise  the 
attributes  of  OPCODEl  and  0PC0DE2,  and  were  made  as  realistic 
as  possible  so  as  to  simulate  the  problem  of  a condenser 
design  and  specification  during  the  early  stages  of  power 
plant  design.  Only  input  parameters  that  would  normally 
be  available  were  used. 

When  comparing  the  results  from  the  different  cases,  two 
cautions  must  be  kept  in  mind.  First,  since  all  the  cases 
involve  four  to  six  design  variables,  the  design  is  taking 
place  in  a four  to  six  dimensional  design  space  and  intuition 
on  how  an  optimized  design  "should"  turn  out  is  not  always 
applicable.  Secondly,  the  percentage  change  referred  to  in 
each  case  is  calculated  based  on  the  initial  design  for  that 
particular  case. 

1 . Constraint  Framework  for  OPCODEl 

In  order  to  simulate  an  actual  trade  off  study,  the 
constraints  for  each  case  study  were  kept  the  same,  even 
though  an  unimportant  constraint  could  become  active  during 
a particular  case  study.  In  this  way,  each  case  study  could 
be  directly  compared  with  all  others. 

The  main  condenser  for  the  CVA  67  was  to  be  designed 
with  a maximum  bundle  diameter  of  ten  feet;  a terminal 
temperature  difference  (pinch  point)  of  at  least  five  degrees 
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Fahrenheit  (°F)  but  not  more  than  35°F;  a cooling  water 
temperature  rise  of  at  least  five  °F,  but  not  more  than  20®F; 
and  a ratio  of  tube  sheet  hole  area  to  tube  sheet  area 
without  the  drilled  holes  of  less  than  0.36. 

The  constraint  on  bundle  diameter  was  chosen  only 
because  it  is  a reasonable  value.  If  more  or  less  vertical 
space  was  available  in  a proposed  machinery  arrangement,  the 
diameter  constraint  would  be  appropriately  changed.  The 
lower  constraint  on  pinch  point  came  from  reference  1 
which  called  for  a minimum  terminal  temperature  difference 
of  five  "F.  No  reference  to  an  upper  limit  on  pinch  point 
could  be  found,  but  several  technical  ma.ni'als  specified 
temperatures  in  the  range  of  20°F  to  45®F.  A value  of  35®F 
was  therefore  chosen  as  the  upper  bound  on  the  terminal 
difference . 

Reference  24  states  that  the  difference  in  temperature 
between  the  steam  and  cooling  water  streams  entering  the 
condenser,  or  temperature  range,  is  ordinarily  kept  to  20°F, 
and  that  the  cooling  water  temperature  rise  is  usually  made 
about  five  ®F  less  than  the  temperature  range.  Since  the 
CVA  67  e.xhausts  steam  with  a saturation  temperature  of  125°F 
and  reference  3 calls  for  a cooling  water  injection  tempera- 
ture of  75°F,  the  temperature  range  used  was  50®F  and  the 
cooling  water  temperature  rise  upper  limit  was  iS^F.  As 
the  guidance  given  by  reference  24  seemed  inappropriate  and 
since  neither  reference  1 nor  reference  2 specifies  limits 
on  cooling  water  temperature  rise,  the  lower  limit  was 
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arbitrarily  set  at  five  ®F  and  the  upper  limit  was 
arbitrarily  set  at  20°F.  The  amount  of  tube  sheet  material 
that  can  be  removed  by  drilling  for  the  installation  of 
tubes  is  specified  at  24  percent  of  the  blank  tube  sheet 
area  by  reference  2.  Since  OPCODEl  does  not  allow  for 
steam  lanes  in  the  design  of  a condenser,  and  since  these 
steam  lanes  would  provide  blank  tube  sheet  area,  the 
24  percent  limit  was  raised  to  36  percent. 

To  ensure  that  an  unrealistic  tube  wall  thickness 
was  not  specified,  a constraint  on  the  wall  thickness  was 
added.  The  wall  thickness  was  calculated  from  the  current 
values  of  tube  outside  diameter  and  tube  inside  diameter, 
and  the  constraint  applied.  Values  of  wall  thickness  in 
the  range  from  BWG  24  (0.022  inch)  to  BWG  12  (0.109  inch) 
were  used  as  the  lower  and  upper  constraints. 

In  summary,  the  general  design  constraints  and  the 
associated  upper  and  lower  bounds  were: 

0.022  £ tube  wall  thickness  (inch)  ^ 0.109 

1.0  £ bundle  diameter  (feet)  £ 10.0 

5.0  ^ terminal  temperature  difference  (°F)  £ 55.0 

5.0  £ sea  water  temperature  rise  (°F)  £ 20.0  . 

These  design  constraints  and  the  associated  upper  and  lower 
bounds  were  used  in  all  of  the  OPCODEl  test  cases  e.xcept 
where  specifically  modified. 
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Design  Variable  Framework  for  OPCODEl 


The  condenser  tube  outside  and  inside  diameters 
were  used  as  design  variables.  The  side  constraints  were 
set  to  correspond  with  the  values  of  normally  available 
tubes  [1].  The  tube  outside  diameter  was  allowed  to  vary 
in  the  range  between  0.625  inch  and  1.25  inch.  The  tube 
inside  diameter  was  allowed  to  vary  from  0.407  inch  to 
1.206  inch . 

The  tube  pitch  is  defined  as  the  center  to  center 
spacing  between  tubes.  The  pitch  to  diameter  ratio  (S/D) 
is  an  accurate  measure  of  how  closely  packed  the  tube 
bundle  is.  Generally  accepted  S/D  ratios  lie  in  the  range 
of  1.3  to  1.7.  However,  to  give  greater  latitude  to  this 
design  variable,  and  since  OPCODEl  was  insensitive  to 
shell-side  conditions,  the  S/D  ratio  was  allowed  to  vary 
within  the  range  from  1.1  to  3.0. 

There  was  no  guidance  available  on  the  range  of  tube 
lengths  that  would  be  applicable  to  a condenser  of  this  sice 
using  0.625  inch  tubes.  Reference  1 gave  a range  of 
recommended  tube  lengths  of  eight  to  fourteen  feet  for 
0.625  inch  outside  diameter  tubes  with  an  upper  limit  on 
heat  transfer  surface  area  of  1000  square  feet.  In  the 
range  of  the  expected  heat  transfer  area  of  12,000  to  18,000 
square  feet,  the  recommended  tube  outside  diameters  were 
from  0.75  inch  to  0.875  inch  with  tube  lengths  ranging 
from  16  to  24  feet.  Since  the  lower  bound  was  not  considered 
to  be  as  crucial  as  the  upper  bound,  it  was  set  at  one  foot. 
The  upper  bound  on  tube  length  was  set  at  25  feet. 
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Cooling  water  velocity  generally  ranges  from  three 
to  nine  feet  per  second  for  all  common  tube  materials 
except  titanium  which  has  an  upper  bound  of  15  feet  per 
second . 

In  summary,  the  general  design  variables  and  the 
associated  side  constraints  were: 

0.625  ^ tube  outside  diameter  (inch)  <_  1.25 

0.407  ^ tube  inside  diameter  (inch)  ^ 1.206 

1.1  £ pitch/diameter  ratio  ^ 3.0 

1.0  ^ tube  length  (feet)  ^ 25.0 

3.0  £ cooling  water  velocity  (feet/second)  < 9.0 

These  design  variables  and  the  associated  side 
constraints  were  used  in  all  of  the  OPCODEl  test  cases 
except  where  specifically  modified. 

3 . Constraint  Framework  for  0PC0DE2 

The  same  constraints  that  were  used  for  OPCODEl 
were  used  for  0PC0DE2  with  several  additional  constraints 
required  due  to  the  physical  calculations  performed. 

To  ensure  that  the  cooler  was  not  designed  to  be 
wider  than  the  void  inside  diameter,  the  ratio  of  cooler 
width  to  void  diameter  had  an  upper  bound  of  1.0  and  a lower 
bound  of  0.1 

As  the  initial  design  with  0PC0DE2  vented  over 
20  percent  of  the  inlet  steam  without  condensing  it,  an 
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upper  bound  on  the  exit  steam  fraction  was  set  at  five 
percent.  The  lower  bound  was  set  to  zero. 

Since  0PC0DE2  was  attempting  to  allow  for  a central 
steam  lane,  the  upper  limit  on  bundle  diameter  was  set  at 
12  feet  and  the  lower  limit  was  set  at  five  feet. 

Bundle  volume  had  an  upper  limit  of  1500  cubic 
feet  and  a lower  limit  of  100  cubic  feet.  These  limits 
had  been  used  satisfactorily  with  OPCODEl. 

To  ensure  that  the  cooler  height  did  not  become 
larger  than  the  band  of  condenser  tubes  from  the  outside 
diameter  of  the  bundle  to  the  inner  void,  the  ratio  of 
cooler  height  to  tube  band  width  was  calculated.  This 
ratio  had  an  upper  limit  of  1.0  and  a lower  limit  of  0.01. 

The  pinch  point  and  the  sea  water  temperature  rise 
had  the  same  range  as  was  used  for  OPCODEl. 

The  tube  wall  thickness  was  not  used  as  a design 
constraint  in  0PC0DE2  as  it  had  been  in  OPCODEl.  Instead 
the  tube  inside  diameter  was  used,  with  the  lower  constraint 
set  at  0.407  inches  and  the  upper  constraint  set  at  1.206 
inches  to  correspond  with  the  values  of  normally  available  [1]. 

In  summary,  the  general  design  constraints  and  the 
associated  upper  and  lower  bounds  were; 

0.1  £ cooler  width/void  diameter  ^ 1.0 

0.0  £ exit  steam  percentage  ^ 5.0 

5.0  < bundle  diameter  (feet)  < 12.0 
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100.  £ bundle  volume  (cubic  feet)  ^ 1500. 

0.01  ^ cooler  height/tube  band  width  ^ 1.0 

0.407  £ tube  inside  diameter  (inch)  ^ 1.206 

5.0  ^ terminal  temperature  difference  (°F)  £ 35. 

5.0  ^ sea  water  temperature  rise  (°F)  £ 20. 

4 . Design  Variable  Framework  for  OPCODE2 

The  condenser  tube  outside  diameter  and  wall 
thickness  were  used  as  design  variables.  The  side  con- 
straints were  set  to  correspond  with  the  values  of  normally 
specified  tubes  [1].  The  tube  outside  diameter  had  a 
lower  side  constraint  of  0.625  inch  and  an  upper  side 
constraint  of  1.25  inch.  The  tube  wall  had  a lower  side 
constraint  of  0.022  inch  and  an  upper  side  constraint  of 
0.109  inch. 

Tube  length  had  the  same  side  constraints  as  were 
used  for  OPCODEl. 

The  tube  pitch  to  diameter  ratio  had  a lower  side 
constraint  of  1.4  and  an  upper  side  constraint  of  2.0. 

This  band  of  allowable  values  was  reduced  from  that  used 
with  OPCODEl  because  of  instabilities  that  developed  for 
values  of  this  variable  outside  of  this  band.  More  prob- 
lems occurred  at  the  lower  end  of  the  range  than  at  the 
upper  end.  These  problems  originated  from  the  tubes  coming 
too  close  together  and  causing  an  excessively  high  steam 
velocity  with  a subsequent  large  pressure  drop  through  a 
row  of  tubes.  Since  condenser  pressure  was  already  low, 
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passage  through  very  few  rows  led  to  a negative  steam 
‘ pressure  and  thus  the  computational  problems  associated  with 

the  logarithmic  calculation  of  thermodynamic  properties 
discussed  in  Chapter  IV. 

The  sea  water  velocity  was  allowed  to  vary  within 
the  range  from  three  feet  per  second  to  ten  feet  per  second. 
The  upper  side  constraint  was  raised  from  the  nine  feet 
per  second  specified  in  OPCODEl  in  an  attempt  to  increase 
the  mass  flow  rate  of  the  cooling  water  and,  therefore,  the 
heat  rejection  rate. 

In  summary,  the  general  design  variables  and  the 
associated  side  constraints  were: 

i 

r 

5.0  £ tube  length  (feet)  <_  25.0 

3.0  ^ sea  water  velocity  (feet  per  second)  £ 10.0 
0.625  ^ tube  outside  diameter  (inch)  ^ 1.25 
0.022  ^ tube  wall  thickness  (inch)  £ 0.109 

1.4  £ tube  pitch/diameter  ratio  £ 2.0 

B.  CASE  STUDIES  USING  OPCODEl 
1 . Case  One 

The  objective  of  this  test  case  was  to  minimize 
the  pumping  power  requirement  while  holding  the  heat  load 
* to  the  condenser  constant.  The  input  parameters  are  pre- 

[ sented  in  Table  IV,  the  initial  design  is  presented  in 

Table  V,  and  the  results  of  the  optimization  are  presented 
' in  Table  VI. 
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These  results  show  an  84  percent  decrease  in  pumping 
power  with  a commensurate  37  percent  increase  in  heat 
transfer  surface  area,  a 55  percent  increase  in  condenser 
volume,  and  a decrease  of  27  percent  in  overall  heat  trans- 
fer coefficient.  The  log  mean  temperature  difference 
remained  constant. 

The  lower  bound  on  tube  wall  thickness  constraint 
and  the  upper  bound  on  sea  water  temperature  rise  constraint 
were  both  active.  .Mo  side  constraints  were  active  and  no 
constraints  were  violated. 

The  decrease  in  pumping  power  is  impressive  but 
so  are  the  increases  in  condenser  dimensions  with  the  implied 
increases  in  weight  and  cost. 

2 . Case  Two 

The  objective  of  this  test  case  was  to  minimize 
condenser  volume  with  the  heat  load  to  the  condenser  held 
constant.  The  input  parameters  are  presented  in  Table  IV, 
the  initial  design  is  presented  in  Table  V,  and  the  results 
of  the  optimization  are  presented  in  Table  VII. 

These  results  show  a 15  percent  decrease  in  condenser 
volume  with  an  unexpected  3.0  percent  decrease  in  pumping 
power.  This  design  can  be  understood  by  noting  that  the 
tube  wall  thickness  was  reduced  from  0.049  inch  to  0.022 
inch  causing  an  increase  in  material  correction  factor 
from  0.90  to  0.99.  This  increase  was  the  primary  factor 
in  increasing  the  overall  heat  transfer  coefficient  9.9 
percent . 
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The  lower  constraint  on  tube  wall  thickness,  and 
the  upper  constraint  on  tube  sheet  area  ratio  were  both 
active.  In  addition,  the  lower  side  constraint  on  tube 
outside  diameter  and  the  upper  side  constraint  on  sea  water 
velocity  were  both  active.  There  were  no  violated  constraints. 

3 . Case  Three 

The  objective  of  this  case  was  to  maximice  the  heat 
rejected  while  holding  pumping  power  constant.  Since  the 
pumping  power  is  a calculated  quantity,  the  method  of  com- 
bining pumping  power  and  the  heat  rejected  with  an  appro- 
priate weighting  factor  in  the  objective  function  was  used. 

The  objective  function  was  therefore: 

OBJ  = QREJ  ^ A * POWER 

and  POWER  was  added  as  a design  constraint  with  the  target 
value  of  68836  foot  pounds  per  second  as  the  upper  constraint. 
With  the  weighting  factor.  A,  set  at  5600,  the  pumping 
power  was  constant  within  one  percent.  Turbine  exhaust 
steam  and  auxiliary  exhaust  steam  were  both  added  as  design 
variables  to  provide  the  driving  force  to  increase  the  heat 
rejection  rate. 

The  input  parameters  are  presented  in  Table  IV, 
the  initial  design  is  presented  in  Table  V and  the  results 
of  the  optimization  are  presented  in  Table  VIII. 
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These  results  show  an  82  percent  increase  in  the 


heat  rejection  rate  with  an  accompanying  88  percent  increase 
in  surface  area.  The  log  mean  temperature  difference 
remained  constant  and  the  overall  heat  transfer  coefficient 
decreased  by  3.7  percent.  The  condenser  volume  increased 
by  ’8  percent  and  the  number  of  tubes  increased  by  96  percent. 

The  upper  constraint  on  pumping  power  was  active,  as 
desired,  and  the  upper  constraint  on  the  tube  sheet  area 
ratio  was  active.  There  were  no  violated  constraints  and 
no  active  side  constraints. 

It  should  be  noted  that  the  turbine  e.xhaust  steam 
rate  increased  by  83  percent,  whereas  the  auxiliary  exhaust 
steam  rate  increased  only  1.0  percent.  Inspection  of  the 
component  parts  of  QREJ  explains  the  dominance  by  the  main 
steam  rate.  The  rejected  heat  was  calculated  by: 


■rej 


■"ms  ‘^^MS 


"’as 


+ heat  from  other  sources 


where  the  subscript  MS  corresponds  to 
subscript  AS  corresponds  to  auxiliary 
other  sources  was  assumed  constant  dur 
The  product  m^jg  ^^MS  several  order 

than  the  product  m^g  ^^AS  there 

optimization  process  as  a change  in  th 
would  yield  a greater  design  improveme 
change  in  the  m^g  design  parameter. 
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In  all  probability,  the  m^g  design  variable  would 
not  begin  to  make  an  appreciable  move  until  the  design 

variable  approached  its  upper  side  constraint.  This  could 
not  occur  as  the  m^^^g  design  variable  had  an  unbounded  upper 
side  constraint. 

No  attempt  was  made  to  balance  the  flow  of  steam 
from  the  two  sources  to  more  evenly  distribute  the  incoming 
heat  load.  Balancing  could  be  accomplished  by  appropriately 
scaling  the  two  mass  flow  rates  so  they  would  have  approximately 
the  same  values. 

4 . Case  Four 

The  objective  of  this  case  was  to  minimize  the 
condenser  volume  while  holding  pumping  power  and  the  rejected 
heat  constant.  POWER  was  again  included  in  the  objective 
function 

OBJ  = VOL  + A * POWER 

and  as  a design  constraint  with  a lower  bound  of  68836  foot 
pounds  per  second.  The  lower  bound  was  the  target  value 
for  pumping  power.  Trial  and  error  solutions  with  various 
values  of  the  weighting  factor  showed  that  a value  of 
A “ 0.0  gave  a value  of  POWER  which  was  constant  within 
0.5  percent.  This  was  interpreted  to  mean  that  the  optimi- 
zation would  have  reduced  the  objective  function  further 
if  the  pumping  power  constraint  was  not  present  and  that 
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the  pumping  power  lower  constraint  was  therefore  always 
active.  With  the  value  of  this  constraint  set  at  the 
target  value  for  constant  pumping  power,  inclusion  of  POWER 
in  the  objective  function  was  unnecessary  and  A = 0.9  was 
appropriate . 

The  input  parameters  are  presented  in  Table  IV, 
the  initial  deisgn  is  presented  in  Table  V and  the  results 
of  the  optimization  are  presented  in  Table  IX. 

These  results  show  only  a 2.0  percent  decrease  in 
condenser  volume.  The  tube  pitch  to  tube  diameter  ratio 
remained  essentially  constant  as  did  the  tube  outside  diameter. 
The  tube  wall  thickness  increased  by  11  percent  with  an 
accompanying  decrease  in  material  correction  factor  from 
0.90  to  0.88. 

The  lower  constraint  on  pumping  power  was  active, 
as  discussed  above.  The  upper  constraint  on  sea  water 
temperature  rise,  as  well  as  the  upper  constraint  on  tube 
sheet  area  ratio  were  active.  The  upper  constraint  on  sea 
water  velocity  was  the  only  active  side  constraint.  There 
were  no  violated  constraints. 

The  results  of  this  case  study  indicated  that  the 
circular  bundle  prototype  was  very  close  to  having  an  opti- 
mum volume  within  the  given  design  variable  and  design 
constraint  framework. 

Comparison  of  Cases  Two  and  Four  shows  that  a 
smaller  condenser  with  essentially  constant  pumping  power 
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was  achieved  with  Case  Two  then  with  Case  Four.  This 
anomaly  can  be  understood  by  comparing  the  constraints 
for  the  two  cases. 

Pumping  power  was  not  a constraint  for  Case  Two 
and  it  was  an  active  constraint  for  Case  Four.  Thus, 
the  Case  Four  optimization  had  to  contend  with  an  additional 
constraint  whereas,  Case  Two  had  more  latitude  within  the 
design  space  and  a better  optimum  was  found. 

5 . Case  Five 

The  objective  of  this  case  was  to  minimize  the 
condenser  overall  length  while  holding  pumping  power  and 
the  heat  rejected  constant.  Pumping  power  was  included 
in  the  objective  function: 

OBJ  = CLOA  + A * POWER 

and  as  a design  constraint  with  a lower  bound  of  68856  foot 

pounds  per  second  as  the  target  value.  With  a weighting 

- 4 

factor  of  3.526  x 10  , pumping  power  was  held  constant 

within  0.5  percent . 

The  input  parameters  are  presented  in  Tabe  IV,  the 
initial  design  is  presented  in  Table  V and  the  results  of 
the  optimization  are  presented  in  Table  X. 

The  overall  condenser  length  was  reduced  by  28  percent. 
The  heat  transfer  area  increased  by  9.7  percent  and  the 
overall  heat  transfer  coefficient  decreased  by  15  percent, 
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while  the  log  mean  temperature  difference  increased  by 
".1  percent.  The  number  of  tubes  increased  96  percent  and 
the  tube  length  decreased  44  percent.  The  tube  wall  thick- 
ness increased  causing  a decrease  of  14  percent  in  the 
material  correction  factor  thereby  contributing  to  the 
decrease  in  overall  heat  transfer  coefficient. 

•Accompanying  the  23  percent  decrease  in  condenser 
length  was  a compensating  37  percent  increase  in  bundle 
diameter  and  an  18  percent  increase  in  condenser  volume. 

The  lower  constraint  on  pumping  power  was  active, 
as  desired.  The  upper  constraints  on  pinch  point  and  on 
tube  sheet  area  ratio  were  act.ive.  There  were  no  active 
side  constraints  and  no  violated  constraints. 

6 . Case  Six 

The  objective  of  this  case  was  to  minimize  pumping 
power  while  holding  condenser  volume  and  the  heat  rejected 
constant.  For  this  case,  volume  was  included  in  the 
objective  function: 

OBJ  * POWER  > .A  * VOL 


and  it  was  also  included  as  a design  constraint  with  an 
upper  bound  of  928.8  cubic  feet  as  the  target  value.  The 
usual  trial  and  error  procedure  with  various  weighting 
functions  was  performed  and  the  best  results  were  obtained 
for  a weighting  factor  of  -2.5.  With  this  value,  condenser 
volume  was  held  constant  to  within  1.3  percent. 
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The  input  parameters  are  presented  in  Table  IV, 
the  initial  design  is  presented  in  Table  V and  the  results 
of  the  optimization  are  presented  in  Table  XI. 

The  pumping  power  was  reduced  by  55  percent.  The 
log  mean  temperature  difference  and  the  heat  transfer  area 
remained  essentially  constant,  while  the  tube  wall  thickness 
decreased  52  percent  causing  a ten  percent  increase  in  the 
tube  material  correction  factor.  The  sea  water  velocity 
was  reduced  by  16  percent,  however,  and  the  overall  heat 
transfer  coefficient  remained  essentially  constant.  There 
were  no  active  constraints,  no  violated  constraints  and 
no  active  side  constraints. 

Volume  was  not  forced  against  its  upper  constraint 
as  planned  because  the  two  members  of  the  objective  function 
were  attempting  to  move  in  opposite  directions.  As  power 
was  decreased,  the  volume  tended  to  increase.  The  negative 
weighting  factor  helped  to  turn  this  process  around  but 
the  optimum  was  reached  before  the  constraint  was  activated. 

7 . Case  Seven 

English  has  shown  [25]  that  by  placing  boundary 
layer  fences  in  front  of  the  injection  scoop  of  a ship, 
the  pressure  coefficient,  and  hence  pumping  power,  will 
increase  by  68  percent.  These  boundary  layer  fences  cause 
shedding  vortices  to  form  behind  them  thus  pumping  water 
from  the  boundary  layer  around  the  ship  into  the  condenser 
intake . 
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This  case  takes  advantage  of  the  increase  in 
available  pumping  power  by  setting  pumping  power  at  a 
value  50  percent  greater  than  the  normal  68336  foot  pounds 
per  second,  and  holding  this  new  value  of  pumping  power 
constant  while  maximizing  the  rejected  heat. 

The  constraints  and  design  variables  were  set  up 
exactly  as  they  were  for  Case  Three  with  the  exception  that 
the  upper  constraint  on  pumping  power  was  raised  to  103,250 
foot  pounds  per  second.  POWER  was  included  in  the  objective 
function  as  before  and  a weighting  factor  of  4095  yielded 
pumping  power  constant  within  0.5  percent. 

The  input  parameters  are  presented  in  Table  IV,  the 
initial  design  is  presented  in  Table  V and  the  results 
of  the  optimization  are  presented  in  Table  XII. 

The  heat  rejected  was  increased  by  107  percent.  The 
heat  transfer  area  increased  by  101  percent,  and  the  log 
mean  temperature  difference  remained  constant.  The  tube 
wall  thinned  down  by  53  percent  causing  an  increase  in  the 
tube  wall  material  correction  factor  of  ten  percent.  The 
sea  water  velocity  decreased  13  percent  and  the  overall  heat 
transfer  coefficient  increased  2.5  percent. 

The  upper  constraints  on  pumping  power  and  the  tube 
sheet  area  ratio  were  active.  There  were  no  active  side 
constraints  and  no  violated  constraints. 
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C.  CASE  STUDIES  USING  0PC0DE2 
1 . Case  Eight 

This  case  was  run  assuming  plain  copper -nickel 
tubes  with  a 14  percent  enhancement  factor  applied  to  the 
outside  film  heat  transfer  coefficient  as  described  in 
Section  IV(C)  of  this  work.  Steam  baffles  were  specified 
and  the  inlet  steam  flow  rate  was  set  at  200,000  pounds  per 
hour  as  only  one  half  of  a condenser  bundle  was  being  simu- 
lated. The  Korodense  equation  from  reference  [23]  was  used 
to  calculate  the  correction  factor  for  tube  inundation  by 
condensate  rain. 

The  objective  of  this  case  was  to  minimize  condenser 
volume.  The  initial  and  optimum  designs  are  presented  in 
Table  XVI. 

The  volume  was  reduced  13  percent  and  the  heat 
transfer  area  increased  14  percent.  The  log  mean  tempera- 
ture difference  decreased  5.6  percent,  the  overall  heat 
transfer  coefficient  increased  11  percent  and  the  heat 
rejected  increased  21  percent.  The  vented  steam  rate 
decreased  76  percent  to  five  percent  of  the  inlet  steam 
flow  rate.  Finally,  pumping  power  increased  65  percent 
and  the  sea  water  flow  rate  increased  35  percent. 

The  upper  side  constraint  on  sea  water  velocity  as 
well  as  the  lower  side  constraints  on  tube  outside  diameter, 
tube  pitch  to  diameter  ratio  and  tube  wall  thickness  were 
all  active.  The  upper  constraint  on  vented  steam  rate 
was  active.  There  were  no  violated  constraints. 
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Case  Mine 


The  objective  of  this  case  was  to  minimize  the 
condenser  volume  when  enhanced  tubes  were  specified. 

The  initial  parameters  were  the  same  as  those  used  in 
Case  Mine  with  the  e.xception  of  the  outside  enhancement 
factor  which  was  set  to  1.8.  It  was  felt  that  this  modest 
enhancement  was  attainable  with  the  augmented  tubes 
currently  being  manufactured. 

The  initial  and  optimum  designs  are  presented  in 
Table  XVII. 

Use  of  enhanced  tubes  reduced  the  condenser  volume 
by  21  percent.  The  heat  transfer  area  increased  2.6  per- 
cent, the  overall  heat  transfer  coefficient  increased  3-i 
percent,  the  log  mean  temperature  coefficient  decreased 
11  percent  and  the  heat  rejected  increased  21  percent.  The 
vented  steam  rate  decreased  "6  percent  to  five  percent  of 
the  inlet  steam  flow  rate. 

Pumping  power  increased  52  percent  and  the  cooling 
water  flow  rate  increased  35  percent.  Tube  length  was 
increased  by  2.7  percent  and  the  tube  inside  diameter 
increased  ten  percent. 

The  pumping  power  calculated  was  based  on  smooth 
tubes.  The  most  common  augmented  condenser  tubes  use  either 
a rope  design,  a twist  design,  or  internal  flow  promoters 
to  enhance  heat  transfer.  Reilly  [26]  has  found  that  the 
pressure  drop  through  enhanced  tubes  is  1.5  to  3.0  times  the 
pressure  drop  through  smooth  tubes  for  the  same  flow  conditions. 
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The  upper  side  constraint  on  sea  water  velocity, 
as  well  as  the  lower  side  constraints  on  tube  outside 
diameter,  tube  pitch  to  diameter  ratio  and  tube  wall 
thickness  were  all  active.  The  upper  constraint  on  vented 
steam  rate  was  active,  and  there  were  no  violated  constraints. 

3 . Case  Ten 

This  case  was  run  using  the  same  conditions  as 
Case  Eight  with  the  e.xception  that  the  original  equation 
from  ORCO.Nl  was  used  to  correct  for  tube  inundation.  Reference 
3 specified  a value  for  the  flooding  factor,  FDAVE , of 
zero  for  pitch  to  diameter  ratios  greater  than  1.4. 

This  value  was  used  during  this  run  of  0PC0DE2. 

The  objective  of  this  case  was  to  minimize  the  condenser 
volume.  The  initial  and  optimum  designs  are  presented  in 
Table  XVIII. 

.A.  more  conservative  design  was  produced  with  the 
Eissenberg  equation  than  the  design  produced  with  the 
Korodense  equation.  The  volume  was  reduced  5.3  percent  and 
the  heat  transfer  area  increased  24  percent.  The  log  mean 
temperature  difference  remained  constant,  the  overall  heat 
transfer  coefficient  increased  by  3.9  percent,  and  the 
heat  rejected  increased  by  33  percent.  The  vented  steam 
rate  decreased  83  percent  to  five  percent  of  the  inlet 
steam  flow  rate.  Pumping  power  increased  75  percent  and  the 
sea  water  flow  rate  increased  35  percent. 

The  upper  side  constraint  on  sea  water  velocity, 
as  well  as  the  lower  side  constraints  on  tube  outside 
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diameter,  tube  pitch  to  diameter  ratio  and  tube  wall 
thickness  were  all  active.  The  upper  constraint  on  vented 
steam  rate  was  active.  There  were  no  violated  constraints. 
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CONCLUSIONS 


VI . 

The  intent  of  this  investigation  was  to  couple  a 
numerical  optimization  code  with  both  a simple  computer 
code  for  condenser  design  based  on  the  HEI/DDS  method 
and  the  more  sophisticated  ORCONl , and  to  test  the  programs 
'.eveloped  with  a variety  of  test  cases  to  prove  their 
viability.  The  results  of  these  test  cases  were  presented 
in  Section  V;  the  resulting  conclusions  are  summarized 
here . 

A.  A condenser  designed  by  the  HEI/DDS  method  is  very 
near  the  volumetric  optimum  in  design,  as  only  a 
15  percent  decrease  in  condenser  volume  was 
realized  when  volume  was  optimized  using  OPCODEl. 

B.  OPCODEl  is  an  excellent  design  tool  for  the  concep- 
tual design  of  a condenser  to  meet  specified  con- 
straints such  as  length,  height,  and  weight. 

C.  OPCODEl  can  optimize  a variety  of  objective  func- 
tions, and  with  some  trial  and  error  application 
of  weighting  factors,  equality  constraints  can 

be  met. 

D.  OPCODEl  will  yield  an  optimum  design  on  which 
experience-based  safety  factors  may  be  applied. 

E.  The  ability  to  enhance  tubes  and  the  sensitivity 
of  the  program  to  shell-side  conditions  make  the 
ORCONl-based  0PC0DE2  very  attractive. 
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F.  The  basic  equations  and  the  research  applied  to 
the  development  of  ORCONl  are  sound,  but  a conser- 
vatively-designed condenser  results.  More  baffles 
are  required  to  more  closely  approximate  an  actual 
condenser  and  to  aid  in  reducing  the  inundation 
effect  on  interior  tubes. 

G.  The  0PC0DE2  designed  condensers  were  more  conserva- 
tive than  those  designed  by  OPCODEl.  This  was 
caused  by  the  detrimental  effect  on  outside  film 
coefficient  by  tube  flooding.  It  is  felt  that  this 
effect  could  be  reduced  by  the  addition  of  more 
baffles . 
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VI  I.  RECOMMENDATIONS 


In  addition  to  the  insight  that  this  investigation  has 
given  into  the  generation  of  automated  design  programs  for 
condenser  design,  it  has  also  generated  an  awareness  o 
this  investigation’s  shortcomings.  Presented  herein  ai 
recommendations  for  improving  upon  and  furthering  de'  lop- 
ment  of  the  OPCODE  family  of  condenser  design  progr' 

A.  OPCODEl  should  be  expanded  to  include  the  ility 

for  multi-pass  condenser  calculations,  g - .al 
options,  material  costing  equations,  stre  .n  of 
material  considerations,  and  an  algorithm  to  compute 
steam  velocity  through  the  rows  of  tubes  to  ensure 
that  the  velocity  remains  realistic. 

B.  Research  should  be  performed  to  develop  enhancement 
factor  data  for  various  tube  types  that  can  be 
applied  to  the  basic  equations  used  in  OPCODEl. 

C.  ORCONl  should  be  rewritten  with  the  intent  of 
applying  optimization  techniques  to  the  new  program, 
and  based  on  the  equations  and  logic  developed  by 
Eissenberg,  Korodense,  and  other  investigators  and 
manufactures,  with  the  option  to  choose  the 
desired  relationship  to  be  used. 

D.  In  an  effort  to  reduce  the  computer  time  required 
for  a run  using  0PC0DE2  (typically  35  CPU  minutes), 
an  investigation  into  using  OPCODEl  as  a pre- 
processor in  front  of  a detailed  analysis  scheme  is 
attractive. 
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E.  Expand  0PC0DE2  to  perform  calculations  for  other 


than  a circular  bundle. 
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FIGURE  2.  Two-Variable  Design  Space  Showing  the  Initial 
Design  at  Point  . 
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TUBE  PITCII/DIAMETER  RATIO 


TUBE  PITCH/DIAMETER  RATIO 


SIDE  CONSTRAINT 


FIGURE  4.  Two-Variable  Design  Space  Showing  the  Second 
Design  Iteration 
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FIGURE  5,  Two-Variable  Design  Space  Showing  the  Third 
Design  Iteration 


TUBE  PITCH/DIAMETER  RATIO 


FIGURE  6,  Two-Variable  Design  Space  Showing  the  Fourth 
Design  Iteration 
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FIGURE  7.  Two-Variable  Design  Space  Showing  the  Fifth 
Design  Iteration 
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TUBE  PITCH/DIAflETER  RATIO 


SIDE  CONSTRAINT 
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FIGURE  3 


Two -Variable  Design  Space  Showing  the  Meed  for 
Optimization  Techniques 
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FIGURE  9 


Two -Variable  Design  Space  Showing  Relative 
Minima  for  a Constrained  Minimization  Problem 
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FLOW 


FIGURE  10,  Circular  Condenser  Bundle  Designed  by  ORCONl 


^ MAIM  ^ 


FIGURE  13.  Flowchart  from  ORCONl 
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CALL 

INPUT 


REMOVED 


FIGURE  14,  Flowchart  Showi.Tg  Modification  of  ORCONl 

to  opcode: 
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IX.  TABLES 


TECHNICAL 

MANUAL 

OPCODEl 

DEVIATION 

C%) 

Heat  transfer  area;  ft^ 

16,011 

16,099 

+ 0.55 

Number  of  tubes 

6,612 

6,648 

+ 0.54 

Bundle  diameter;  ft 

N/A 

7.26 

- 

Heat  rejected;  BTU/hr 

4.035x10^ 

4.035x10® 

0.00 

Overall  heat  transfer 
coefficient • 

BTU/(hr) (ft-) (“F) 

635 

637 

+ 0.31 

Log  mean  temperature 
difference;  °F 

39.7 

39.3 

-1.01 

Sea  water  temperature  rise; 
op 

19.9 

19.9 

0.0 

Terminal  temperature 
difference;  °F 

30.53 

30.2 

-1.08 

Sea  water  flow  rate,  gpm 

40,565 

40,645 

+ 0.20 

Tube-side  pressure  drop; 
ft . w. c . 

16.  3 

12.4 

-23.9 

TABLE  I.  Verification  of  the  CV.A  67  Main 
Condenser  Design  Using  OPCODEl 
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TECHNICAL 

MANUAL 

OPCODEl 

DEVI.ATION 

C%) 

Heat  transfer  area;  ft^ 

6,600 

6,527 

-1.11 

Number  of  tubes 

3,892 

3,850 

-1.08 

Bundle  diameter;  ft 

5.9 

4.6 

-22.1 

Heat  rejected;  BTU/hr 

2.05x10® 

2.047x10® 

-0.15 

Overall  heat  transfer 
coefficient • 

BTU/ [hr) (ft-) C^F) 

6 3 5 

637 

+ 0 . 31 

Log  mean  temperature 
difference;  °F 

50.0 

49 . 2 

-1.6 

Sea  water  temperature  rise; 
°F 

16.9 

17.4 

+ 2.96 

Terminal  temperature 
difference;  ®F 

41.9 

41.0 

-2.15 

Sea  water  flow  rate;  gpm 

23,900 

23,539 

-1.51 

Tube-side  pressure  drop; 
ft . w. c . 

11.5 

8.21 

-28.6 

TABLE  II.  Verification  of  the  DE1040  Class 
of  Main  Condenser  Using  OPCODEl 
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INITIAL 

VALUE  AT 

CHANGE 

VALUE 

OPTIMUM 

(%) 

Heat  transfer  area;  ft 

16,008 

.Number  of  tubes 

6,612 

Bundle  diameter;  ft 

10.23 

Heat  rejected;  BTU/hr 

3.241x10 

Overall  heat  transfer 
coefficient 

BTU/ (hr) (ft2)  (“F) 

562 

Log  mean  temperature 
difference;  °F 

36 

Sea  water  temperature 
rise;  ®F 

15.6 

Terminal  temperature 
difference;  *F 

25.6 

Sea  water  flow  rate;  gpm 

40,341 

Tube-side  pressure  drop; 
ft.  w.c. 

12.3 

Tube  length;  ft 

14.8 

Tube  outside  diam. ; in. 

0.625 

Tube  inside  diam.;  in. 

0.527 

Tube  wall;  in. 

0.049 

S-W  velocity;  fps 

9.0 

Pitch/diameter 

1.6 

Exit  steam;  % of  input 

21 

Pumping  power;  ft'lb/sec 

35,409 

Condenser  volume;  ft^ 

775 

18,290 

6,612 

8.95 

3.916x10 

625 

34 

13.7 

24.6 

54,479 

16.9 

0.625 


TABLE  XIII.  Case  Eight  Initial  Design  and 
Optimization  Results 
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INITIAL 

VALUE 


VALUE  AT 
OPTIMUM 


CHANGE 

(%) 


Heat  transfer  area;  ft 

16,008 

16.431 

Number  of  tubes 

6,612 

6,612 

Bundle  diameter;  ft 

10.23 

8.95 

Heat  rejected;  BTU/hr 

3.241x10® 

3.924x10 

Overall  heat  transfer 
coefficient, 

BTU/(hr) (ft^) (°F) 

562 

751 

Log  mean  temperature 
difference;  °F 

36 

1 

32 

1 

Sea  water  temperature 
rise;  °F 

15.6 

13.1 

Terminal  temperature 
difference;  °F 

25.6 

21.6 

Sea  water  flow  rate;  gpm 

40,341 

1 

54,479 

Tube-side  pressure  drop; 
ft.  w.c. 

12.3 

1 

1 

Tube  length;  ft 

14.8 

15.2 

Tube  outside  diam. ; in. 

0.625 

0.625 

Tube  inside  diam. ; in. 

0 . 527 

0.581 

Tube  wall;  in. 

0.049 

0.022 

S-W  velocity;  fps 

9 

10 

Pitch/diameter 

1.6 

1.4 

Exit  steam;  i of  input 

21 

1 

S 

Pumping  power;  ft*ib/sec 

35,409 

53,674 

Condenser  volume;  ft^ 

775 

609 

TABLE  XIV.  Case  Nine  Initial  Design  and 
Optimization  Results 
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INITIAL 

VALUE 

VALUE  AT 
OPTIMUM 

CHANGE 

(%) 

Heat  transfer  area;  ft^ 

16,008 

19,808 

+ 24 

Number  of  tubes 

6,612 

6,612 

0.0 

Bundle  diameter;  ft 

10.23 

8.95 

-13 

Heat  rejected;  BTU/hr 

2.937x10® 

3.910x10® 

+ 33 

Overall  heat  transfer 
coefficient : 

BTU/ (hr) (ft^) (°F) 

505 

550 

+ 8.9 

Log  mean  temperature 
difference;  °F 

36 

36 

0.0 

Sea  water  temperature 
rise;  ° F 

14.4 

14.2 

-1.4  i 

Terminal  temperature 
difference;  °F 

25.2 

26.3 

+ 4.4  1 

! 

Sea  water  flow  rate;  gpm 

40.341 

54,479 

+ 35 

Tube-side  pressure  drop; 
ft . w. c . 

12.3 

16.0 

+ 30 

Tube  length;  ft 

14.8 

18.3 

+ 24 

Tube  outside  diam.  ; in. 

0.625 

0.625 

0.0 

Tube  inside  diam. ; in. 

0.527 

0.581 

+ 10 

Tube  wall;  in. 

0.049 

0.022 

+ 11 

S-W  velocity;  fps 

9 

10 

+ 11 

Pitch/diameter 

1.6 

1 . 4 

-13 

Exit  steam;  % of  input 

29 

5 

-83 

Pumping  power;  ft*lb/sec 

35,409 

62,059 

+ 75 

Condenser  volume  ; ft^ 

775 

1 

734 

-5.3 

TABLE  XV.  Case  Ten  Initial  Design  and 
Optimization  Results 
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APPENDIX  A 


DEVELOPMENT  OF  THE  ANALIZ  SUBROUTINE  FOR  OPCODEl 

OPCODEl  was  based  on  the  HEI/DDS  method  of  calculation 
of  heat  transfer  and  thermodynamic  parameters. 

The  following  equations  are  used  in  the  design  of 
a condenser  using  the  HEI/DDS  method. 


^c  = 

F,F2F3C/7 

CA-l) 

Q = 

CA-2) 

Q * 

W Ah 

(A- 3) 

Q » 

500  G (t^  - tp 

(A-4) 

^0  - ^ 

in 

S 0 

(A- 5) 

* 

L N s 

(A-6) 

G - 

Ng  V 

(A-7) 

k - 

s/g 

(A-8) 

The  ANALIZ  subroutine  of  OPCODEl  was  divided  into  three 
sections:  input,  analysis  and  output.  The  correct  section 
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of  ANALIZ  was  entered  by  testing  for  the  value  of  ICALC 
that  was  passed  from  COPES  to  ANALIZ. 

When  ICALC  = 1,  the  input  section  of  ANALIZ  was  entered 
and  the  input  variables  and  problem  identification  were 
read.  The  input  section  had  the  capability  to  default  the 
latent  heat  of  vaporization  to  950  BTU  per  pound  as  speci- 
fied by  reference  [2]  and  to  calculate  either  the  saturation 
temperature  or  the  saturation  pressure  depending  on  which 
value  was  not  read  in  on  the  data  card. 

IVhen  COPES  set  ICALC  = 2,  the  analysis  section  of 
ANALIZ  was  entered  and  the  calculation  of  the  outside 
area  of  a condenser  tube  per  foot  of  tube  length  was  made: 

2 

s = IT  D ■ 

The  rate  of  flow,  in  gallons  per  minute,  through  a 
condenser  tube  at  a velocity  of  one  foot  per  second  was 
then  made: 

G1  - (J)  (d^)  (60  x 7.481)  = 448.46  (J)d‘ 


The  value  of  C in  equation  (A-1)  was  assumed  to  be  a 
constant  270  when  in  fact,  C was  weakly  dependent  on  tube 
outside  diameter.  With  this  assumption,  the  uncorrected 
heat  transfer  coefficient  was  calculated  as: 
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u 


C/V  [ 


BTU 

hr-ft“. °F 


] 


A call  was  then  made  to  subroutine  NIATFAC  to  find  the 
new  value  of  the  material  correction  factor,  F2,  based  on 
the  current  value  of  tube  wall  thickness.  The  temperature 
correction  factor,  F^,  was  found  with  the  function  sub- 
routine TEMFAC . 

The  data  from  Figure  1 of  reference  [2]  was  implemented 
in  TEMFAC  and  a value  of  F^  as  a function  of  injection 
temperature  was  retrieved  using  subroutine  I.NTRPL,  a systems 
supplied  interpolator. 

The  fouling  factor,  F^^ , was  set  to  0.85  in  accordance 
with  references  1 and  2 during  the  input  section  of  ANALIZ. 

The  corrected  value  of  the  overall  heat  transfer 
coefficient  was  calculated: 


Uc  = F^  F^  F3  U 


The  original  form  of  equation  (.A-4)  was 


Q = 60  W C G (t^  - t.)  . 


For  sea  water,  the  following  value  of  60 W Cp  was 


calculated : 


60(W  Cp)s.;y  * (60)  (8.  55)  (0.94)  » 482 
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while,  for  fresh  water,  60 


W C was  calculated  as: 
P 


60(W  Cp)p,,^  = (60)  (3. 33)  (1.0)  = 500  . 


Reference  [2]  recoioinends  the  use  of  60  W C = 500  in 

P 

keeping  with  industry  standards.  This  approximation  will 
induce  an  error  of  approximately  one-half  percent. 

Be  setting  equation  (A-2)  equal  to  equation  (A-4) 
and  substituting  for  values  of  G,  A^^  and  the  following 

relation  was  obtained: 


where 


temperature 


t 

r ^ 

t . 

1 

U - 

^0^ 

s 

written  as 

t 

- t . 

s 

1 

t 

- t 

s 

0 

. '^c 

L K 

- 500  V 

(A-9) 

was 

^0 

=*  t 

L K 

500  V 


(A-9) 


^s  - ^i 
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The  sea  water  temperature  rise  and  the  pinch  point  were 
calculated,  as  was  the  initial  temperature  difference.  All 
the  factors  were  now  available  for  the  calculation  of  the 
logarithmic  mean  temperature  difference,  using 

equation  (A- 5) . 

The  input  heat  load  was  divided  into  the  auxiliary 
steam  flow,  main  steam  flow,  and  heat  from  other  sources. 

The  auxiliary  and  main  steam  flows  were  multiplied  by 
their  respective  changes  in  latent  heat.  Ah,  and  summed 
with  the  heat  from  other  sources  to  yield  the  total  heat 
input,  Q. 

If  an  input  value  for  Ah  is  not  read  in  for  either 
auxiliary  or  main  steam  flow,  the  default  value  of  950  BTU 
per  pound,  specified  by  reference  [2],  is  used. 

With  the  value  of  Q,  equation  (A-4)  was  utilized  to 
calculate  the  required  flow  rate  of  cooling  water,  G.  With 
the  value  of  G known,  the  heat  transfer  area,  Aj^ , was 
found  using  equations  (A-6),  (A-7)  and  (A-8): 

Ah  -t—  • 

The  heat  rejected  to  the  cooling  water  was  found  by 
using  equation  (A-2). 

The  next  portion  of  ANALIZ  was  devoted  to  the  calcu- 
lation of  bundle  geometry.  The  call  to  subroutine  GEOM 
yielded  a condenser  design  of  circular  bundle  cross  section 
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with  a 12  inch  diameter  void  along  the  longitudinal  center- 
line  to  serve  as  a header  for  the  removal  of  nonconden- 
sable gasses  by  an  air  ejector.  The  rows  of  tubes  were 
filled  from  the  inner  void  and  were  concentric  to  the 
center  void.  Partial  tubes  were  permitted  in  order  to 
simplify  the  algorithm  and  the  final,  outermost  row  was 
only  partially  filled.  The  bundle  radius  was  calculated 
with  no  allowance  for  circumferential  steam  lanes,  since 
it  is  beyond  the  capability  of  the  HEI/DDS  method  to  provide 
for  steam  lanes. 

It  was  assumed  that  the  waterboxes  were  hemispherical 
caps  on  the  ends  of  the  cylindrical  bundle.  The  waterboxes 
may  not  be  any  deeper  than  45  inches  [2]  and  logic  was 
included  to  ensure  that  this  specification  was  met. 

The  ratio  of  tube  sheet  material  removed  for  tube 
installation  to  original  undrilled  tube  sheet  area  was 
calculated  and  used  as  a design  constraint  during  the 
optimization  process.  Reference  [2]  requires  that  the 
ratio  calculated  must  be  less  than  or  equal  to  0.24  to 
allow  for  adequate  tube  sheet  strength. 

A hotwell  capacity  capable  of  receiving  the  condensate 
from  one  minute  of  full  power  operation  is  specified  by 
reference  [2].  This  value  was  calculated  as  the  product 
of  the  specific  volume  of  the  incoming  steam  and  the 
steam  flow  rate. 

The  calls  to  subroutines  FRIFLAC  and  PRSDRP  found  the 
factors  required  for  the  calculation  of  tube  side  pressure 
drop  and  pumping  power. 
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Total  condenser  volume  was  the  sum  of  the  tube  bundle 

I* 

I volume,  the  waterbox  volume,  and  the  hotwell  volume  with 

the  individual  watervoxes  treated  as  spherical  caps. 

The  final  step  in  the  analysis  section  of  subroutine 
ANALIZ  was  to  define  the  objective  function  with  any 
necessary  weighting  factors  as  described  in  Chapter  II. 

COPES  repeatedly  enters  the  ANALIZ  subroutine  with 
ICALC  = 2 during  the  optimization  process  — often  on  the 
order  of  hundreds  of  times.  Therefore,  no  WRITE  statements 
were  included  in  the  analysis  section.  All  output  was 
performed  within  the  output  portion  of  ANALIZ,  a region 
that  is  entered  only  when  the  optimization  process  is 
' terminated  and  ICALC  is  set  equal  to  three. 
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APPENDIX  B 


DEVELOPMENT  OF  SUPPORTING  SUBROUTINES  FOR  OPCODEl 

In  this  appendix,  the  subroutines  that  were  utilized 
during  the  execution  of  ANALIZ  in  OPCODEl  are  briefly 
described. 

FRIFAC 

This  subroutine  iteratively  solved  the  transcendental 
Colebrook  equation  [17]  : 


/F 


7 1 r2.51  ^ 

-2  log  ( + 

Re/T 


1771- 


for  the  value  of  the  internal  friction  factor  for  tube 
flow.  FRIFAC  was  valid  for  the  range 


0.01  < f < 0.10  . 


This  method  was  chosen  over  the  simpler  Blasius  equation 
because  of  the  inclusion  of  the  desirable  dependence  on 
the  relative  roughness,  e/d,  in  the  Colebrook  equation. 

GEOM 

This  subroutine  calculated  bundle  geometry.  The 
condenser  bundle  was  assumed  to  have  a 12  inch  central  void 
along  the  longitudinal  centerline  with  the  rows  of  tubes 
in  concentric  rows  about  the  void.  The  void  served  as  a 
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header  for  the  venting  of  noncondensable  gasses.  The  rows 
were  filled  from  the  inside  out;  partial  tubes  were  permitted 
and  the  outermost  row  was  unfilled. 

PRSDRP 

This  subroutine  calculated  the  cooling  water  pressure 
drop  through  the  condenser  bundle.  The  development  of  this 
subroutine  was  thoroughly  discussed  in  Chapter  III  and  will 
not  be  repeated  here. 

DENSE 

This  subroutine  calculated  the  density  of  water  as  a 
function  of  temperature  with  the  relation  from  reference  [27] 

p = 63.8  - 0.01781  t 1.132  x lO'^  t^ 

- 6.786  X 10'^  t^  . 

Based  on  the  assumption  made  for  determining  the  coeffi- 
cient for  equation  (A-4)  in  Appendix  A,  the  relation  given 
above  was  used  for  the  sea  water  as  well  as  the  steam  and 
condensate  densities. 

For  the  calculation  of  cooling  water  density,  the 
numerical  average  of  the  cooling  water  injection  temperature 
and  the  cooling  water  overboard  temperature  was  used. 

PSATFN 

This  function  subroutine  calculated  the  saturation 
pressure  of  steam  as  a function  of  the  saturation  temperature 
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I 


in  degrees  Rankine,  using  the  relation  from  reference  [3]; 


TEMFAC 

The  retrieval  of  the  temperature  correction  factor, 

Fj,  was  performed  in  this  function  subroutine.  The  data 
presented  in  Figure  SF-2  of  reference  [1]  was  implemented 
in  tabular  form  and  was  retrieved  as  a function  of  cooling 
water  injection  temperature  by  a call  to  the  library 
supplied  interpolation  subroutine,  INTRPL. 


TSATFN 

This  function  subroutine  calculated  the  saturation 
temperature  as  a function  of  saturation  pressure  by  solving 
equation  (B-1)  for  T^ . 


NtATFAC 

This  subroutine  retrieved  the  material  correction 
factor,  F,,  as  a function  of  wall  thickness  and  tube  material. 
The  data  presented  in  Figure  ST  - 1 of  reference  [1]  was 
in  tabular  form  and  was  retrieved  using  the  library  supplied 
interpolation  subroutine,  INTRPL. 


APPENDIX  C 


USER'S  MANUAL  FOR  OPCODE  1 

This  Appendix  describes  the  data  cards  required  for 
the  use  of  OPCODEl.  A complete  optimization  run  is  included 
as  Appendix  D. 

The  data  is  divided  into  the  COPES/CONMIN  program  section 
and  the  HEI/DDS-based  condenser  design  program  section. 

The  COPES  data  is  segmented  into  "blocks"  for  convenience, 
.-^ll  formats  are  alphanumeric  for  TITLE,  END,  and  STOP  cards, 
FIO  for  real  data  and  110  for  integer  data.  Comment  cards 
may  be  inserted  anywhere  in  the  data  stack  prior  to  the  END 
card  and  are  identified  by  a dollar  sign  ($)  in  column  1. 

The  COPES  data  stack  must  terminate  with  an  end  card 
containing  the  word  "END"  in  columns  1-3. 

The  analysis  data  is  also  segmented  into  blocks  for 
convenience  and  must  immediately  follow  the  "END"  card. 

No  comment  cards  are  permitted  and  the  analysis  data  stack 
must  terminate  with  the  word  "STOP"  in  columns  1-4. 

It  should  be  noted  that  only  the  information  for  using 
OPCODEl  for  either  a single  analysis  or  for  optimization 
is  included  in  this  User's  Manual.  Information  pertaining 
to  the  use  of  the  sensitivity  analysis  and  the  two-variable 
function  space  features  of  COPES  can  be  found  in  reference  [15 


114 


DATA  BLOCK 


A 


DESCRIPTION:  COPES  Title  Card 


FORMAT:  20A4 


TITLE  CARD 


REMARKS 

1)  Program  is  terminated  by  the  word  'STOP'  in  columns  1-4. 


US 


DATA  BLOCK 


B 


DESCRIPTION : COPES  Program  Control  Parameters 
FORMAT:  7110 

1 2 3 4 5 6 7 8 

NCALC  NDV  IPNPUT 


FIELD 

1 NCALC : 

0 - 

1 - 

2 

2 NDV: 

5 IPNPUT: 

0 - 

1 - 
2 - 

REMARKS 

ri  Field  1 determines  program  execution. 

2)  Fields  3,  4,  6,  7,  and  8 to  be  left  blank  for  the 
OPCODEl  application  of  COPES/CONMIN . 

t 

A 

I 
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CONTENTS 

Calculation  control 
Read  input  and  stop.  Data  of  blocks 
A-B  is  required.  Remaining  data  is 
optional . 

One  cycle  through  program.  Data  of 
blocks  A-B  is  required.  Remaining  data 
is  optional. 

Optimization.  Data  of  blocks  A-I  is 
required.  Remaining  data  is  optional. 

Number  of  independent  design  variables 
in  optimization  or  optimum  sensitivity 
s tudy . 

Input  print  control 

Print  card  images  plus  formated  print 
of  input. 

Formated  print  of  input  only. 

No  print  of  input. 


1 


DATA  BLOCK 


C 


DESCRIPTION : COPES  Integer  Optimization  Control  Parameters 
FORMAT:  8110 

j 

1 2 3 4 5 6 7 8 

I PRINT  UMAX  ICNDIR  NSCAL  ITRM  LINORJ  NACMXl  NFDG 


FIELD  CONTENTS 


1 


2 

3 

4 


5 


6 


7 


IPRINT: 

0 

1 

2 

3 

4 

5 


ITMAX: 
ICNDIR: 
NSCAL : 


ITRM: 


LINOBJ: 


NACMXl : 


Print  control  used  in  optimization  program, 
CONMIN. 

No  print  during  optimization. 

Print  initial  and  final  optimization 
information . 

Print  above  plus  function  value  and 
design  variable  values  at  each  iteration. 
Print  above  plus  constraint  values, 
direction  vector  and  move  parameter 
at  each  iteration. 

Print  above  plust  gradient  information. 
Print  above  plus  each  proposed  design 
vector,  objective  function  and  constraints 
during  the  one-dimensional  search. 

Maximum  number  of  optimization  iterations 
allowed.  DEFAULT  = 20. 

Conjugate  direction  restart  parameter. 
DEFAULT  = NDV+1. 

Scaling  parameter.  GT.O  - Scale  design 
variables  to  order  of  magnitude  one 
every  NSCAL  iterations.  LT.O  - Scale 
design  variables  according  to  scaling 
values  input.  DEFAULT  = No  scaling. 

Number  of  subsequent  iterations  which 
must  satisfy  relative  or  absolute 
convergence  criterion  before  optimization 
process  is  terminated.  DEFAULT  * 3. 

Linear  objective  function  identifier. 

If  the  optimization  objective  is  known 
to  be  a linear  function  of  the  design 
variables,  set  LINOBJ  =»  1. 

DEFAULT  * Non-Linear. 

One  plus  the  maximum  number  of  active 
constraints  anticipated. 

DEFAULT  » NDV+2. 
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rj 


DATA  BLOCK 


C (Continued) 


FIELD  CONTENTS 


8 


NFDG:  Finite  difference  gradient  identifier. 

0 - All  gradient  information  is  computed 

by  finite  difference. 

1 - Gradient  of  objective  is  computed 

analytically.  Gradients  of  constraints 
are  computed  by  finite  difference. 

2 - All  gradient  information  is  computed 

analytically. 


RE.vtARKS 

1)  The  value  of  NSCAL  = 5 is  suggested  and 
ITRM  = NACMXl  = 0 should  be  used. 

) The  value  of  IPRINT  may  be  reduced  when  the  user  is 
familiar  with  the  optimization  output. 
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DATA  BLOCK 


D 


DESCRIPTION : COPES  Floating  Point  Optimization  Program 

Parameters 


FORMAT:  8F10 


_1 2 3 4 5 6 7 8_ 

FDCH  FDCHM  CT  CTMIN  CTL  CTLMIN  THETA  PHI 


Note:  Two  cards  of  data  are  read  here. 


FIELD  CONTENTS 


1 

FDCH: 

Relative  change  in  design  variables  in 
calculating  finite  difference  gradients. 
DEF.AULT  = 0.01 

2 

FDCHM: 

Minimum  absolute  step  in  finite  difference 
gradient  calculations. 

DEFAULT  = 0.001. 

3 

CT: 

Constraint  thickness  parameter. 

DEFAULT  = -0.1. 

4 

CTMIN: 

Minimum  absolute  value  of  CT  considered 
in  the  optimization  process. 

DEF.^ULT  = 0.00  4. 

5 

CTL: 

Constraint  thickness  parameter  for 
linear  and  side  constraints. 

DEFAULT  = -0.01. 

6 

CTLMIN: 

Minimum  absolute  value  of  CTL  considered 
in  the  optimization  process. 

DEFAULT  =0.001. 

7 

THETA: 

.Mean  value  of  push-off  factor  in  the 
method  of  feasible  directions. 

DEF.AULT  = 1.0. 

8 

PHI : 

Participation  coefficient,  used  if  one 
or  more  constraints  are  violated. 

DEFAULT  = 5.0. 
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FIELD  CONTENTS 

1 DELFUN:  Minimum  relative  change  in  objective 

function  to  indicate  convergence  of 
optimization  process. 

DEFAULT  = 0.001. 

2 DABFUN:  Minimum  absolute  change  in  objective 

’ function  to  indicate  convergence  of 

, the  optimization  process. 

DEF.AULT  = 0.001  times  the  initial 
objective  value. 


120 


DATA  BLOCK 


E 


DESCRIPTION : Total  Number  of  Design  Variables,  Design 

Objective  Identification  and  Sign  on  Design 
Ob j ect ive . 


FORMAT : 2110,  FIO 


1 NDVTOT ; Total  number  of  variables  linked  to  the 

design  variables.  NDVTOT  must  be 
greater  than  or  equal  to  NDV.  This 
option  allows  two  or  more  parameters 
to  be  assigned  to  a single  design 
variable.  The  value  of  each  parameter 
is  the  value  of  the  design  variable 
times  a multiplier  which  may  be 
different  for  each  parameter. 

DEFAULT  = NDV. 

2 lOBJ:  Global  variable  number  associated  with 

objective  function  in  optimization 
or  optimum  sensitivity  analysis. 

3 SGNOPT : Sign  used  on  objective  of  optimization 

to  identify  whether  function  is  to  be 
maximized  or  minimized.  +1.0  indicates 
maximization.  -1.0  indicates  minimization. 
DEFAULT  = -1.0. 
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DATA  BLOCK 


F 


DESCRIPTION : Design  variable  bounds,  initial  values  and 

scaling  factors. 


FORMAT:  4F10 


_J 2 3 4 5 6 7 8 

VLB  VUB  X SCAL 


Note : Read  one  card  for  each  of  the  .NDV  independent  design 

variables . 


FIELD 

CONTE.NTS 

1 

VLB: 

Lower  bound  on  the  design  variable. 

2 

VUB: 

Upper  bound  on  the  design  variable. 

3 

X: 

Initial  value  of  the  design  variable. 

If  X is  non-zero,  this  will  supercede 
the  value  initialized  bv  subroutine 
ANALIZ. 

4 

SCAL : 

Design  variable  scale  factor.  Not  used 
if  .NSCAL.GE.O  in  Block  C. 
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DATA  BLOCK  G 

DESCRIPTION : Design  Variable  Identification 

FORMAT : 2110,  FIO 


1 2 3 4 5 6 7 8 


Note : Read  one  card  for  each  of  the  NDVTOT  Design  Variables. 


FIELD  CONTENTS 

1 NDSGN;  Design  variable  number  associated  with 

the  variable. 

2 IDSGN:  Global  variable  number  associated  with 

the  variable. 

3 .AMULT ; Constant  multiplier  on  the  variable. 

The  value  of  the  variable  will  be  the 
value  of  the  design  variable,  NDSGN 
times  AMULT. 

DEFAULT  = 1.0. 


DATA  BLOCK  H 

DESCRIPTION : Number  of  constrained  parameters. 

FORM^:  no 

1 2 3 4 5 6 7 8 

NCONS 


FIELD  CONTENTS 

1 NCONS:  Number  of  constraint  sets  in  the 

optimization  problem. 


i REMARKS 

1)  If  two  or  more  adjacent  parameters  in  the  Global  common 
block  have  the  same  limits  imposed,  these  are  part  of 
the  same  constraint  set. 
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DATA  BLOCK 

DESCRIPTION : Constraint  Identification  and  Bounds. 

FORMAT:  3110 


Note : Read  two  cards  for  each  of  the  NCONS  constraint  sets. 

FIELD  CONTENTS 


1 

ICON; 

First  Global  number  corresponding  to 
the  constraint  set. 

2 

JCON: 

Last  Global  number  corresponding  to 
the  constraint  set. 

DEFAULT  = ICON. 

3 

ICON ; 

Linear  constraint  identifier  for  this 
set  of  constrained  variables. 

LCON  = 1 indicates  linear  constraints 
DEFAULT  = 0 = Nonlinear  constraint. 

DATA  BLOCK  (continued) 

FORMAT;  4F10 


1 

BL: 

Lower  bound  on  the  constrained  variables 
Value  less  than  -l.OE+15  is  assumed 
unbounded . 

T 

SCALl  : 

.Vormalizat  ion  factor  on  lower  bound. 
DEFAULT  = Max  of  ABS(BL),  0.1. 

5 

BU: 

Upper  bound  on  the  constrained  variables 
Value  greater  than  l.OE+15  is  assumed 
unbounded . 

4 

SCAL2  : 

Normalization  factor  on  upper  bound. 
DEFAULT  =«  .Max  of  ABS(BU),  0.1. 

REMARKS 

1)  The  normalization  factor  should  usually  be  defaulted. 
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DATA  BLOCK  P 

DESCRIPTION:  COPES  data  'END*  card. 
FORMAT:  3A1 


1 The  word  'END'  in  columns  1-3. 


REMARKS 

1)  This  card  must  appear  at  the  end  of  the  COPES  data. 

2)  This  ends  the  COPES  input  data. 
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OPCODEI  Analysis 


Data  for  the  condenser  analysis  follows  the  ' END ' 
card  in  the  COPES  data  deck.  If  the  general  design 
capability  of  COPES/CONMIN  is  not  needed,  the  condenser 
analysis  portion  of  OPCODEI  can  be  run  in  a stand-alone 
mode  by  using  the  following  main  program: 

C MAIN  PROGRAM  FOR  STAND-ALONE  CONDENSER 

ANALYSIS 

C INPUT  SECTION 

ICALC=»1 

CALL  ANAL  I Z (I  CALC) 

C EXECUTION  SECTION 

ICALC=2 

CALL  ANAL  I Z (I  CALC) 

C OUTPUT  SECTION 

ICALC=«3 

CALL  ANALIZ  (ICALC) 

STOP 

END 
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DATA  BLOCK  ^ 

DESCRIPTION : Condenser  Analysis  Title  Card 


FORMAT:  20A4 


DATA  BLOCK  ^ 

DESCRIPTION : Objective  Function  Weighting  Factor 


FORMAT;  El  2. 5 


DATA  BLOCK: 


CC 


DESCRIPTION : Condenser  Tube  Input  Data 

FORMAT : 6F10,  12 


_J; 2 3 4 5 6 7 8_ 

SDO  SDI  XW  SDD  ALGTH  ROUGH  NPASS 


FIELD 


CONTENTS 


1 

2 

3 

4 

5 

6 
3 


SDO:  Tube  outside  diameter;  inch. 

SDI:  Tube  inside  diameter;  inch. 

XW:  Tube  wall  thickness;  inch. 

SDD:  Tube  pitch  to  diameter  ratio. 

ALGTH:  Tube  length,  feet. 

ROUGH:  Tube  inside  absolute  roughness,  foot. 

NPASS:  Number  of  tube  passes. 


BLOCK 


DD 


DESCRIPTION : Sea  V/ater  Information 

FORMAT;  8F10. 


1 2 3 4 5 6 7 8 

VELC  T1  


FIELD  CONTENTS 

1 VELC;  Sea  water  velocity  in  the  condenser 

tubes,  feet  per  second. 

2 T1 ; Sea  water  injection  temperature,  °F. 
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BLOCK 


EE 


DESCRIPTION : Inlet  Steam  and  Heat  Load  Information 

FORMAT;  7F10 


_1 2 3 4 5 6 7 8 

WMS  HFGMS  WAS  HFGAS  PSAT  TSAT  TRNSHT 


FIELD  CONTENTS 


1 

WMS  : 

Incoming  main  steam  rate,  pound  per  hour 

-> 

HFGMS: 

Latent  heat  of  condensation  for  main 
steam,  BTU  per  pound. 

3 

WAS: 

Incoming  auxiliary  steam  rate,  pound 
per  hour 

4 

HFGAS: 

Latent  heat  of  condensation  for 
auxiliary  steam,  BTU  per  pound. 

5 

PSAT: 

Saturation  pressure  of  incoming  steam, 
pounds  per  square  inch  absolute. 

6 

TSAT: 

Saturation  temperature  of  incoming 
steam,  degrees  Fahrenheit. 

7 

TRNSHT : 

Heat  load  from  other  sources,  BTU 
per  hour. 

REMARKS 

1)  If  no  value  for  HFGMS  or  HFGAS  is  read  in,  these 
parameters  will  default  to  950  BTU  per  pound. 

2)  Either  PSAT  o_r  TSAP  is  to  be  specified. 


DATA  BLOCK  FF 


DESCRIPTION:  Tube  Material  Identification 


FORMA.T:  110 


1 IDNIATL:  Material  identification  from  the 

following  table. 


TUBE  MATERIAL  CODE 


ADMIRALTY  METAL  1 

ARSENICAL  COPPER  1 

ALUMINUM  1 

ALUMINUM  BRASS  2 

ALUMINUM  BRONZE  2 

MUNTZ  METAL  2 

90-10  CU-NI  3 

70-30  CU-NI  4 

COLD  ROLLED  LOW 
CARBON  STEEL  5 

STAINLESS  STEELS 

TYPE  410/430  6 

TYPE  304/316  7 

TYPE  329  8 

TITANIUM  9 
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DATA  BLOCK 

GG 

c- 

1 

DESCRIPTION: 

STOP  Card 

FORMAT:  4A1 

STOP  CARD 

STOP 

FIELD  CONTENTS 

1 The  word  'STOP'  in  columns  1-4. 


REMARKS 

1)  This  card  must  appear  at  the  end  of  the  analysis  data. 

2)  This  ends  the  analysis  data. 
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APPENDIX  D 


SAMPLE  OUTPUT  FROM  OPCODE  1 


ccccccc 

c 

c 

c 

c 

c 

ccccccc 


cccccoo 

C 0 

0 0 

0 0 

0 0 

0 0 

CCOCOOG 


ppppppp 
p p 

p p 

ppppppp 

p 

p 

p 


EEEEEEE 

E 

E 

EEEE 

E 

E 

EEEEEEE 


sssssss 

s 

s 

sssssss 

s 

s 

sssssss 


CCNTRCL  PROGRAM 
FOR 

ENGINEERING  SYNTHESIS 


TITLE 


CVA  / CASE  ONE  /OEJ=PGi.£R/MIN:  POU6R  WITH  CONSTANT  CPEJ/ 
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CARC  IMAGES  OF  CONTROL  OATA 


CARO 


IMAGE 


li 

i\ 

A) 

i) 

6) 

7) 

81 

IC) 

il! 

I :•) 
U) 
IS) 
lei 

17) 

18) 
IS) 

IV> 

ill 

1^1 

2!) 

281 

281 

I*** 

So) 

31) 

33  1 
3A) 

1 « I 

5)1 

381 


CVA  / CASE  LNL  /CuJ*(>C.»ER/MlN:  PCWER  i*lTM  CONSTANT  CRSJ/ 

1 AO  S 

18  -l.O 

S TUBE  C.O.  . INCFES 

c.tis  i.i« 

I TUBE  I.O..  INChES 
C.AC7  l.iCa 

S PITChi/OUMETER  RATIO 
1.1  3.C 

1 Tube  lencim,  peet 
1C  2 5 C 

s'^s-A  velccIty  in  Tuaes,  fps 


3.0 


S.  C 


t TUBE  AALL  THICKNESS,  INCHES 
7 7 

C.022  0.109 

i CC.NOENSEfi  OIAMETER,  FEET 
S B 

1.0  10. 0 

i PINCH  POINT,  F 

12  12 

S.C  3S.0 

s s-»  temperature  rise,  F 

13  13 

5.0  20.0 

S RATIO  OF  TUBE  hCLE  AREA  TO  TUBE  SHEET  AREA  WITHOUT  ORILLEC  HOLES 
lA  IP 

C.I  0.3P 

ENO 


IS 
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CASE  ONE  /QBJ-PQWER/MlN: 


POWER  WITH  CONSTANT  CREJ/ 


CCATPCC  PARAMETERS: 

CALCLLATtCN  CCNTRCL.  NCALC  » 2 

NUMBER  OF  global  OfcSICN  VARIABLES,  NOV  « 5 

AUMcER  CF  SENSiTlVI’Y  VARIABLES,  NSV  • 0 

number  of  FLNCTlJNi  IK  TwO-SPACE,  N2VAR  « 0 

input  information  print  coos, 

SEKSniVlIy  PRINT  CCCE, 

TtaO-SPACe  PRINT  CODE, 

CEEUG  PRINT  CCOE, 


CALCULATICN  control,  NCALC 

value  meaning 

1 single  analysis 

2 OPTIKUATIcN 

2 sensitivity 

A T’.C-VARIAoLE  plnction  space 

GICBAL  variable  NUMEtR  CF  OSjECTIVE  - 18 

multiplier  (NEGATIVr  INOICATSS  MINIMIZATION!  - -O.IOOOE  01 


CONPIN  parameters  (IF  ZERO,  CONH  IN  DEFAULT  WILL  OVER-RIOE) 


IPfllNT 

1 

ITmax 

AO 

ICNCIR  NSCAL 

0 5 

ITRM 

0 

LINOBJ  NACNXl 

0 13 

NFDG 

0 

focm 

0.0 

FOCmn 

0.0 

CT 

0.0 

CTMIN 

0.0 

CTL 

0.0 

^TLHIN 

Yheta 

0.0 

PHI 

0.0 

OELFUN 

0.0 

OABFUN 

o.c 

• 

IP.NPUT  » 0 

IPSENS  » 0 

IP2VAH  ■ 0 

IPOBC  - 0 
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SUNCUE  GEOMETRY  INFCRHATION 

ROW  NR.  RAOtLS  (IN.J  NR.  CF  TUBES 
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t .JO 
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e.37 
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.02 

5 

9 . 9 fc 
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191. 

08 

21 

23.32 
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52 
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29.19 
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97 

23 

25.35 

157 

.91 

29 

2§.92 

162 

85 

25 

2t.  78 

166. 

29 
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21. o5 

173 

73 

*7 

28.52 

179 

17 

28 

29. 3£ 

169. 

61 

29 

30.25 

no 

» 0 5 

30 

31.11 

193 

50 

31 

3 1 .9  £ 

200 

99 

32 

32  .35 

206 

.38 

3 3 

33.71 

211 

82 

>9 

39. 5£ 

‘i7 

26 

35 

35  .9*1 

ih 

70 

3t 

6t  . J 1 

226. 

1 9 

37 

37.18 

233. 

.5£ 

38 

38  .09 

239 

03 

39 

3£.9l 

299. 

97 

90 

39.77 

29» 

.91 

91 

90  .6  9 

2 35 

35 

Hi 

91.31 

260. 

79 

•*3 

92.3  7 

266 

23 

99 

•«3.29 

113. 
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eUNCLE  CEOMETRt  INECPHATICN 


ACW 

NR. 

RADIUS  (IN.) 

NR.  OF  TUBES 

1 

6.00 

31.65 

1 

7.03 

37.09 

i 

£.06 

4’.  5 3 

•f  .39 

4^.57 

m 

10.13 

53.41 

6 

11.16 

56.05 

7 

ii.i9 

64.25 

8 

13.23 

69.74 

9 

14.25 

75.  18 

10 

15.28 

80.6  3 

11 

16.32 

86. 06 

12 

17.35 

51.  50 

13 

18.38 

5t.54 

14 

15. *.1 

102.  •'8 

15 

2 C .■»4 

107.  82 

16 

21  .-.7 

113.27 

17 

22.51 

118.71 

18 

23.54 

124.  15 

24.57 

125.59 

k 

25.60 

ls5 . 03 

?i 

k.o3 

140. 47 

2 7 .6  6 

145.91 

23 

28.  70 

151 .35 

24 

a5.72 

156. ec 

25 

30.76 

162.24 

31.75 

lo7.68 

27 

ii  .S3 

173.12 

<£ 

j3  .85 

178.50 

34.35 

184.00 

lo 

35.52 

1?5 .44 

31 

36.55 

154. 33 

32 

3 7.58 

200. 33 

33 

35  .01 

2C5. 77 

34 

40.34 

211.21 

35 

41  .03 

2 lo. 65 

'6 

42.1  1 

222.09 

j7 

43.14 

227.53 

Ss 

44.17 

232.57 

aS 

53.20 

238.41 

40 

46.2s 

243. 65 

41 

47.26 

249. 30 

42 

48.30 

2 34 . 7 4 

43 

45.  js 

260.  1 8 

44 

50.36 

365.62 

45 

5i.39 

271.00 

46 

52.42 

276. 5C 

47 

53.45 

2?'.  .94 

4C 

54.45 

287. 38 

49 

55.51 

J52.62 

£C 

5c. 53 

k8.27 

5i 

57.58 

iOs. 71 

52 

5 8.61 

309 . 13 

the 

OUTER 

RCli  HAS  A PITCH  OF 

1.19  INCHES 

146 


vPPEN'DIX  E 


OPCODEl  PROGRAM  LISTING 


OOOOOOOOOOOOC OOOUOOOOCOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOO  T/OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOUO'->0'-Jw^OO‘J<^<wOO'^OOOUOOOoOOOOOOOOOOOijOOOC>^OUO 

a.a.a.a.a.a.Q.a.aa.aQ.a.aaaa  Q.:^a^a.a.Q.CLa.o.a.a.Q.cLaQ.Q.o.acLQ.ao-a.aQ.a.a.Q.cLa. 

oooooooacoocciocoaccacoocucooooacoaacoooooooococo 


(/> 

< 

z 

a»" 

UJ 

* 

h“ 

^LJOO 

X 

iUt-Z  ^ <X 

*^/)  GL  3 


• 

1 

U Ql-4»-ioo 

CD 

u 

t 

0t/>O3t/)^Z 

••u 

uy 

1 

X » 

c^ 

1 

XXdX  UJ 

«» 

1 

U^x^zx< 

-^oO 

00 

1 

z 

QO 

O' 

1 

XXZ  — Z 

Zu 

(Ti 

1 

LJUCX(/)UJ01 

CQUJ 

ir» 

1 

Z^LJXXLJ 

••X 

OJ 

1 

^OH-  ao^x 

^ •> 

1 

A,*“» 

» 

1 

'^Oa.  «j 

••M 

o 

1 

u>- 

•• 

1 z 

OIZXXOOZU 

Q< 

1 

o«*-*xz<  • 

tOiX 

•M 

1 o 

CLa:^-XuJ  >z 

•*< 

rn 

1 

OXxcO 

o •• 

00 

1 

ooxctozxct^ 

-JcO 

<M 

1 

^ 30  X 

xuu 

• 

1 t— 

</>XZQOl-^^- 

^OJ 

o 

>o 

1 

•• 

1 a 

X x*-«x^xx 

x»- 

(J' 

1 

uy 

1 UJ 

zaox^as  z 

o< 

o 

^4 

1 

<^o  OIX  -o 

XO 

QQl 

<r 

1 00 

wO'S*X<CJ-i</) 

<x 

p--oa 

1 

XXl/)X^UO.^ 

QX 

• 

1 

IQLUX.  ox 

X • 

1 1— 

»-  Q.  uj<ia 

0<M 

o •» 

V. 

1 

QOZUJ 

.^H-O 

>-^ 

o 

1 3 

^XOO</) 

•o 

> 

1 

Qir> 

< 

1 a 

•^ujQ  ^aji<l  •■  • 

o o 

1 

i/iXXuj  •at—  •an-rD 

U-  X 

1 z 

UJ  P-.  ^ OOOO  M < O 

h-QO 

•M/) 

X < 

1 

aCH-  <lQlh-3 

«^i/x 

XX 

*0 

i 

4Q.u.CUUi<lO^ 

••X 

1 

XOO  t-tNI<0O(vtQ  .. 

Z)a-^>r\j  > 

1 

XUJ*-IQ^ 

^Q(X 

Q<-0®'-'-' 

1 

XUJZ<XXO  • 

X(/)UJ 

o^cc 

—Ht-O 

1 

<XOa£<i->H-^-<'>»3 

■v'viTta.zm 

1 

ot  t- "»  o QC »- </i  r- Z X o 

^('J'V  •*U  •• 

1 

a >— a*iu.au><4act 

^ Jtt— |-1(JCJ 

1 

aix<ixQ.oa.r-t 

CO  . 

ooaQ.i-«<\j 

1 

o^aua 

LUO-I 

ouoa 

1 

Z_lO 

_J_H-'»-ZO 

1 

® ^(/)uj4<a: 

— o> 

2Q(S 

• •O^UJ 

t-0<a.o0  2HH<aD32X<Z2:^-^Q.oO  z 
aCOUJoiZXOOh-XOO 
(/)o  «joui  oja^XQ.  •>xx<<aiH-H* 
^UUJ<ZXUJ^COXO<^*XX»-H-a.  z 
X ac  X ^ l;  X ^ oj  ju  o «.j  ^ Lj  ^ 
h- a.  H- < o h- UJ  t/)  o JB  ^ o W O Q o o o 


a 

> 

< 

o 

X 

X 

I- 

il 

X 

GO 

/> 

< 

z 

GD 

1-0 

01 

u. 

1-0 

X 

<l 

•» 

H- 

•-* 

z 

<. 

o 

l/> 

• 

X 

1- 

o 

3 

o 

< 

Oh- 

OC 

GO 

>o 

• 

a. 

^ • 

X 

m 

Cl 

>— 

GO 

xco 

Cl 

x<a. 

X 

CJ 

II  3 

rg 

< 

X 

GO 

X 

X 

z • 

» 

o 

»» 

CIO 

X 

Q 

GO 

xz 

31- 

GO 

< 

x< 

— • 

X'J 

•• 

3 

• 

3 

3 

•■  a 

X 

«< 

X 

GO  h- 

oo 

•» 

Ci  ► 

» 

X 

• • 

O H- 

UOX 

— <0  a 

H-  Z 

-o 

Q»— X XO 

<^X 

U3 

f— 

t/) 

-Xt- 

• • 

Xl-O 

uo 

z 

ujx 

GO  01 

UJ 

O-fooZO— 

— • 

«G0 

Q 

QUJZQ  •(- 

I— 

• • 

H— 

<(J0>3  — O< 

< 

(M— OO^-hCM® 

GO 

(/) 

u.uJOrr><noOO 

X 

Goto 

*— 'VOvO'O'O'O'O'O 

oooooooaz 

z 

oo 

OXH-h-l- 

•uM/>tnir\uMn  •x 

X 

LLu-t/Kioaxa 

333 

rvj  • 

» •» 

•>  » «•  h" 

h— 

xx<u.xooaooa 

-^zzzzzzi— «iru<uj 

Iw'-pJKh- 

P-H^P-4 

1 II  — 

^ VI 3 oOl/)  p 

IQ— 

m'>«.»wGOX 

t_2CQcDt03:  II 

Q.  11  a II  '-'<<z 

OXXXXUJX 

IQLJQOOOU^^I 

«IZ 

II 

^ (J  Ui  uj  UJ  UJ  UJ  UJ  Up  I O ^ U.  X V * ^ OL  OC  OC  01 QC  oc 
X ^ OL  q:  (X  QC  01  QC  Qp  U>  r- O ^ W UO  ^ -B  :x  i -B  i 


uuuuuowo 


o 


o 


147 


r 


I 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOCJOOOOCJOOtaOOOOOOOOO 
O' o 'O  O •-•C'J  ni  >}■  IT*  vO  r- CD  Qv  O If* -O  00  O' c ir\  <0  ® CT*  O ^ f ' "S’ If' >0 

•4'ir\u'uMr<u'Mr'U'uMr\u'  >w><j%4j'0'0'ti'0'U'W  ^r-r-r^r-r-r-p^r-r^i^co®  ® ® w®  ® ® ® ® O'O'  y'O'O'O'O' 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOtJO 

U O O O O <_)  O <-J  O l_>  O (J  O O O <_)  O O O O O O O O O O O O O O O 'J  O O O <->  UO  O <-)  O o O O O (JO  o 
a.a.a.a.QLa.a.a.a.a-a.(XQ.a.a.a.Q.a.a.a.Q.Q.a-Q-a.Q.a.Q.a.a.a.a.c.a.o-a.a.a.Q.Q.a.Q.a.ixa.Q.a.Q. 

oooococuQOOcaaoaooooaooooooccoQOoooooooQcoocoooo 


1 

1 

1 

1 

1 

1 

1 

\ 

1 

1 

1 

1 

1 

( 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

t 

U 

1 

1 

1 

t— 

1 

1 

o 

1 

1 

►- 

1 

OJ 

1 

X. 

t 

1 

00 

! 

1 

z. 

1 

1 

a 

cc 

1 

1 

M 

-(O 

< 

h- 

1 

CO 

1 

1 

\ 

— t 

X 

OOOT** 

a 

(/) 

1 

M 

1 

-4  00 

< jcm 

oox 

1 

1 

_l< 

^Z 

JO  • 

-SO>- 

1 

OO 

1 

M LUU. 

n X 

u-QCO 

oO 

•Ki.< 

1 

1 

Ojh- 

oO*-^Z  ^ 

u 

^Xoo 

1 

> 

1 

X » Q. 

< > 

LU  *•  II  • 

< 

£ OIL 

1 

1 

••  •X 

U-oO 

LJ  OUJl/VCM 

a. 

LU 

1 

-J 

1 

sO 

^ OLU 

'SX 

> 

1 

1 

vf 

1-  ox 

-J-O 

UJ 

«> 

«> 

®o< 

1 

< 

1 

• 

< 

uu. 

> Q-fraa:— 

•iM 

<0.^ 

1 

1 

00 

X ^h- 

<IX 

N*  1— (MOO^OO 

u. 

J6X^-  < 

1 

1 

'I- 

LJ 

Oli-  » — 

X -S 

1 

t 

'!■ 

MM  MCO 

• '-'oorsj 

X — -1* 

1 

< 

1 

O » 

u •■  U.I— 

(MOXO 

•«■  0<Qi^«-*QCC 

0'®vfif\®r~(r>c'-iNt(r( 

1 

1 

• rg 

_if\j  «•  1 

xC3< 

(j-m-coc  z< 

000000(Mro®o<^ 

1 

1 

(M  ^ 

UJCL  X«-4 

— -J4.U- 

*LU<i— O •*—(:£ 

0®'0'0'0>0'0'0'0'0 

1 

t 

O'v  ♦ 

>—*01-1-  b 

H-<loO* 

X>UJ*iX  a(>0QC 

* • 

1 

1 

Q—  (_) 

M *»  uLvr^uj;3^<^ 

h» 

1 

1 

• tOOM  M 

o»— ro-jMH-h-H- 0^0  •>r<-04-<OM  •M;aa 

:J3003r>3333 

1 

1 

(Mr'J*Qi-0  K- 

cc<'-a.<  + 1 1 

1 0(NJU.O 

•0<JI-®  Q*C 

oaoaoooooo 

1 

1 

* 

OU.ij'IHH-H-H*H--fOXO®<I— I-*XO<  ‘Z 

MM  1-41^1^1-4  UJ 

>*^Q  t iCL'T 

i/>H-<irsio<-o<<^<-fr  xo<QC(Ma5 

w 

0'—(/)0  ♦o 

* <U.U.0(/)00>'^l— u.(/1  — 

'T  II  OiJ  II  UJ_j  II  II  II 

zz 

QQ  II  Uo0»-® 

•XX*''I-|-|- 

1— w II  <xx<  II  j(<0(ju.ozi 

UJ  LU  UJ  ui  UJ  UJ  LU  UJ  LLI UJ  QC 

(/!(/!  11/1 X * 

oo  uj<-«  II  II  II  II 

II  II  Qizoovo  II  UJ  -^<l<lh- 

»— H-»— h-3>- 

II  II  O— i — CS.v^r-OI-U.^I-rgfn 

roi-  II  — 

II  QCLU^^^mQCmO. 

M 

C9U»~  II  II  a.  II  It  II  u)^<^ijo^xz^^xajuju-j  c^cjlj 

(Xcr(xcx(xacix(xar.(xuji^ 

II 

M II 

^M  t| 

u.h-^cxx<iu.^^a» 

H*  a.  ^ o Cl- ^ o u.  3 u.  ix 

LU. »- <i  c-#o  o X h- <3  < uj  cn  ca  ^ 

a 

ru 

OOOOU 


148 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
00  CT' o ^ 'O  QO  O' O fvj  ro  tn  sO  oo  (jN  o fNJ  ro  m 'O  GO  O'  o ^ >0  00  O' O 

O' O' tr  o O o o o o o o o o ^ ^ ^ ^ <M  <N*  ^ <NJ  csi  (\j  cN  fM  ro  n 1 m n > r ) m i t > fti  m >r  >j- vj- ^ 


VJOOOl->OOOOOOOOOUOOO^OVJUOOOOOOOOOOOOOOOOUOOOOUOOUO 
Q.a.aQ.Q.a.Q.aNQ.a.Q.Q.o.a,a.Q.CL^a.Q>CLa.Q.o.Q.Q.Q.Q.a.Q.Q.a.a.a.a.Q.Q.Q.a.a.a.Q.a.<xo.Q.a.a. 
OOOOGQOOOOOOOOOOOOCOCICOUOOOOOC  000000000000000000 


# * 

« ^ 

♦ ♦ 

♦ 

it-  ♦ 

♦ ♦ 

♦ 2 •» 

♦ O ♦ 

<— li 

♦ *-to>  ■» 

♦ 1-^-  •» 

1 

1 

o 

•n-  u*-  •■»■ 

1 

1 

u 

« Z^UJ« 

1 

1 

OJ 

♦ 

1 

1 

•»  a-OZ» 

1 

1 

to 

-to 

■» 

I 

1 

2 

(NJ 

« UULUlX« 

1 

1 

N-» 

« > LU<ff> 

1 

I 

CO 

♦ -•zx* 

1 

1 

Jt 

LJ 

1 

1 

to 

•-to 

n-u;  oc.* 

1 

1 

QL 

« LUU. 

1 

1 

UJ 

> 

L-) 

“7*— fco.* 

1 

1 

z 

< 

♦ OJ  u* 

1 

2 1 

>JJ 

h- 

OJ 

♦ u ••  ♦ 

1 

1 

to 

w 

* 

* <o* 

1 

O 1 

o 

UJ 

o 

* 2 Ztf 

1 

1 

00 

• 

■It-  < -(—.It 

1 

to—  1 

LU 

2 

•»  (X  -It 

1 

1 

> 

LU 

+ 

■>*•  i/lCoi'lt 

1 

►-  1 

X 

to 

o 

f\J 

^ <>-UJ-»- 

1 

1 

1-0 

to 

♦ 

♦ 

1 

O 1 

1— 

Q.-to 

tO-T 

♦ 

♦ 

1 

1 

X 

OJ^to 

in 

<l  • 

X 

« UJtU< 

1 

LU  1 

to 

aj< 

rs. 

a.>r 

LU 

♦ oO  X-fr 

1 

1 

2 

JLCL 

• 

20 

> 

0. 

» OOU-lt 

1 

y)  1 

rvj 

a 

to< 

«T> 

It* 

Q 

♦ 2 « 

1 

1 

LL 

< to 

II 

OfM 

rsi 

312 

♦ O^LU* 

1 

1 

to 

to 

LUO 

X 

♦ 

ZX 

« i^h-UJ« 

1 

h-  1 

OO 

to 

xo 

K» 

I— 

■»  — Xi-o.)t 

1 

1 

zx 

l/)X 

<uO 

a 

-o 

o 

» o 

OJO 

1 

3 1 

X33 

SOh- 

y)  to 

Q 

00  -J 

Mito 

X_l 

* B—  » 

1 

1 

»o 

•-X  < 

i-< 

JO 

UJUJ-T 

h-< 

♦ LU  •■if 

1 

0.  1 

o 1/1  Xi/) 

-o 

i— 

CQ>  • 

o* 

a* 

It  i>ozt--» 

1 

1 

Q •■ 

toJ  31  to^ 

XX 

Onj 

o*  -r 

-j>r 

Qf'i-J 

♦ to-i  OJ-lf 

1 

H-  1 

h“ 

tOi 

LU  .*tO  to 

<u 

•i^  • 

a3> 

CD-fr  U 

* <r  * 

1 

1 

2 

> *•<!— 

— to 

ir><\i*  < 

-O'*. 

xa3 

z»  > 

♦ QC  * 

1 

O 1 

LU 

OO 

to  O00< 

f'-*  a.UJ 

— 

•»  X X 

# UJUJO-fr 

1 

u 

O-J 

^ c^to 

•i<Qt 

UJH- 

lO(X 

~<x 

# 

1 

U 1 

— 

ux 

1-  Xil- 

03^ 

o* 

■«■ 

1 

1 

Q.Q</) 

a.— <iof-_iQ. 

• 0-1 

H-  ^ ^UJ« 

1 

1 

*to-to<«-»toto 

-to -to -to 

-to -to -to 

• QZt- 

a.LLl/lh--lOC-J 

^20 

« <003 

CL 

1 

1 

OOfNJ-^ 

►_a3£t)^ 

toxh-Oix 

>^ia> 

* iz  * 

LU 

1 

1 

Oi'OfOO 

of^OOoOOiM— 

OJC*  </) 

Z — >-_l> 

•—I*  Q 

* Z»-</)-lt 

X 

1 

1 

>0<3O< 

'jj  'JJ  > •j'l  <4j 

• -^nTUJ 

♦ ^ ir 

^ ir  OCX  >rx 

♦ -'♦ 

OJ 

1 

1 

to  to  to  * 

to  to  « 

- to  to  to 

xx®ooa<ac*  azo^uiT^*  co  + 

*1-1— ►-■» 

« 

1 

1 

b b b b b b b b b o b b b b b b 

h-h--*x  OJ>a.QU-  z 

^-ICO'N.O*— -i 

* 2* 

< 

1 

1 

oooooouooooooooo 

0_j  II  II  i/ioontoc  II  II 

ZOH-OJIM  II  > 

# -7UJ<  ♦ 

-1  o 

£Q<<Xa  Z II  XQlI/1</1  II  I—  II  Z II  _ia3 

♦ (Xcnix* 

022 

^ II  UJM 

II  _j  OOi/i  II  acO_iOZ 

>(X-' 

UJLUUJUJiiJUJUJUJUJbUUJUJUJlUUJUJ 

u 4 ^ a.  ^ -J  ^CZ3CJ  ^ ^ 2 

II  >CJ-J 

* 1— xz* 

tototototo—M 

toHtotoi-M 

tototototoHMtototo- 

u. < QC  ^ <4  ^ ^ ^ Uj  OJ  ^ O 

W UJUJWir 

■UUJOJ 

(xacix.a:crixccaiuL(XQLOcacQCcxcx. 

O ^ < J(  X ^ U H- UJ  JB  o Q.  X 00  ^ 

♦ OOQCVJ* 

UO-O 

xxxxxxxxxxxxxxxx 

o 

(T\ 

OOOOOOO  ooouo 


149 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOQOOOOOOOO 

in  >o  CO  O' O cvj  t’O  irv  00  o>  o ^ ® O p-N  rsi  rn 'T  iTN  ^ r*  CO  c^  o ^ >0  ^ CO  O' o CNJ 

>r  >**  u nn  u > u ^ m u ^ u Mn  m m >o 'O  >o  >0 'O 'O  ^ 'O 'O n- ^ ^ r- 00  00  00  00  CD  00  00  CO  CO  00 1/' li' O' 


UOOUOOOOOOOOOOOOOOOOOUO'-^OOOOOOOOUJOUUI. 

a.  ^ (X  a.  ^ Q.  a.  Q.  a.  a.  a.  a.  Q>  a,  o.  ^ a.  a.  a.  a.  ^ o.  ^ Q.  Q.  a.  a.  a.  a.  a.  a.  a.  Q.  a.  a.  a.  o- Q.  a.  a.  a.  a.  a.  a.  o.  Q.  CL  a. 

coooocoQooooaooooocoococooocoooooauQUoocooooaooc 


• 1 

• 

x-i 

• 

♦ ♦ 

ft 

• 1 

1 

03 

X 

• 

♦ # 

ft 

I 

^ * i 

1 

X 

Z3 

uO 

♦ ♦ 

V 

• 

H* 

• 1 

X 

Ujh- 

» 

X 

* * 

ft 

0 

1 

.» 

•> 

X 

X 

CO 

If  ♦ 

1 

• 

1 1 

X 

X 

»— 

♦ ♦ 

1 

X 

1 1 

X 

1 

< 

0 

< 

* * 

1 

XX 

X 

1 ••  1 

X 

X 

XX 

X 

# ♦ 

1 

XX 

1 0 1 

h- 

^ J) 

X 

♦ ♦ 

1 

XX 

X 

1 *-  1 

<& 

<✓1 

«'<  z 

X 

X 

If  ♦ 

1 

z</> 

X 

1 1 

a 

X 

If  * 

1 

1 X 1 

z 

zx 

X 

♦ ♦ 

1 

X 

a 

1 X 1 

<1 

X 

• 

♦ * 

1 

xo 

X 

1 X 1 

Ui 

< 

'dX 

X 

♦ ♦ 

1 

OJ  J 

> 

1 ►-  1 

UO 

X 

X 

X 'T 

* ♦ 

1 

X 

1 X 1 

X 

.fl 

X 

X 

X 

♦ ♦ 

1 

Oil 

X 

1 0 X 

0 

nj 

0 

0 « 

u 

0 

♦ ♦ 

1 

•0 

uin 

1 OJ  r-4 

rsi 

» 

(M 

oum 

fNJ 

iM 

♦ ♦ 

1 

■Oro 

ruf- 

1 p>  ru 

» 

m 

•a* 

«. 

ft 

» » 

1 

X » 

ft  •> 

1 0 

0 

h- 

0 

0 

0 

« « 

1 

•0 

O''. 

t >r— «> 

>r 

>r 

— >r  ► 

sT-^ 

>r 

♦ ♦ 

1 

>r 

^>r 

XI— >r»— 

>s 

r—  fMh— 

♦ ♦ 

1 

x^fM^>r 

» 

«» 

*>  « 

• 

4f  ♦ 

X ^ 

• ft  • 

f-I'TO's. 

0 

z 

00*^0 

r^onj 

* « 

00 

OJ 

0^0 

• 

u 

• 

^4  • 

• 

— *«■  ♦ 

X • 

P-4 

lOOXZ 

0 

« 

Ox  ox 

ou.  oo»  * 

» 

xoxox 

^ *0 

K- 

ft 

P-*  ♦ 

•X 

X 

< 

X 

• X 

X 

X 

— ♦ ♦ 

►- 

^x 

• X • 

» 

z 

• »x 

•a 

• 

•-if  # 

ft 

X 

• ft  • 

► • f<t 

• 

X 

• 

a X 

•a 

f 

•*  » 

s. 

X • • 

Z • ‘S 

• 

a 

• 

X >N^ 

•X 

• 

•*lf  ♦ 

ft 

<1^->1-  • • 

'.J  • tx 

X 

X 

• 

"vU-O 

H-OOC 

— X 

• • 

u 

oc. 

OJ 


gLOj 

jj 

X 


1^  <1^  •«■  «■ 
U-r-  ♦ 


*>  ta. 

M t<i  ••i  •< 


0 

X ^ 

X 

X 

33  m-t-  . V/l  .<H-  !»•  * 

<U.  OJO-  *0 

u 

t/i 

X 

oc  • 

X 

JZZ*-  i/IOCa-Z*-# 

Z •■  Xf—uj  r>>> 

I— 

XX 

ft 

0 • ft^ 

X 

X 

.QUJO.  UJ(X  OJ  » * 

acDOuja«  X 

X 

XX 

x<^xx 

H- 

f- 

•<u  UJ  ZU  *0*  ♦ 

LjniotZ  *0 

X 

fMti 

INI 

r^xi- 

cO 

Zl-ZX  XuZ  — ♦» 

< 

me 

X 

0 

-ftft-  << 

z 

<i/1<  0 < « # 

Z •Qti—Ou.— 

ft 

X •* 

X 

cO 

<X 

0 

hft 

'j-l  JJ  • *?  •UJ_l«  » 

-HCJlilOlC  ^ 

X 

<3C 

ft 

3 

•xO 

ftft 

X 

1- Jl— -i— 0 J_i<*  ♦ 

•►—  ac  » 

0 

XX 

< 

ftft 

OUJ 

.ift 

>f  £Ox< 

>— 

X 

y3<00»  :OIX.  # 

luOXluujujlu 

ft 

X 

0 

♦-X 

fU 

<ix  :ox 

X 

z 

>—  t_)_J  t-  ♦ 

-ju.'-i-jr-a.oa 

rvjO 

X 

< 

< 

O^xOi/) 

X 

•0  •O'S  •<U:UJ*» 

a « 00^3 

X 

XX 

< 

OC 

oo. 

ft 

<MXa3>^  X 

-) 

X 

X^-X<3>/>Z3^^-^»  ♦ 

ZfNiUJZ3—H— 

oc. 

< ^ 

ft 

Xs/1 

••XOX-* 

X 

n 

3IDXt-0QX3<** 

x<Mi3iax 

X 

xo 

X 

ftft 

X 

• 

»nft^f-xrsi 

fM 

a^-^-anjxj 

>^ 

«*v  — 

0 

^ ftXi^  • 

X 

Q 

t/'V 

XrMXcM  XfNiXX** 

X -XXfMXfM 

or^  iTs  _ 

00  O' .-«<Mftl  0 3 P.J 

x*^ 

tn 

X 

CM 

U t 
►- 

ij  ft<J  ^ ^VJ  ft\.^\^  tC  W 

(Njinrumxnjin'NJfMX  # 

(M  fM  lf\ <Ni  ITi 

• 

^ •• 

• ft/^  ft 

•ftft  * 

• ftO* 

0 '0>o>o>o>0'Oa  >0 

O'O 

X » 

0^f 

-ift  • •«* 

in^N^u^f— 

lAN^'s. 

inh-mf-zmH-^ninm'*  # inuMninH-vnh- 

ft^  H-  •» 

"•^GOOCa  • • ••■sT  » 

*0 

OU.U.  (-lUJ  X X XO  X 
<Mooma3«4}»-4^ 


•nj 

•O 

xox 


0000000  II  o 00 

^ ^ ^ ^ ^ QJ  ^ yj 

a«aJUiujujujijjOii4'^'^uj^QC<< ^ << <<<  •<  I <</)< 

•t  i zclx 

m^^P~ccaL<x(x,a.cca:cc  •ol  » ocilo: 
(xacuc(xuLac.Qc.uicx,u*^ocLjujooi-JOi-iLjuuj  •u/  1 u 
j:Lj^(^>-«^oa.LLiLa.a.iLu.u.ix  *0.111 


• V*  *0  X # » » • •►Cnj 

• *0  • •^♦♦O  O • • 

I xoxoxxoxox# ♦ xxxxoxo 

^ ♦ ♦ #H^<i  I 

i ^ k 

I < •<  0<  •<  •<»*<<<<rsi<  • 

tx  •x>>rx  .x  •x^^xxxx^fx  • 

I (X  •Qca^-a  •x  •x#  ♦xxax^-acH* 

LJ  I l-.*XLJ-J  •ui  •u  •0’»  ♦ UIUCJU^UO. 

»u.  i i-L xi^<(X^ix  *u.  •ii.'»  ♦ xa.xu.xu. 

H ^ ^ *^00  •M  • ^nj  ^ •-• 


<>4 

.H  ^ M t ' l >r  X ^ (^J 

m 

>!• 

ai 

iij 

X 

u\ 

a 

0 

00000000 

0 

0 

0 

0 

0 

0 

0 

pH 

.0 

>U 

-0 

•a 

0 

>0 

>0 

ISO 


OOOOOOOOOOOC  OOOOOOOOCOOOOOOCJOOOOOOQOOOOOOOOOCOOO 
<»> -f’ (f\ ,0  00  o «M ' ' ' ir>  >o  ® (7' o rsi  rf' if\ >o  CO  O' (j  t-t  ni  vj- If' vO  oo  O' o •-<  CM  ' nJ- ir>  uo  O' o 
O' O' O' O' O' O' CJ  O O O O O O o O O >-1  ^ CM  CN  • N«M  04  CM  CM  CM  CM  <Ni  cn  m c T»  rn  m fc  c cri  rc  1 o t 'T 

CM  CM  *M  CM  CM  CMCM  CM  CM  CM  <M  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  eg  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM 


O <sj  VO  W O O O c^  C_>  VO  O C_)  O O c^  O O CO  <->  O O O O C_l  V_>  O O O U O O O O c^  C_)  O O C_)  c_)  O C_)  c_J  (_)  c_4  O c^  <_)  c^ 
ci.a.a.a.aa.a.^a.a.a.a.a.a.a.a.a.a.a.Q.a.o-c3.Q-Q.Q.a.a.a.aLaa.a.Q.cxQ.caLaQ.<3.Q.aQ.a.a.a.Q.Q. 

oocooccaooooooooQooooocziooooacccooaoooocoooooocoQ 

I 


r*. 

• » 

• » 

1 

r<\ 

• 

(X  ex 

• UJ 

1 

• 

UJ  UJ 

® aC 

1 

• 

h" 

jC  ^ 

U4  .0 

1 

• 

• 

0 u 

Xi  oO 

1 

» 

• 

•— 

X a 

X ® 

1 

X 

• 

0 

UJ 

1 

CJ 

• 

LL 

(J  0 

QL 

1 

z 

z z 

X Q. 

1 

u 

X. 

P-4  ^ 

u 

t 

—a  X 

• 

» 

•»  1 

oe: 

H- 

M.X  X 

• ^ 

■v 

1 

0 

< 

A 

• Xl  J 

QC  < 

•*— II  1 

X 

a. 

vU 

<x 

(/IX  X 

Z c/i 

♦ ♦ 

«'  * 

It  # 

» ♦ 

0 

OL. 

> 

X 

UJX  X 

X X 

♦ 

•» 

♦ X 

< 

a 

C3 

cao  0 

0 0 

* ® 

♦ 0 

X 

u. 

X 

(M 

XCM  CM 

'M  ru 

♦ <M 

z 

0 

►—  » » 

^ » 

# 

it  m 

<J 

CM 

—0 

0 0 

0 0 

♦ 

it  0 

z 

* 

X -r--*  'T 

U) 

^ X 

♦ U\ 

CJ 

0 

Cj^-rj  ^ 

* 

♦ 

it  h- 

►— 

«»  • » 

♦ 

uu 

-»  X 

ir  ► 

X 

l-CMCI-  CMOtJ 

fViCJ^ 

# 

» 

# V* 

a 

a. 

• ♦ •^1— 

« 

• 

♦ 

3 

it  ► 

<—« 

X 

iTiOI-OU-t- 

ZOLLI-O 

ou-o 

« 

■» 

X 

•.-<U.~4  •< 

X*-4 

«P4 

-J 

« X 

♦ OJ 

3 

a.'vMu.*  a.  •X 

CMU.  •Xu. 

U.  'X 

fr 

♦ 0 

•‘CL  •‘  •QL 

^ • •» 

♦ 

< 

» ac 

* a 

z 

Q 

1 uu 

• X .c/)0 

» ► ♦ •X 

•c 

■» 

•fr 

♦ 0 

0 

Z 

j » 

• ^ •XXv^LTV  *> 

• ^x 

-a* 

> 

it 

•ip 

< 

1 • 

•0 

♦ 

it  <✓> 

it  ^ 

»— 

1 • 

•-JX  ^ 

XOCQUX 

•* 

♦ 

-fr 

it 

X 

a 

1 f 

1 •KJJO  •’QL 

• XUJ 

♦ 

H- 

it 

UJ 

1 • 

• OJ  •XJ  V 

1 -•  o<i 

h-  X 

♦ 

♦ 

it  QL 

•1 

0 

1 • 

• Xujt-ac: 

1 • V'UJCO 

UL  •X 

♦ 

X 

« 01 

It  UJ 

>s 

X 

1 a 

1 zujccf-oc 

XI- 

♦ 

♦ 

* h- 

» » ^ ^ 

iO>Z 

1 UJ 

• ^•-4  OCX 

1 — H-X  •“XI 

•mX*I 

X 

> 

it  < 

Q'v 

X 

•< 

1 ^ 

• “cX  X 

1 '-<01VJu. 

xax 

•» 

* X 

» «l 

•»»»*»  lUj  •»  » » 

oa 

1 0 

•0  xo 

1 Q^<oOUJ(/) 

H-OUJ 

♦ 

z 

* 

^^.-iCSjCUrvirOx^  jiA  ■xx^v.ox*— 1 

1 UJ 

• O'  •Otu 

1 (/>  UJ<-OZ 

0 Q- 

■fr 

♦ LU 

X ► ► p. 

X 

IL. 

1 “> 

•»-Q._JO 

1 

z •x 

♦ 

♦ < 

♦ CD 

• • 

•u 

1 LUO.  'Xu. 

1 — OQ.U.[XCMUJ_JUJ 

•» 

♦ 

» X 

• t 

• UJ 

•< 

1 IX 

a^uj  UJ 

1 Q-J  ®*-* 

X«  K 

# z 

* 1- 

•a: 

• 

• -> 

•Ui 

1 

• XH*  jK-4 

1 <x.<^  « 

♦ 

♦ 

♦ 

XUJ 

• •UJ  • • •QUJX  • • 

• 

• X 

• I 

1 1 a 

1 ac  1 ♦ ►-KUj<  • 

♦ 

it  <l 

* X 

<0. 

• CO^sl  • • •XUJUJ  • • 

• 

•0 

• X 

1 <H-XX'/1Z 

1 i^c/i(-<i  u-mxt— 

# 

•!»  0 

H-O. 

•</)Z  • • •OH* UJ  • • 

• 

•X 

• r^ 

1 UJ 

x>-  1 x;o 

1 <M  1 XU.lijX:3X< 

•fr  ♦ 

Osl 

LUO 

• <X  • • •OU^H’XJJ 

• 

• r- 

•m 

1 x-j®v/^voca 

1 ^uOOXXiTV^-Ot/) 

xixxxxxx  •• 

xvj 

• XX-J—* •—X  c/lOI-g 

• 

• «. 

• «• 

c 


I ZXIXc-cCX 
I nO'OO  •'O- 
I — *CM  •'CMCC^cm 
I *0  *0*  “ 
I ITHfi-^cnC— COC 


■—■v  I »— t'-f-H-  »C—  I 
, j • «CM  • I 


I -XCMOX  -XCMI 

I 00  coir.o  *0 

C (MCMin  •CM>— CMlfCCM 
I ^ •0"irv  » • •'f'*  • 
I ••me— >tcn'^ir\H-tfc 


'T  O' '»•  o- '!•  O' o- o- O 

cn  pry  ,0  rn  co»  rc>  1 cc  u> 

*ir>mir>inifimcrcifc  *■ 
»'j'o-o'>r>r'ro"OM, 


I 


OC  * •CM  •—4  •»y  »»»•►► 

Z • • t .0 

lOIXOXOXOX  X X 

o 


UJ»- 

I l-^c-l- 

•h^XuL*^ 

•J— H* 

•H* 

h- 

®< 

1 <x<< 

• << 

•< 

• < < 

< 

<X  1 X 

•r  r 

• xrMH-a: 

•rr 

•r 

• r <r 

r 

XOL 

1 QL 

•OLOL 

• QC  ► •QC 

• OCQC 

•QC 

• QC  i-«QC 

oe 

XU 

1 U 

•UO 

• uu«u^u 

•UOU-U 

•UVIU 

u 

rukL 

1 U. 

• XUL 

•LLH-^X 

• UpXXU. 

•L^a.x 

LL 

• cDcD<ZZ>'4Z  Orn  • « • O 

>-_l  • C-  I I •<Oc/1^>s  • *101  < 

C-<XXXuj3X  •nac/100'0'X>-  U- 
-JO^^XXVOvJOcXUJ'-'OCmX  • >--• 

<'-'ZZZ  u'\< -JO-roiOC*— V.  Qi 

V,  OCZ'— 'MVJVJC—VJZ  Z^  u. 

. i—UJXXXt— -4ro  • •— iiJUJUJ<  • 

IC/130XZ  I I 3t<Q.Q.Q.i-  UU 

X •Oci _J_l-l300XOt->>.>'->X  Z 

— 'mXXIIXXXX  — XXXXXXw  t- 
I-  •<JOVJOCOOC30>-0'0®®®<^C-  X 
< CMCMCMCMCMCMrvJCM<tCM^4r-I^^CM<  Q 

XI  ••>••••••  0£ 

OC  I vjoovjooc3C3aeoocMrMrMoix.Qm 

VJ  I UtU>U<UCUCUIUCUC(_)C>CU>U>UtU>U1VJZ— c 
a.  I u-h-p-c— U.UJC/) 

•M  -gcM  •-tcM  -MCMrrc>rmvof~®0'  -gcMm^if\so 


•0 

P*p 

XU' 

rj 

X 

rvj,T> 

'T 

(TV 

U 

X 

tvj 

m 

^r 

X 

X 

ru 

ojrsj 

CNJ 

Csi 

m 

•n 

m 

ro 

<Ti 

0 

•0>0 

0 

>0 

>0 

>0 

>0 

>0 

>0 

*0 

151 


OOOOOOOOOOOOOOOOOOOOOOCJOOOOOOOOOOOOOOOOOC30000000 

i-t  f\j  m ® O' o •-«Ni  'f  vr>  >u  oo  O' o ' 'T  If' >0  r~  00  O' o «-•  f\i  1 •4' lo  vO  00  O' O ^ ovj  fo  m >0  r- CO 

>r  >»■ 'J' >r -J" '!■ 'T  if"f' K' u ' u ' uMr\  ir>  m u \ .O  vO -o -o  >o  -o  >0  •o '0 r~  r~  r~  r~  r'- f'- 00  00  00  00  00  00  OD  00  CO 

U O O O O <J  O U U O U O U OUO  U U O O O O O O O »_)  o <->  (JO  O VJ  O O O O O O O kj  (_)  u u u o uo  o 
a.  a.  a.  a.  a.  a.  a.  a.  o- a.  a.  a a.  a.  a a.  a.  Q.  a.  0.  Q.  o.  cx  Q- o- o- a.  a.  a.  a.  a.  a.  a.  a.  Q.  0.  Q.  a.  a.  a a.  a.  0.  a.  Q.  a.  a.  o. 
oooacooocaoaaaQoouaaaaooooooooooacjooQoooooc/ooooo 


z 

XX 

XX 

^CL 

•» 

» 

H-O 

-> 

XX 

• 

UJ 

UJ 

X 

GO  CD 

cr 

X 

OC 

o 

►- 

o 

-o 

QC 

U, 

X'O 

a 

• 

X 

M 

(S 

u 

D 

o 

1— <. 

h- 

X 

•UJ 

u> 

z 

oo 

X • 

C/' 

nZ 

3 

<>— 

— » 

<x 

XU 

X 

• 

MUO 

D 

X 

QO 

QC 

OO 

O' 

X 

zo 

o 

zu 

in 

►-Z 

z 

OJUJ 

<r 

1- 

XX 

»— • 

z 

a 

O 

•■X 

X 

ax 

M 

X 

Cl 

^ - 

f- 

x^ 

■M* 

^L-» 

UJ 

X 

JVX 

• 

^x 

•iM 

o 

zoo 

►- 

LJ»- 

Z 

Ol- 

•» 

X 

LO 

Q< 

< 

o 

o< 

O' 

XX 

QC 

t/)QC 

M 

M 

kTKZ 

CD 

►- 

H* 

n 

-3 

X 

O 

z 

X 

o 

o •* 

as 

XI- 

-JoO 

in 

X 

M 

Ouo 

fSI 

X 

z 

UJUJ 

mO 

OC 

— XX 

M • 

Zoc 

M 

>cu 

z 

h-rv« 

X 

x>x 

C30 

^x 

UJ 

CL 

u\  - 

»— X 

(✓) 

X^- 

U 

OC 

-iXf- 

M<7V 

LO 

X 

H»  •* 

<✓) 

COK- 

u 

XH-  •• 

M^ri 

-•o 

X 

o< 

z 

x^ 

X 

— o< 

O-^ 

Ox 

3 

xo 

Q 

< • 

DO 

-xo 

acc^ 

H- 

1- 

C£UJ 

♦ O 

z 

X<X 

xD 

Q 

»— 03 

< 

X 

u -o 

OO  • 

XX 

X 

X - 

X 

MM 

o 

M 

<I  - 

I-X 

o 

Or\i 

o 

UJ>- 

0'<) 

xonj 

O -"u 

Q 

►-0 

> 

X 

1 

o 

lOXO 

>— o 

oOX 

• 

m 

>-< 

Ox 

a 

a.^< 

•o> 

QC 

CLX 

o 

oo 

X 

o 

•.XX 

OirN< 

uO 

z 

•»  » 

UOOJ 

X • 

t/1 

M^O 

• 

a?o 

M <1m 

o 

iMcn 

1— •.•x 

X 

u 

QU 

• 

>p-. 

O ♦ 

o 

o 

OQO 

-uOX 

Xh- 

z 

X 

M • • 

z 

XO  - 

• 

< 

•X 

oz 

^ rsjl- 

M 

X 

^ -X 

CjM>r 

X Z 

CMX 

o - 

Ov. 

00*-* 

n 

(NJ 

mO  - 

ou> 

oO  QC 

• 

VH- 

ao: 

O'O 

QO  • 

0> 

M 

in 

aoL 

0<D 

O tx 

as 

X 

Q— 

'/)LU 

^ » 

UJQ^ 

M J 

O' 

o 

*» 

£v/^X 

qCQCm 

. u^l— 

• 

3 

mQ^m 

Z<D 

».XO>X-^ 

O' 

O 

X 

Ov3C 

:bh- 

o 

X 

OQ>^f^ 

rQ 

•■MX, 

>0 

o 

xxo 

- 

zo< 

OC 

MMX  — 

otx 

H-O 

• UJ 

* 

<u 

oox 

xQca 

X 

M 

UXQCX 

OJ  . 

ID< 

• x,ij»acuj 

oo 

•i 

D - 

oox 

X 

X 

M>MQC 

oa 

oo 

OX 

xo 

Q 

o 

xoo 

003 

toxa 

1- 

O M XD 

-IQ 

-J*-* 

0^ 

M 1 

•o 

O 

X 

zxu 

XXh- 

Zm<1 

M 

X 

>M  II  3 

o> 

CD  * 

CDZ 

x*x  H ♦ 

•X 

M 

MO> 

DO  - 

XQCX 

o. 

H- 

# •HmI- 

-OOVZ 

=?o 

^ laoi-i  .oOQCQCOa 

M 

»l/1's.'>,l-i 

Xl-3 

11 

< 

O — — -4  II  -4 

2X<Z« 

</)HM 

O OO  >-40UJUJ4-ZZ 

Zl- 

azx<zzQw 

— 20 

X 

X 

II  a.L0— X II 

ot-3a 

•ODOO  II  <« 

I—Ii.— 0£LU0t< 

001-300 

0X2 

X 

3 

O II  3O3C0 

ZQ. 

►K 

oOO 

O II  O .-m-w  M 

— iLi-ai-oz 

OCZQ. 

•rz< 

x< 

< 

X 

<X— QCh- JC 

ZLJUl^h- 

MM 

II  QCQCf'l  1/1  OC 

X2-* 

CL 

X 

XQCOXXX 

QOJXO^ 

xcx 

LUX  II  II  XXXXXXXXQCXXM^— ^OX^OX< 

XXQC 

LO 

u-M^^^ut 

h-X 

U.M<a30e*MM 

M X o JB  DC  X X u ) u :s o o o 

^Xl- 

u 

U 

— XQCH-^2 

#-^rsj 

• • 

XCNJ 

• • • 

• 

• • 

o 

• • • 

• 

• • 

ooo 

o 

• • • 

• 

• • 

D 

• • • 

• 

(_)000«-)  UOUOO  WJOU 


152 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOCOOOOOOOOOOOOO 

00  O' V/*  O' (J' U' (T' Cr  O' VT  O O O O O O O O O O ^ ^ ^ ^ ^ ^ *Ni  ^ <N  ^ ^ ^ f’ ' 

<vu>jCNCvjrsicacsj(Nj<Nirsjrsjtt)<t>o>(t^rO«»>(t^fO,t^o>ounro«^f»)(^tnfn(nfnfni^nifniii<r>ffi4t^<r>.nnifr>«»»fT>roo»fn 

UOOUOOUO(^OOUUOO^OLJOOOO<^WUCJUOOC^OWO>I^<^OW<^WOU^OOC/O^UJU^ 

Q.o.a.cL^a.^^a.^&^^Qua.o.^a.a.o.Q.a.a.^a.^^Q.^a.a.a.^a.^a.a*a.^a^a.a.a.a>&a.^a. 

oooooooocaooaoooccaooocQCjooQCoocQcooooocococaoccj 


LU 

UJ 

UJ 

a. 

— 

03 

u 

ITS 

X 

•J 

UJ 

> 

o 

< 

<l 

X — 

X 

M 

Ui 

X 

Ul 

UJ 

03 

wo 

'T— 

-D 

z 

M 

<x 

• 

X 

X 

LJ 

u 

«»Mcn 

• 

O 

u. 

—mo 

O 

UJ 

rn 

a 

— O' 

(X 

z 

mx 

o 

M 

UJ 

l~ 

H- 

h^mm^V  *€0 

X 

UO 

UJ 

<xo  ••  .'v^r-r'v'x^ 

h- 

a 

UJ 

o 

Xh-  •<^ 

o 

h- 

< 

•fc  •»  •o^Of^  • 

U. 

(X 

UJ 

•— • 

•» 

CO 

r-m'O  » • • •»  •*>r 

U) 

u# 

— — u»>o  •ujr»uju\'4joo 

h- 

UO 

« 

tNjO'  •c/''Taj  • •'4jm'^  • 

— 

UJ 

O 

> 

> 

JJ 

03 

ac 

<yy 

H‘^0'  • • •••  ••m 'v 

► 

► 

o 

o 

'^<'J’0  *ir\f^aoO'OirvrOi-i 

X 

X 

X 

h- 

u 

o 

— •H- 

xxoof^®  • •r^vO'OOp^ 

> • 

► 

UJ 

« 

■M 

X 

-•X 

^4 

Q. 

cc 

z 

► 

o 

— X 

X r**'*^m  ► •o®0''00'^®^ 

<M 

o 

z • 

h- 

• 

• 

UIUI 

< 

h" 

o 

o 

o-i 

o 

X — 

X •j.j  • • • •f'»m 

< 

w 

ZCU 

-J  — o 

03  UJ 

^ — H^^cUW'CNim'^ 

COZ 

X 

z 

o— 

U.  ^ ^ 

3Q 

^ <^CO  ^ • •Q0f'»p^<\|® 

► » 

► 

• 

0-1 

^ w • 

►-00 

aj 

X 

UJ 

-*•— XO®V-J 

1 ♦ 

LL  ^ ^000  ►OUJ  ••  ► teO  » 

. 

UJ 

•J 

1— 

wO<^ 

^•-•CJXSUJ 

OJH- 

w • • •r-O'coinvj'Ai  •m 

•-m 

H- 

w 

• 

H»  03H~ 

3QCO 

u)— r^p**  •>>^^cJ'  • •oop'»p'*  ►GO 

-X 

X 

z: 

3 

-J< 

i/>i/1>-<3Z  " 

l-<0 

cf\»- 

1— 

3 

30a)-30 

XXlO 

X — ••  ^ — H ^ 

► •> 

► 

O) 

U-  • 

z 

U.UJ 

o 

•-.•-•'s  •^</)0 

•5—  ^t- 

t—  o-JOOOO''7'Ooo(«'flooiri 

N. 

1 

— Qt 

ac 

o 

QOOtD  1 

1 UIX 

<Ui-it—  • • •</«  t •floooi^  •y' 

l/l 

a 

<<X«J«JI  -V 

'-'UJUJ’Q. 

•0*-*  • • •'s.  t 

m 

Hi 

wuu. 

z 

a 

» 

XX  — H-UJ— *oo 

1 QC<Q. 

X 

(13 

Ui 

UJ 

oz 

ru 

II  * OZCDZ2 

00  — X II 

ujoc  ••^•-«»M<rt>f  in.of^®  zy* 

(NJX 

a 

o 

<jOO 

ac 

u< 

H 

— < II  33^0 

UIOOO  — 

Z Z-.0-l_)-l-l_l_J-J_(0_l 

•*x 

X 

UJ 

X 

M 

•-•  X ^O  ►- < X Ui 

(D— ^ — 

— 1— — H“ 

•-«»— 01-0 

II 

c»- 

CO 

o>- 

— 3—  It  H ZZ3^3  II  II  •— 

(_•-•— WX<<<<<<<<X<3 

^Z<m2cm 

z 

t-Z 

03 

o 

O 

v<ot-— cOZo  II  Z II  t-UjUJ—z 

30</)'Ot-ZZZZZZZZI-ZX 

►H 

M 

< 

.HCO 

O 

3 II  o->“ai:oo«-< 

0 II  ooooc 

O'-'ZO  II  0 c 

1 0 

< 

►-0 

UJ 

>-•-■0  oc  t- •«x  — i t- <t- < < ac  3 

OC_JUJI-<  <<<<<<<<<<  — 

P-  ^1- JH- 

X. 

coo 

UJOJ 

It  OXUJXX  LJ^  JXXX®^U)®XX<H>h-»— ^ 

1 ^ 

UJ 

UJO 

U. 

UJ' 

•^<4  — ^WUJU.Xc^c^  ^ VICO^UJ^  — /X'*^X  — 

Ui^OJ^LJ 

at: 

>-< 

O >-i  OC  O »- >-0  oc  ►->  z o 

•-iQ.<tQ.I— aCU.!U)*—lJ  -OOQOCQQfcJLjaOXOOOOO 

ooo 


KJ 

o 


153 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
^ 00  C^  O ^ ^ ^ ^ <30  O' O ^ <M  < ^ ^ IT' ^0  00  O' O ^ J rn  IT' ^ ^ 00  <7' O ^ rO 'T' in 'C  00  O' O f*  <M  rO  nT 

it  1 rn  rn ^ 'T ''T ^ « » u'^  m m m u > in  m in  *n  'O  »<;  ^ 'O r»  r- r»  CD  00  w CO  CO 


OOl->OOOOOOOOUOUOOOOOOOOOOOOUOUOOUUU<->U<-?Ol->UOOOOOUl-^l-> 

o.  a.  Q.  a.  CL  Q.  Q.  Q.  a ^ Q.  a.  Q.  a Q.  Q.  Q.  a.  a.  Q.  a o.  Q.  a.  Q.  a.  Q.  a.  a.  Q.  o.  a.  a.  a.  o.  Q.  o.  a.  o.  a a.  a.  o.  a.  Q.  Q.  Q. 

ooooocooocooaooooaooooaoooaoQooQoaoQaooaaooouooQ 


3 

O 

OJ 


GD’-'O 

JK>u^ 

Q— ‘CO 
^-40' 
ujwfn 
»^xm 


OL 


>>oo 

OOfOO' 

Q.UJCO 

2Q<M^  *0  • 

•*  • •onj 

• •'O  •“ 


<1 

QL 

LU 

O. 

r 


A 

LJ 

M 

> 

>• 

> 

> 

> 

•KNjin  ► *0  • 

1/)  >• 

>• 

> 

>- 

►- 

•> 

•> 

» 

» 

•• 

X 

LOX^O*^  *0 

uo  » 

«> 

» 

< 

X 

X 

X 

X 

X 

X 

«k 

LU_J^  • “O'  “ 

3 X 

X 

X 

X 

X 

» 

•» 

• 

XLUi^icaorgin 

-J  X 

X 

X 

X 

3 

*■4 

m4 

•» 

. 

X •• 

* 

•k 

•> 

•• 

* 

•l 

O' 

J—  *<0  *^00 

UJ  -4 

< 

lO 

in 

O 

r— 

oo 

u 

•— .>%o  “ • • 

♦ » 

» 

«✓> 

u 

<vOO  •.O'OCO 

LO  -4 

(M 

—4 

H* 

t— 

H- 

h~ 

< 

ujw>coinFH<N 

-.00  a. 

a. 

lL 

LL 

< 

< 

< 

< 

< 

< 

z 

^ • • • 

(MfMO  -I 

o 

a. 

a. 

a. 

a. 

a. 

•• 

<0.0^000 

« « OJ 

OJ 

OJ 

OJ 

» 

•» 

9k 

O' 

oO^l^  • » » •> 

♦ ♦ z o 

o 

Q 

o 

z 

X 

X 

L-LUX'O'TT-fM 

<-•01  • 

•k 

o 

o 

u 

u 

o 

o 

» ^OJO^ 

oO  —4 

•.4 

^4 

•— 

»-« 

w«« 

« 

•-9 

1-^ 

• • •vicsj 

OJ#  * 1/)  ^ 

X 

X 

X 

X 

X 

X 

X 

X— > •ooOLu* 

ffl#  — ^ < OJ 

UJ 

OJ 

OJ 

X 

o 

H- 

H- 

h- 

h~ 

h- 

q^'OCD'T  *,  *>  ^ 

1 <1/1  X > 

> 

> 

> 

z 

» 

•» 

•» 

» 

» 

0«— I—  ••>Ot-n»'3QLUOt— wOOZ  “ 

•k 

» 

3 

r'- 

r'- 

r- 

m 

go  O.  O O O O H- »-• 

tUJ<3  “ ^ 

-O 

-a 

■c 

OJ 

IL 

lO 

O.lU^fNicJOO^i**  a 

— 1 Z><J 

< 

U 

•J 

-j 

-j 

>X^VX's.LUsJ-'^'»  O—*  IM_J3 

-i-J 

-J-* 

OJ 

a 

Q. 

a 

a 

a. 

o. 

a 

OJ  3 ^(N^'^ocxin  •►-in  •huq.«-« 

CL -4 

OJQ.<^ 

u 

to 

QC 

OC 

QC 

X 

a 

X 

X 

ZZI-— •xxx<ia3tu 

• ^<r'-0>Xi-iX'-" 

>X>-4X  — 

< 

h- 

O^- 

Oh-G 

H- 

«l- 

«a  “-J— i-ui/i^-xo—o 

►-^OCOK  — 

►-> 

Z«MZrM2(vjZ(M2»'iZ''iZ  — 

f-i-ii— UJLULUUJt-X<  It  It  -lO— XZ>Z>'<t-*2>Z>- 

o 

> 

M 

— ^z 

3Lox>Qaa  II  II  </)i/)(/)x  II  II  »>— • II 

II  >->>*-  II  z 

a 

3 C 

c 

3 

o —x 

OZ  <<  II  i/)L<0  II  _IO»“  II 

1-0 

— • H 

►-QC 

H- 

_!►- -It- _!►- _J»- -it- _)*-_)>  X 

X LU  < < <<  < LU  UJ<  O O X UJ 1- ^ _l  2 3 3 1- _l  Z _J  3 3 

o 

»jO 

-J 

-J 

-J 

_j  II  h-tjcnxh-i— t-t-t-xxi--J— iq:>  ^3>-i_iu 

-joi-az 

OJ 

u u)  vj  o oo  o o o ou  o <-t  u.  (X  uj  i/t  a a a a (_)  o OJ  >/ 1 CO  uj  LU  LU  jc  u X u ^ o jc  o X o 2 u 2 cc  LU  u. 

a 

o 

o 

o 

o 

m 

•r 

in 

<} 

X 

0'<>< 

rg 

fO 

154 


OOOOOOOOOOOOOOOOOOUOOOOOOOOOOOOOOOOOOOOOUOOOOOOO 

m -o  00  O' o "O  O' o «»•  vr»  >o  00  O' o r<i -a"  if\ >0  00  O' o (M  r 1 m -T  If' -o  00  O' 

0000a»00  000'0'O'O'0'0'O'0'O'0'OOOOOOOOOO'^<^'^'-**^'*''^"^'-<*^0'J<'J0'<t'JOtJOOOOOCJO 

o I (T ) f 1 m f' 1 on  ro « »' fo  <»i  fo  fn  ro  fo  ^ ^r '»■  sr  ^ -T 'T  "T -r  O' '*■  O' "T '»■  ■j’ ^r 'T  o- >r 'T  o o o o o o o o o 

>-»  O U U O O O O O VJ  O O (_)  O <_)<-»  O <-J  O U «J  O O'-*  O O <-)  O U O <-)  U O O OO  O O ►->— t- f- h- H- f- 

a 0.0.0.  a 0.0. 0.0.0. 0.0.  a ao.  a.  o-Ol  a.  0.0.  a.  0.0. 0.0.  a.  a 0.0. 0.0.0. 0.0.  a.  a.  0.0.2  ZZ2ZZZZZ 

aooooooooaoooooaooooocooooooooouoococoo^— 


f- 

o 

o 

CL 

X 

UJ 

z 

— • 

u 

ro 

(NJ 

0. 

•If 

• 

£ 

1— 

♦ 

•If 

o 

u 

X. 

a. 

o 

z 

♦ 

lU 

o 

SJJ 

W 

l/i 

u. 

O 

fM 

's,  • 

o 

1 

• 

z 

u 

lU 

tn 

o •* 

M 

UJ 

sO 

.n 

om 

u. 

D 

30 

ir» 

— P-40 

r* 

r* 

z 

< 

• 

.■o 

UL. 

> 

o 

» 

>-o  • 

I— 

1 

<I 

INJ 

UJ(/> 

1 

i/) 

UJ 

^UJ 

♦ 

>• 

eg 

— 4^ 

OCX 

« 

o: 

*^o  • 

> 

• 

<x< 

> 

zo 

X 

in 

<< 

UJ 

ex 

» 

IM 

X ^Cf' 

X 

m 

» ♦ 

z 

u 

i/)U^ 

tr> 

>0 

X • • 

X 

m 

OO 

o 

u 

o 

•OO 

• • 

UJ 

<u 

1 

• 

■i^f^  » 

mCM 

o 

K- 

z 

UJ 

•M 

ro  •m 

o 

-— 

< 

a.a 

eg 

in 

w •00 

ac 

# 

CL 

o 

u 

O 

rO 

□CO  • 

o 

< 

o 

o 

«j 

o 

a<oo 

< 

o 

lf\ 

z< 

• 

V-  • 

o 

►-»— 

1/) 

•— 

LJO 

1 

O •!« 

U'O 

UJUJ 

X 

W>, 

•f 

<7^ 

<Or^ 

a. 

f-H  • 

OO 

r- 

•1-0 

I--J 

X 

U.IO  • 

ro 

> • 

< 

♦ 

^ m. 

a. 

» 

zo 

-JZ 

h- 

O 

X— • • » 

x: 

avin  1 o 

OCQCO 

X 

X* 

oa 

GO 

'M' 

m 

— •oooin 

UJ 

00-4- 

o 

30 

•• 

a 

z 

z • ••»  •O 

• 

eg 

QO 

U 

z 

ZQ 

LL 

• 

«*a  *•  • 

» 

Min'^*  OO 

4.  1 O 

^in 

Q. 

a 

LULU 

o 

h— 

U.Z  *0 

CO 

p>^^^rg>o 

Z-— 

(X 

KWO 

• 

< 

XOJO"* 

'm0^ 

^h-  1 f>j  • • 

rgrg  • 

O'^ 

»— 

Z<1 

1 

00 

uJH-moc 

oOO-^'OOH- 

sg>0O 

«» 

co-r 

z 

< 

-4  CO 

CO 

a. 

Q. 

1-  'sO 

f-OQ.in*^o 

ininm 

o 

■«xo> 

tai4 

o 

• 

X 

ZO.I- 

a> 

• o 

• • 

X 

>-«iz 

•-4 

• • 

ro 

z 

•UJ 

ZUZU 

►-> 

Z •org-*^ 

rsiiNj^  •> 

• X 

<iox 

UJ 

u. 

LU 

~o 

OO  II 

0"-aj<xz  II 

OZOU>f-  — 

inm>o  II 

O M 

0^0 

z 

M 

z 

II  z 

1*4 

■ozz 

•—<-01-0.  II 

—oz 

Z'T'*'  1 'OZZ^OZZI- 

zo 

K 

< 

uj(X 

f-'fu.of. 

h-Z  - 

<CL 

H-  < '.•ou4oe-0'OX«»'u.oc-fu.QC< 

xo 

Z 

z 

cos 

U •^1-3 

O'J0<<-4-ju.O 

0<>-  II  — ►-O  1 1 — 1 1-0  1 l-OZ 

oz< 

UJ 

UK-  ^ 

^QCLJLJX^Q 

o 

O 

• 

II  (XUi  II  II  II  4/#UJ  II 

(X 

< 

OQCUJu.Qca.acuju.LJocjxo^<xuJu.o<o«-^2a&x>i-«i>^cxx^QLULUjaj^ouj 

a 

o 

• 

o 

o 

o u 

o 

o 

IM  in 

•o 

u u o o o •»*o  u <-> 


155 


OOOOOOOOOOOOOOOOOOOCJOOOOOOOOOOOOOOOOOOOOOOOOOOCJO 
o ^ (M  (T)  r- CO  (7>  o CM  ro  vT  m "O  r- 00  O' O «M  cn 'T  If*  >0  00  O' o <NJ  ro  ^ m oo  O' o cj  if» -o 

t-4 1-4  0*4  »M  <\j  O'!  tNj  iM  (M  I'j  fvj  CM  ft  > f * a ' ro  o 1 m r<  n"  n * o 1 v»- .r  «r '*■  >r '*■ 'T  -T  u ' ifMfMn  ifMf*  u ' m 

OOOOOOOOOOOOOOOOOOOCJOOOOOOOOOOOOOOOOOOOOOCJOOOOOO 
acacacai<xccac.a^(X.o^ac.a^.ce.<x.accLoc.oLaLac.aicco^BCocacaC(Xoe:acaLccaL(Xa^aLcx:a^a^al:a^acc£o^a:aia^a: 


z 

UJ 

Zc^ 

z 

0 

XU 

z 

u 

z 

*-z 

z 

a 

4 

U-^ 

o 

ZQ^ 

z 

LUZOU 

z 

z 

uu 

00  — 

h~ 

z 

I«J. 

u 

V)UJ 

00* 

Z'T 

4 

z 

h~  wO>— 

UJU) 

UJ^ 

zz 

U* 

U 

u 

z<u 

ru 

3t/) 

04 

oz 

UU 

0 

z 

.-J< 

—a 

44 

z 

z 

U (/IJ 

UJ 

4 

<z 

>z 

z 

z 

t/lMU. 

z 

^Z-* 

> 

4z 

Z 

z 

4 

UOf— u 

»-ls^ 

z 

ZOJ 

X 

h- 

X « 

>0- 

U 

z 

ZQ 

z 

V^ 

> 

zu 

zz 

^z 

X 

3 

Zh-4 

ZH- 

Zoo 

x> 

u 

H» 

X 

XZUJ 

zo 

Z4 

^ <1 

f-^Z 

h- 

>z 

XU 

<H/>> 

0 

Zw 

0 

KZ 

*.1— UlI 

VJ 

UiZOJ 

0 

z 

8 

OCSUJ 

h- 

Z Z 

ZV) 

ZZ 

Z 

Z 

•-•4> 

xo 

ZV) 

CC  z 

z 

u.  <a 

Im4 

ZKH- 

zz 

uz 

OK 

z 

z 

t->z 

0 

0<N-» 

OH-. 

— 

>-z 

xo 

u 

z 

H-Q 

1-0 

x^ 

oO*-* 

3u 

z 

IJOQ.XI/) 

V^ 

</)Z 

Xoo 

0 

OIZ 

< 

z 

zz 

V) 

Zz 

>— 

-4D 

-i< 

>Q 

z 

0 

• • 

4 

Z-* 

f- 

4 

zu 

z<l 

p- 

~ju.  a. 

UJ 

U) 

zz^-» 

Z4 

V) 

zz 

z 

u 

uC 

UU»-U 

z 

* '**' 

uz 

V>I-N 

— z 

z 

z 

U.  OJ 

X 

v)*-^ 

h- 

00  V) 

II  z 

Z 

tXc/)  1/11/1 

00  UJ 

t/)h- 

z^ 

ZV) 

z 

zz 

< 

CJ 

OJUJ  UJ 

V) 

z 

ZXA 

ZO 

«UJ 

zz 

z 

ZK 

ZK 

z 

>-»-Z3 

z 

z 

UJH- 

zz 

Ol- 

zo 

4 

z< 

• 0 

z 

Z<UJ_t 

UJ 

Z — 

zz 

z< 

Z 

— z 

00  0 

Q 

— 2X 

> 

— z*-« 

z 

oz 

4 

V)  Z 

< 

U 

•> 

LLl 

z 

QO  + 

0 

Q 

Z 

Q 

z 

UJQtJ 

z 

0 

z 

00 

z 

zz 

U Z 

Zo;  UJ 

•» 

< 

Z 

zx 

z 

0— 

►— 

00 

z 

UJ 

«-o<x 

z 

z 

z 

Ooox 

Ol- 

zo 

00 

0 

z 

<✓> 

< 

K-  H- 

• • 

4 

UJ 

< • 

Zw 

>4 

z 

> 

a 0 

u 

X 

n> 

3<jOot 

K 

>- 

Z 

CO 

>00  • 

>z 

CO 

400 

K 

<V) 

z 

QC 

u 

r- 

o<o  •• 

Z 

z 

<00  • 

<a 

Zz 

0^00 

0 

ZZ 

h" 

M 

u* 

QC  u.-» 

Ui 

X 

f- 

z^  • 

z 

oo 

z — 

U 

zo 

•• 

— U 

i. 

z 

ajz  c/) 

Z 

z 

zuz 

z— 

zu 

4U 

4U 

00 

11  r-z 

3UJOH- 

UJ 

-J 

z 

<00  UJ 

<o0 

< 

00 

z 

4 

z 

i-iH- 

• • 

«-o>ZZ 

z 

UJ 

coo 

z 

z> 

ZO 

X 

Z> 

0 

U«i^Z 

•• 

z 

Z 

►-<»- 

< 

«J 

X 

z<z 

ZK 

X 

x< 

K 

X 

4 

II  0 

ac 

UJ 

0 

.0  3 

►- 

z 

z 

<^a 

<< 

»-> 

h-^ 

i-> 

00 

UJ 

H- 

z 

£ uja  • 

VI 

z 

UJ 

z 

Z 

00 

uot 

z 

Z 

QtZZ  — 

z 

z 

z 

•»  0 

z 

UJ 

z 

1-0<QK 

0 

z 

H- 

1 

1 

1 a 

1 

1 

z 

1 

X 

♦ z 

< 

Z 

Z*-»_IUJw 

z 

ac. 

X 

♦ z 

CC 

Ui 

•• 

•• 

M 

UJ 

0 

3 

z 

z z 

0 

UJ 

UJ 

0 « 

Z 

-1 

X 

>— 

z 

0 

> 

u 

a 

z 

V) 

ujZVi/1  II 

•• 

UJ 

0 

z 

z 

< 

a 

X3  1 UJ 

UJ 

< 

< 

X 

z 

z 

a. 

a 

i-a.xa>* 

0 

u 

0 

X 

ru 

f*i 

'4- 

in 

< 

0 

z 

z 

< 

• 

• 

• 

• 

• 

• 

UJ 

CJ 

z 

< 

«i4 

z 

• 

IM 

• 

0*>*OOU<-)  VJ<_IO<->UOU«-)UOc_)  V_>00*-'U000U00U000<~>(->V-l*_)<«4(_)U0U0<->U<-> 


1S6 


VALLES  ARE  PRINTED  kiPERE  FCRMATS  ARE  INCICATEC) 


CJOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOC3000000000 

oo  O' o ro  ^ iri  "O  00  O' o r<^  >r  if> 'O  00  O' o ^iM  -»•  ur\  00  tr  o 'Tin  r- 00  O' o --xM  tn -T  ir\ 

u ' m >u  ^o  ■<J  "O  >o  >*> -o  vo  >*;  r- I'- r- r- r- r~  r~  w flo  oD  oa  00  00  00  uo  uo  oj  O' O' O' O' O' O' O' O' O' O' o o cj  O o O 
OOOOOOOCOOOOOOOOOOOOOOOUOOCJOOOOOOOOOOOOOOUO’H^'-*'-'-*'^ 
ao  a:  Of  QC  a:  a:  oo  oo  oo  oo  a:  oo  oo  ex  a:  oo  oo  oi:  a a:  a:  oo  a:  oo  oo  o' oc  o£  oc  oc:  a:  oc  cc  o£  OL  o:  a:  o:  a o^  oo  a a£  a a a 

H- h- »- K- H- h- 1— h- *- 1— »— t— t- 1— (— f- *- K- h- H- H- K- H- ►- 1- 1- 1— t— h- 1— t— H- 1— (— h- H- 1— t— f- 1— »— t— 

ZZZrZZZZZZZZZZZZZZZZZZZZZZZZZZZZZHZZZZZZZZZZZZZZ 


u. 

LU 

o 

z 

z 

QCUJ 

-^UJ 

■-•tJZI- 

z z 

<V1  •* 

•-HLU 

Qa<< 

< OS  < 

Zuj<v) 

1 

o—  o 

00  UJUJX 

^a.aj-j 

UJ^ 

I-IQOU 

1—  oux<-» 

^o.  <l 

.3  • 

QCLUZf- 

00  z^-  • 

^UJ 

XtX3  • 

uja.LU3 

• 

' J 3 U >- 

o ^z 

<u 

Q.  Q. 

o 

XH-ZXI-3 

-1  O 

V*-* 

LUXZ 

UJ 

3 U 3 

uj^»-Z 

OJ 

UJ 

«30>-' 

QC 

30.30034 

^UJZ 

u<t 

4Z3333 

QC^  -J 

XCJ 

lo  xlulu 

1— 

:d  u 

K-UJ 

I--S  X 

CJI 

— >Q3I- 

</> 

XcOfVJ 

Z l-H- 

UJ 

h-VO 

UJ<K 

QC 

— OCQI-QZ 

<ZUJU. 

o 

w 

LO  O) 

004  4UU 

— - IbJ 

z 

UJOILJ 

lU 

3</>h-^  zr- 

1— UJbJ 

OOOCliJ 

• 

u. 

</>3*-  3 

<a.bQ»* 

s 

S uu 

s 

U.LUI_1X. 

VI 

«4. 

333a.  l/)4 

^ ^UJ 

•J 

^z 

-J 

LU^ULU 

UJ 

X33“-4 

H-O^VO 

a. 

Q. 

QC.Q«l> 

1— 

(/) 

1-333X00 

UJ 

ro  Qto 

ro  ex. 

u •-• 

>• 

Z 

Z 3 3 

K- 

• 

• • H- 

Q UJO 

CD 

<t 

>Z3  31— 

<<< 

z 

(Nj  Z<X 

ujrsj  z 

UJOICO 

QC 

Q03X(/3OZ 

3-JQCU. 

1— K-  X 

OD 

O 

3 |_40« 

030 

UJ  h- 

ZUJ 

<ZOCO 

O 

o 

i/iz  Z3q0 

<0.1— 

A 

O 

O' 

QC 

D4Q  CL 

QC<0 

II  r-  « 

j II 

UU3Q 

(NJ 

u. 

3 

J^LU^LU 

• QC 

1-1 

xa.i_i-i 

UJ 

1/1  •»— 34UJ 

wb-^vn 

wUJ 

UO— 

UJ— •— »uj 

QC  X3 

Z> 

3>-<a.Q0Z 

II  Z 

OJ*-*  If  ZuO 

II  ^ 

lU<ILOL_I 

• • 

00 

Z3Z3I-  — 

a. 

MP« 

Ow  ^ 

l-h-  I 

301—003 

ILZ 

z»- 

JXZI-Z 

LLX^>~ 

Z<LOO>  • 

Q. 

Z 

3X  — zo 

Z300 

• o 

< U"— 

a :3 

— <OI-  LU 

QC 

u 

X001>— 33 

— z o 

I/)  CJ 

^ o 

a 

Z01Z 

H- 

3*— 3 X 

<UJ 

OO  QC 

OCLU 

I-  QC 

UJ  _i*— H’  — 

z 

h— 

30X33 

Q -JJZ 

LU 

X a 

"D 

O OZl- 

U 

3l>0  OX 

UJZXO 

-J  z 

z 

o z 

U.Q.'-'S 

z 

o xa30 

i/1iU3m 

-J  -•»- 

ao  oo 

OC 

o 

033  3 

— > 3 

CL 

< 1/5 

VO 

1—  <0.0^ 

o 

lb 

33«—  O 

>— 

O u 

U 03 

ui  u 

1— 1—  03 

u. 

3<-iXl'iZi/) 

303^ 

UJ 

«-  UJZ 

3 UI 

ZOJ<<3 

QC 

Z3  23<3 

Q 40 

o >- 

1—  f— 

-J  1- 

• • 

UolQI-'/) 

1- 

o 

3m  O O 

3003 

— o 

z— o> 

u 

-«  < 

z 

MO  •■4  ^4 

v)XO 

II  r^uA 

ujr~f^uj^ 

1— <3QO) 

UJ 

(aO 

Za3V'3</> 

M333 

Ql-i  — ►-QC 

h“ 

z 

O < — 

X 

UJ 

3 X— ‘'I '/’</) 

4l/) 

Z^uj 

— -«UJOt 

lb 

o 

zoz^rr 

UJ 

z 

> 3m3 

3lZ3 

II  Q 

II  II  o < 

II  II  Q 

Z 

3ZO  •>•*- 

jC 

Zi/I  1 ZoZ 

Z3  Z 

I 

u.<«— _J 

1-^ 

H- 

o< 

— 003 

-iOC 

«-«3(XX 

«-JQC 

►- 

« 

I-  o 

O 

3Xa:aoaLao 

33Z3 

♦ LJ 

* o 

♦ u 

UJ 

LUZ'-'U.I— 

a 

u 

3333 

3 43 

# QC 

* oeuj 

♦ QC 

QC 

UJ 

xooo 

LiJ 

cc 

3X300300 

OX  *- 

* OC 

* QCI 

♦ QC 

< 

a 

a oi 

QC 

CO 

X43a03ao 

30X3 

S lij 

S oil- 

S UJ 

£X 

•• 

o 

3Q0Z3X3 

03 

lb 

i/» 

Z(_>  OIZ 

UJ 

v> 

04  4 

• « 

00OO4 

o 

HI 

X 

XZOu0<-i 

CL 

QOOUO  00  00  3 

Q 

mqOU 

JC 

I3Z  a 

a 

o 

Ca04O<3 

O 

xxzz 

fNJ 

UJ 

< 

zxh-<a. 

o 

z 

333334 

X 

331/10 

o 

z 

z 

H- 

a 

LU 

• 

• 

• 

• 

UJ 

X 

VI 

QC. 

OJ 

o 

bJ 

X 

OUUOUOVOOOOOUOUOOOOOOUOOUOOOOOOUOOUOOOOOOOOUOOOOO 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OD  O' ^ <30  O' O >rin ® O' O vT  lA 'C  00  O' O ro 'T  m 'O  00  ^ 

>r •^'TiAiAUNiA 

cx.xac.c^cLOLctaiac.cc.cLcLaicLai.ixccai.acaic^cx.cLaccx.acci.oLacacccci.cx.aca'acccccccci.cLacoLix.a^cccccc 

►- H- ^ »— H- ►- f- 1— t- 1— H- e- 1- K f- ►- F- >- 1— »— t— f- h- h- f- ^ I— ►- f- h- 1— ►- 1— h- H- 1— f- 1— *— ►- 1— ►- 1— H- 


U.  UJ 

• 

• 

\ 

UJO  Z xz 

• • 

• 

uja.^-^ 

o 

• 

•-«uj  UJO  < 

UJ 

• 

ooaojK  >-Q^ 

Ul^ 

• 

U^-VJ 

• 

uj^^uiuju.  • 

h— z 

-J 

• 

0(/l  UJH*^>  • 

-i  < 

• 

O » < ^00 

ZZI- 

<J  <_> 

• 

OUi>aJ2^/)  jQ 

I-4U4H- 

z ► u 

• 

>— lOJ  — i 

• 

LL 

H*-J(X£ar-o  j. 

O ^u. 

oo 

• 

z u 

U <i-)aL<z^ 

ZO' 

-.»/1  z 

• 

o 

r—  »00(-J<-J^UJ 

Oojuj 

fXW)  Ul 

• 

>-•031/3 

UJ 

oa  o^za. 

uu^>  • 

U"-  •• 

• 

KZlU 

X >■ 

LUQ^UJV^O 

LL  Z< 

QQ 

• 

U<*.Z3 

►- 

-JX  X z<-i 

UJOX 

OJCVJLUUJ 

• 

z-i-i 

Q 

OJ^f— K-UJ'-'  < 

UJOO^-* 

lr~_n/) 

• 

OCL<I 

U UJ 

^UJZ  X LJO 

UJ<  z 

»— U'»— < 

• 

u.  ^ 

V) 

U)  LJ  h- 

a.x^ 

• 

K 

ULi 

UJ  UJ  < 

^ ouj  UJ  r- 

Z 1— 

V)  1— 

• 

U3  1 LU 

-J 

^ ^ ^ ^ 

-juoQ.  >v)a:< 

o<rzu 

• 

XXX 

-JO  -1  ac 

a.*-*  zoj^ujx 

uxux 

_jLuaJZ 

• 

K »— 

<c  »-•  < a.uj 

(XXZO  QCuOUJ 

oX(/) 

0.03  >~ 

• 

z 

> 

> 1-  > Qt:>- 

<H_uj  •-a<ox 

a<ooo 

O K- 

• 

LL  4 - 

< UJ< 

>»—UJQ.^f- 

o cn 

OK-J-t- 

• 

o — 

X 

>•  -J  X »— -J 

qz^gc:zxo< 

XujC»»« 

UJO*"*'  ' 

• 

ZlO 

a ZQ- 

z^oa*-»a  r 

V-XZX 

lyiOO'U. 

• 

ZlO.-.< 

UJoo 

Ujono.  UJ 

< axuuo 

UJt~<X 

O 

• 

OLU  ✓) 

rt- 

Xf-C^  X ^ 

X cx  ^cx 

X ^ 

ujujuj  • 

t 

>-•  JLOl/) 

l-z 

H>ZUJ  h-  UJLJ 

» •O^UJ  UJ 

U.ZUJ 

• 

h—  JK  — 

•ii- 

h-</)<  ^ •LJX 

UlOUJ 

XKO.(xa. 

• 

o<.zu 

OUJ 

OtJZUJ  OH-  H-UJ 

(/)  H- UJ  O UJ  t/>  O ►- 

X • 

K IX_Ji-i 

• 

Z >— >LO  LO 

za. 

ZCL^-«ac  ZZ  oj 

OZ  ZOf-XO 

l-iZ  03^ 

• 

— O U03  t- 

1—4 

•-4  1—4  *- 4— • UJ 

ZK 

<<C0 

Qt“  » < 

• 

>U.XQ.<  • Z 

x< 

QC4XL0  OCU  OCU 

a<>-‘-0'-»aj> 

UZ^UJ 

O XX 

• 

- o — . — 

OI- 

QKOLU  QIO-  LUK 

►-Q.  z^uxia 

UOCOQC 

uzof-  — 

• 

3cjQt<i  X a 

K4— iQ  1—  X 

< ouj  a. 

OQ-X 

jLU<OX 

• 

-UJLLt— I/)—  0. 

-^OOQ 

lOQX  —I/IQ  2tLU 

ZLUtOuO  OOZ 

ZOQ^O 

O Oi/I 

• 

ZO  <tLLI> 

a 

Zl/>QC  LU  CX 

nijZujH»uj»*«  X 

0*)UJUJ 

—oixo 

• 

•-J  •‘OX  <UJ-JH--* 

_JK  — LUZOC  H< 

UJ>*-cXZXX< 

teH  ^“O 

UJ<JI</lQi 

• 

XoO  ^11  H-f-  Oa 

XK  K -- 

aj^XCL«h->-CX 

H-  «20 

I t-  — 

• 

••>ilJK<  <<ZO.UJZCU<-«.<Z-0  Z — 

QCOOOUJCJ  CJ 

4.uj>-iat 

Kiyi  Q j. 

• 

X 1 H-0>>  UUJOZQOZ  LUUJUJUi  0<^ 

X UJQ0Q.K> 

Z3  Q. 

•— U.Z 

• 

(X  — •— Qt  — 1-i</>K0£  — C >—03 

h-UJH-  <(Z3UJ 

<1/HJL 

a.  OO 

• 

-j-j— Jzxziu»-O<>o  a</)  k4o</>  lok 

XUJ^Z 

-JU1U-J 

UX  CX3 

• 

—oo—  i-iq:3 

ZU. 

ZU.ZZ  ZILL0Z< 

UJ^O^UJUOOO 

Q.»-  <1 

<H/JZ 

• 

.JZ^  U.N*^a.(XUJUOUJCJ**^*-<GLUJO**<UJZ(/) 

UJ  i-.<>UJUJX 

X QO 

ZQCZQ  « 

• 

a«ccu.o*-  ZOX  ZX  OQQX 

a:u.«/)X«Q-ZH> 

uuOOO 

QOO*-»t 

• 

QCt/^UJC  0</)^ 

— — Q-Lt  — — 0L>-<QZ 

oo-ooo*- 

I^X_i 

•— 'O'— 'K*/! 

• 

^ Z QC 

iNfOi/lOOLO  O—OLOujOat  LU 

UJ  2 -*<Z 

QO**- 

k_QiK<lu 

• 

Z<Zt/)UJOUjX 

<z 

lull—  <k  ox 

0</lUJ>(X«/^»-< 

OI^UJZ 

0.0.  <..JCiC 

• 

— ->UJ00LLKOUJU.l/)UJU.Ka  LLLLl/ILULL— UJ 

^^-JO  OJX 

-i  a-u 

« U<_l_l 

• 

X fo  X 

ffl(jtoocj<i  uja30</>i.cj  ^ 

u.<a:o  UJOH- 

•"QC 

a:z«-*a.Q 

• 

ox<zoxa 

•—to 

ZiXO  —<  U0< 

a>:30Lax 

<0120 

o<Za:uj 

• 

UJ  ZZUJlLKIUJt— >.0<>->->LU_4h->Oat:>UJK 

CXO  •-•I^UJLU 

)— (QUJUJ 

lyOQCOUJU 

• 

Z Z '"•-'>  O Qt  00  «/1  < ijO 

4 O J3  4 (jO  < 1/3  < < O 00 

XUJ  K >> 

UJOZ>/) 

UUKXKO 

• 

— • (_l  K a — i 3 QC  03  Z QC  Q1 X > -3  a 03  a 0£  _J 

Oh-ujoa.Ocxcx 

Ok  < 

ooexzQC 

• 

t— — <3(^(J<Q.OXQi<'-' 

izox  x.T^<  a:<z 

<ZI<  ZOD 

o<co 

OO'-'Q. 

• 

Oh-OO  > z 

<wz>w<-4»-<>o 

UJ  1*4  h»  UJ  < < o o 

<OS 

<ULU: 

• 

0<<X  < 1 K- 

X -4 

• 

q^^CDcj^  UJO  (f 

II 

It  II  II  Q.  H H- 

• 

a3u:3<(X-Ja 

H-  ^ 

• 

•JU.v/1  LJOZ 

QC 

• 

u^Qt  Z<^z^-J 

X 

> Z X o>  <l 

• 

UIVIUJ  M 

• 

h"  •**•><*)  to  UJ 

UJ  o 

• 

ZX»^^  X 

X UJ 

• 

• 

uJ 

UO«-)OOOc;*-'‘-J(-)OOOOUOOOUO*-)«->0'-’*-’  0000u00u<-j00<^0*->000<_)0000 


158 


000000000000000000000000000^)00000000000000000000 
'T  m r*- 03 O' o ^ ^ >1’ ^ CO o o rsj  t'o  ^ ^ CO  O' o Cl cf  1 >r  tr> 'O  h- GO  cr  o ^ CV4 m m CO  ^ O 

u > m in  ir\  in  m <3  o >1)  ^ o '4J  ^ 'O 'O '0 c- c- r- r- CD  o)  ® ® GO  CD  oj  00  ® CO  O' O' O' O' O' O' O' O' O' O' o O 

p (i  4#  HfM  ( I ( • i i^CNi  f\J 

QC  QC  o^  o:  QC  (X  a (X  QC  a QL  Q^  QC  QC  o:  cx:  oc  QC  oe:  a 0^  o:  (X  Q^  ct  o:  o:  a ct  a c£  o:  iX  Qc  oc  Q^  o:  QC  (X  0^  q:  (X  Q^  Q^  o:  (X  Q^ 

►- H*  h- *— H- h- H- ►- K I— I- ►- 1— I— h- ^- b- ^ I— I- 1— I- H- H- I- >- 1^  H- h- 1- >- h- )- H- 1— 


:zzz: 


in 

£ 

inni 

<i.Kj 

«>  ^ 

mm 

x^ 

» » 

xu^ 

r > 


-»rsj 

X •- 

^>J* 


< 

> 


>'CL 

<i. 


-*3 

H- 

Q<>0 

< » 

z 

ZL/7 

> ► ..tM 

*1^ 

<LU 

— fN<3 

O 

n 

o 

U 

•^rr)ir>X  " 

(7>0' 

O' 

O-CMIMCJ 

rsi 

fSJ 

b- 

XLU 

Z>Z  »AI 

•k 

r 

— > ►Zv 

C30 

m 

QOOO 

O 

o 

U 

> 

OO'fZ'' 

>—  ►— 

O' 

LU 

b- 

1— 

b- 

(XLL 

o 

*ox^  • 

X 

<»-i 

rsi-i 

uo 

— QUO 

U 

o 

U 

oO 

II  II 

_j  »ro  •■< 

0(JJ 

ooOOO 

O 

o 

O 

L^X 

— — .Jt—oO 

LU 

UJLU 

>T0  ••><  •• 

o 

OX 

►XIMQ^O 

LUb> 

— -Z  -VOZ 

LU 

Z 

-JO  — 

X — - 

«>k 

UJ 

nvi 

»>-4 

1-  o-^ 

rsj-* 

z 

ab- 

XwXw-i^t/1 

-J— J 

z 

o< 

UJ 

o UJO— 

> 

X 

•—• 

-1 1 

Z 1-  • J«.x 

XX 

• 

>Q. 

— *0. 

u.i  UJ  a 

o— 

•Mi^ 

• o • • 

X • 

b- 

Q.*-* 

OH- 

^ ^ 

a o ac 

• Or\j>^UJ 

^ ouJujK' 

»4UJ 

^ o 

o 

o o 

• 

• ViH 

zz  z a> 

►-  • II  1 3 

II  —O  kO-l 

♦ o 

(\J^  • 

m 

m m 

a 

Q.UJ 

OO  1 

LJOI  UJ 

•-<  1 

•JUi 

—— Z 

CLXX-jru  • • 

ox  • 

■b  X 

II  LU 

UJUJ*^ 

— _)  _I  >■ 

1 

O — ZZXru 

«JX  Zhn 

iQ'-ZXO 

0^0 

••-•ao 

• • w 

l/1<  < u- 

ozo 

<M  • 

XI- 

O 30_i33  II 

II  *-'3  II  1-  II  —Z 

>- 

H-a»- 

II  i-i- 

->-)X 

^>  > < 

II  M II  z-JZ-4a.o  II  I 00 

XO<-J  « II -^ni^  M 


XO^ -*ZO  jO  II  Z'^^  rH  -J  (M  w>  II 

X lLLJ.aXX  II  U.XOXU.  II  o II  O M Q II  U.ULa.oO  MlLU.m 
^UJnni>-4  — “3*-<»-*X 


a.WXA.U.VJU.i^O 

U.UIUI 

z 

z 

ULUOCUJ 

« (J 

z 

u 

r“ 

UJi^ 

oo 

<>1 

u\ 

I'^UJU 

.;) 

qt^ 

LJCM 

nj 

nj 

IMCSl 

IM 

iNJ 

«>JXn» 

u 

O. 

X 

X 

O 

X 

o 

u 

u 

o 

uo 

159 


OOCJOOOOOOOOCOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOCJOOOO 

(Nj  fo in CO  o 'O  ^ ® O' o ^ rj  VO  o r*  CD  O' o 'O  00  O' o 00  O' 

nj  Psi  rv  cvj  fvj  cvj  osj  eg  f>j  rsi  rvi  rsj  rsj  rsi  <Nj  <N  c>4  rsi  (N  CN4  cNi  rsj  nj  nj  eg  <\j  <\j  f>j  csj  CNj  ru  f\i  csj  fv  cNj  (Nj  <vj  csj  f\i  CNJ  cvi  cj  fvj  eg  CNi  cvj 
<X,QLaLQLaLQC^XQLCX,CCQLaC<X.aLOCQCCCC£,CLQCQLCLQCOCcLQCCCQC.(:£.^aCCCCt:cCOLQCQCCCa:CCCCaiQC^GCcCOCai 
h- ^ h- H- h- 1- H- 1- h- h- r- »- 1- h- H- ►- J- H H*  h- h- I- I- 1- 1- 

ZZ22Z2ZZZZZZZZZZZ22ZZZ2:Z2ZZ2Z2ZZZ2ZrZ7:ZZZZZZZZZZ 


eg 

>4“ 

X 

1 

-n 

<3 

(ij 

sT 

n 

vi->r 

X 

vj- 

4- 

•n 

in 

m 

in 

oo 

CD 

n 

Q 

a 

O 

O 

O 

O 

H-H- 

X 

II 

eg 

h» 

H- 

►- 

►- 

oo 

a 

o 

u 

o 

a 

a 

o 

oo 

o 

X 

(j 

o 

20 

o 

-CO 

o 

rO 

<— 

'vO 


eg 

< 


>1“ 

< 


eg 

I 


• fooeg— eg«j— 

■I  <^>uj  •egojegv  •lArom 

——XI  *0  1 I X I -fx  I ar'rx^r 


in 

4 


— m H- 

— — m>  z 

m I fO^X  I rriuj 

— I ->  a.z  I — a.(X 
I -j— eg  I uj 

•«*-->-f'-CM  •eg-^O^tL 

ujx  I ^xuj  -)inxi^ 

I ni  ^ ♦ — 


I ^reguJ-)■^  I I inoj  ♦ j • ni  — 

•w^4’0>  •w^in>  •eg>oeg  •w>o>t'Q 

x>xw«i-)x>x— -5X>x--^H-a.»-xx-)X«i.-»-x->X'-*~a. 

II  M II  il  ww  M ii  II  II  w ti  II  II  II  ^ II  II  II  ^ ••  --  ••  >1  >•  • 

^ m fniL  u.  eg  rsjfvjnj  tt.  in  in 'T 'T  li.  o O eg  >r  LL 
x>*^x— X^<  X—XXX—OXOXX^ 


m 

X 

.r\  '^•eg  O 

^ XX  • 
I t v-> 
eg^mm  • 
X •xx>ruj 

# Oww^Z 


II  ii  w II  II  II  ^ II  II  II  _ 

cj  ^li.  in  in  o < 

<1 X u’x  — X ox  o 


— O 
^XX  • 

Q.  I I O 
-j^eg  • 

• XXr^UJ 
O'-'^XZ 

•ininoegiutoi/i  ♦ •ui 

tcDcoeg^  • • •i  •cOifirnix  • 
*--•  < < ^ (/)  O o <-4  w — < < 2 V/)  o 
w it  it  t|  ^ II  II  II  11  w It  It  II  ^ It 
U-eg^n  iu,  egni  ii5iu-fn  ^ iix.  rt> 


(M 

'T 


fn 

>r 


a% 

'T 


>r>r 


GC 

lU 

UOXO 
'T  OUl 

z 


INJ 


160 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOCJCOOOOOOOOOOO 
o ^ (Nj  ro  >»•  in  CO  O' o ir\  >0  00  O' o ^ hJ  00  O' o — < c\j  1’ > >j- ifv  <o  r~- 00  O' o ^ m ^ irv 

m ir\  u Mninmin  m u I m -o 'O -o 'O  ■<J  <0 'O 'O  >0 'W  r- r~  p“  r~  f~- f'- r- lo  00  CO  DO  00  00  oo  oo  CD  00  O' O' O' O' U' O' 
ni  <\j  c\i  (\j  rvi  ro  (\i  (\i  oj  <\i  (\j  tvi  <sj  t\j  (Ni  ro  INI  <\j  <\i  r'j  <\)  rvj  (\i  <Nj  {\i  tNi  <\i  <\j  t\i  f\i  (\J  (\i  rva  fNj  cNj  (\i  c\i  pg  CJ  f\j  (\i  rvj  <\i  c\j  c\i  CNI 
oo  QC  Qi:  oe:  o:  a:  of  oc  a:  a:  a:  a:  o:  (z.  at:  oc  0^  Qc  0^  0^  0^  a QC IX  a:  QC  Qc:  oc  0^  a:  0^  c£  Qc:  QL  0^  00  oo  00  Qc:  oc  Qt: 

I— ^ 1— t— »— I— t— I— t— I— ►- 1— I— I- 1— ►- K- 1- H- I— ►- 1- p- 1- 1— I— t— I— I— ►- >- 1— t— >- 1- 1— ►- 1— 1— I- I— 


zzzzzzzi^; 


uj  a. 

W vx. 

Z H- 

— UJ  Z 

< 

< 

-> 

<y) 

• o 

uOUJ 

< 

< 

uj(y)— . 

«y> 

00 

.Z)  m 

JU,  • UJ 

>> 

'S* 

• •^LJ»'U  z 

^ i-H 

OOl/l  h“UJ  h- 

X. 

ujujxx  •X-J 

i 

1 

—• 

-i-t  U it  r^o 

INI 

-J 

s: 

X 

tO 

I'O 

QCQi<<-0-^  •* 

o 

h- 

UOUUJH^  II  z 

o 

* 

Zro^ 

# 

— — 

UJ<rO 

X 

o 

m 

>r 

— 'v<_) 

Q 

Z<XXQ 

►- 

< 

< 

o— * < 

II  II  UJ>'OlP|LU 

1 

1 

CM 

— • 

a • ..»- 

a 

fM 

tr> 

u-H- 

— o 

_JZ*-xxxo 

ISO 

<l 

KL 

u-  1 — U 

•w 

— z 

OOUJ 

O) 

Lur'l's.Z 

# 

X •» 

V* 

X 

— «• 

rvi 

'T 

W nT-J 

o 

♦ ♦ ♦ * 

o- 

< 

< 

♦ * # 4 

X 

1 

1 

UjrO  + Q. 

♦ •»  •»  •«■  » -GO 

« 

nj 

a 

It  it  o 

0-  — 

X 

1-  1 h-OJ 

ao 

♦ 

> o- 

\r\’¥X, 

o 

m op-O' 

</)XXXX^.^a^ 

Q. 

>r 

•H 

u.x<Mr- 

X 

cr* 

U'  U'U'V/' 

UJ 

m -J 

X 

X * 

a — X 

Q 

o 

o 

o ooo 

Zrsjcsjrvjin 

Z • 

fO-—  o 

«rfM 

« 1 1X0^ 

<NJ 

<M 

<NJ  ^N*<M(\J 

» OJ 

<*  -o 

■<HKX 

zomocLO 

*> 

xxxxxxx-o 

mo^nz 

♦ m 

•^»r>  t 1 

a ‘X  1 

CF 

• •'TP'l  1 Z^  II 

wCT' 

wO' 

O-H  — — 1 

i-<ox><x>-^-<ox>-<XK«-'-ao 

Z<l-»-l-KI>-l- 

It  It  It  W 

II  II  II  II  II  II  II  n 

II  II  It  II  II 

It  <11  II  >»  II 

UJ  o UJO  UJOUJ  UJUJ  ^^- < < < < < < 

-T  *^u.n»<'rrrirornfOCJ'r<.r(rrirninfnZtMiT»>-X~3Xh*>-t-f->-t->-i-l— Oi/l3:zZXZZ 

—•  •->  «•-•*—(—  OOQOQOaCQOaO  OJ 

X I—  uj  aoLjaoLjaou;xaoaoujr-Lj(_]i_iLj(_ii_i  ^ 


CL 

LU 

a 

a 

X 

m 

r- w 

xcto 

(XO 

4n 

xo^^uO'or^O' 

UJ  AJ 

or- 00 

xwr 

O' 

U'U'U' 

UJ  v' u' v' y' 

O 

u 

UJ 

U.OOOOOO 

r\irsi(M(sj(M<N< 

U 

u 

o 

O 

161 


BIBLIOGRAPHY 


1.  Heat  Exchange  Institute,  Standards  for  Steam  Surface 
Condensers  , Sixth  Edition'^  1970  . 

2.  Department  of  the  Navy,  Bureau  of  Ships,  Design  Data 
Sheet  DDS4601-  1 , 15  October  1953. 

3.  Oak  Ridge  National  Laboratory  Report  ORNL-TM-  4 248  , 

ORCONl : A FORTRAN  Code  for  the  Calculation  of  a Steam 

Condenser  of  Circular  Cross  Section,  bv  J.A.  Haf ford , 

July  19^3. 

4.  Eissenberg,  D.M.,  Ai  Investigation  of  the  \'aria~le< 
Affecting  Steam  Condensation  on  t!'.e  ufside  o:  a 
Horizontal  Tube  Bundle,  Ph . D . TTTesTs  , 'ini'^er-.* 

Tennessee , December  T??  2 . 

5.  Search,  H.T.,  A Feasibility  Study  of 

Improvement  in  Marine  Steam  Condenser:.  ''  ■■ 

Naval  Postgraduate  School,  December  i-" 

6.  Stoecker,  W.F.,  Design  of  Thermal  5 :.’::^,: 

Hill,  1971. 

7.  Vanderplaats  , G.N.,  Automated  Design  jr  r ' i : j • . 
class  notes  for  a graduate  course  of  t;c 
presented  at  the  Naval  Postgraduate  Schoo.,  M*  . - : r . 
California,  March-May  1977. 

8 . Fox  , R . L . , Optimization  Methods  for  En^jineering  "e  i ;n  . 
.A.ddison-Wesley  , li57 1 . 

9.  NASA  Ames  Research  Center  Technical  Memorandum  NASA 

TM  X- 6 2, 2 82,  CONMIN  - A FORTRAN  Program  for  Constrained 
Function  Minimization  User's  Manual  , by  G.N.  \'anderp laa ts  , 
August  1973 . 

10.  Vanderplaats,  G.N.  and  Fuhs , A.E.,  "LASTOP  - .A  Computer 
Program  for  Laser  Turret  Optimization , " NPS  Report 

NPS  69-77-0004,  Naval  Postgraduate  School,  1977. 

11.  Fletcher,  R.  and  Reeves,  C.M.,  "Function  Minimization 
by  Conjugate  Directions,"  British  Computer  Journal, 

V.  7,  n.  2,  p.  149-  154,  1964"! 

12.  Zoutendijk  , G.G.  , Methods  of  Feasible  Directions, 
Elsevier,  Amsterdam"!  i960 . 


162 


15.  Vanderp laats , G.N.  and  Moses,  F.,  "Structural 

Optimization  by  Methods  of  Feasible  Directions," 

Journal  of  Computers  and  Structures,  v.  3,  p.  739-755, 

19 

14.  Aerodynamic  Analysis  Requiring  Advanced  Computers 
NASA  SP-547  Part  11,  Application  of  Numerical 
Optimization  Techniques  to  AirfoilPesign,  ^ 

G . .N . Vanderplaats  , R . N' . Hicks  and  E.M.  Murmaa,  p.  749- 
768,  March  1975. 

15.  Vanderplaats,  G..\'.,  COPES  - A User ' s Manual , prepared 
for  a graduate  course  on  "Automated  Design  Optimization" 
presented  at  the  Naval  Postgraduate  School,  Monterey, 
California,  March-May  197"'. 

lb.  Harrington,  R.L.  (ed)  , .Marine  Engineering,  p.  4"", 

The  Societv  of  Naval  Architects  and  Marine  Engineers, 
1971. 

17.  Baumeister,  T.  and  .Marks  L.S.  (eels).  Standard  Handbook 
for  Mechanical  Engineers,  'th  ed..  Chanter  3,  pT  55-6'' , 
McGraw-Hill , 1567 

18.  Crane  Company,  Technical  Paper  No.  410,  1976. 

19.  Naval  Ship  Systems  Command  0 946-004  -7010,  .Main  Condensers 
(CV.A  6 7)  , January  1969. 

20.  Westing ho use  Electric  Corporation  Steam  Divisions 
Technical  Manual  1440-C72  (NAVSHIPS  346-0268),  6,600 
Square  Foot  Surface  Condenser  and  Main  .Air  Ejectors, 

May  1964 . 

21.  Holman,  J.P.,  Heat  Transfer,  4th  ed.,  .McGraw-Hill,  19'^6. 

22.  United  States  Atomic  Energy  Commission  Report 
OR.NL-TM- 297  2 , Computer  Model  and  Correlations  for 
Prediction  of  Horizontal  Tube  Condenser  Performance 
in  Seawater  Distillation  Plants,  by  D..^(.  Eissenberg 
and  H.M.  Noritake,  April  1970. 

23.  Universal  Oil  Products  Company  Bulletin  .No.  4020, 
Wolverine  Korodense  Tube:  General  Information  and 
Heat  Transfer  Design  Notes,  10  January  19'?3. 

24.  Fraas,  A.P.  and  Ozisik,  M.N.,  Heat  E.xchanger  Design, 
Wiley,  1965. 

25.  English,  J.W.,  Hydrodvnamic  Considerations  in  the  Design 
of  Condenser  Cooling  Water  Systems  for  Large  Ships, 
paper  presented  at  the  Royal  Institution  of  Naval 
Architects,  London,  England,  12  April  1973. 


163 


J 


26.  Reilly,  D.J.,  An  Experimental  Investigation  of  Enhanced 
Heat  Transfer  on  Horizontal  Condenser  Tubes,  MSME 
Thesis,  Naval  Postgraduate  School,  March  197 8 . 

7.  Westinghouse  Electric  Corporation  WAPD-TM-213, 

STDY-5  A Program  for  the  Thermal  Analysis  of  a 
Pressurized  Water  Nuclear  Reactor  During  Steady-State 
Operation , by  Pyle,  R . S . , June  1960  . 


I 

i 

f 


4 

r 


164 


INITIAL  DISTRIBUTION  LIST 

No.  Copes 

1.  Defense  Documentation  Center  2 

Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0142  2 

Naval  Postgiaduate  School 

Monterey,  California  93940 

3.  Department  Chairman,  Code  69  2 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Office  of  Research  Administration,  Code  012A  1 
Naval  Postgraduate  School 

Monterey,  California  93940 

5.  Professor  Paul  J.  Marto,  Code  69Mx  10 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 

6.  Dr.  Gary  N.  Vanderplaats , Code  1613  5 

David  Taylor  Naval  Ship  Research 

and  Development  Center 
Bethesda,  Maryland  20084 

7.  Professor  Paul  Pucci,  Code  69Pc  1 

Naval  Postgraduate  School 

Monterey,  California  93940 

8.  Lt.  Charles  M.  Johnson  4 

3784  Hyak  Way 

Bremerton,  Washington  98310 

9.  CDR  N.P.  Nielsen,  USN  1 

Naval  Sea  Systems  Command  (033) 

2221  Jefferson  Davis  Hwy , CP#6 
Arlington,  Virginia  20360 

10.  Mr.  Charles  Mil?er  2 

Naval  Sea  Systems  Command  (0331) 

2221  Jefferson  Davis  Kwy , CP#6 
Arlington,  Virginia  20360 


165 


1 


11.  Mr.  Frank  Ventriglio 
Naval  Sea  Systems  Command  (0331) 

2221  Jefferson  Davis  Hwy , CP#6 
.Arlington,  Virginia  20360 

12.  Mr.  .Arthur  Chaikin  1 

Naval  Sea  Systems  Command  (0331) 

2221  Jefferson  Davis  Hwy,  CP#6 
Arlington,  Virginia  20360 

13.  CAPT  J.  K.  Parker,  USN  1 

Naval  Sea  Systems  Command  (PMS-301) 

2221  Jefferson  David  Hwy,  CP^^6 
Arlington,  Virginia  20360 

14.  CDR  D.  W.  Barns,  USN  1 

Naval  Sea  Systems  Command  (PMS-301. 3) 

2221  Jefferson  Davis  Hwy,  CP*6 
Arlington,  Virginia  20360 

15.  Mr.  Walter  Aerni  1 

Naval  Ship  Engineering  Center  (6145) 
Washington,  D.C.  20362 

16.  Mr.  Wayne  L.  .Adamson  1 

•Naval  Ship  Research  5 Development  Center 

(2761) 

•Annapolis,  Maryland  21402 

17.  .Mr.  Gil  Carlton  1 

Naval  Ship  Engineering  Center  (6'23) 
Philadelphia,  Pennsylvania  19112 

18.  Dr.  David  Eissenberg  1 

Oak  Ridge  National  Laboratory 

Post  Office  Bo.x  Y 

Oak  Ridge,  Tennessee  37830 

19.  Mr.  Joseph  A.  Hafford  1 

Oak  Ridge  National  Laboratory 

Post  Office  Box  Y 

Oak  Ridge,  Tennessee  37830 

20.  Mr.  John  W.  Ward  1 

Marine  Division 

Westinghouse  Electric  Corporation 
Hendy  Avenue 

Sunnyvale,  California  94088 

21.  Mr.  Henry  Braun  1 

De  Laval  Condenser  and  Filter  Division 
Florence,  New  Jersey  08518 


166 


1 


: 22.  Miss  Eleanor  J.  Macnair 

Ship  Department 
Ministry  of  Defence 
Director  - General  Ships , Block  B 
Foxhill,  Bath,  Somerset 
E.NGLAND 

23.  Mr.  Kurt  Bredehorst  1 

MAVSEC  6147D 

I Department  of  the  Navy 

Hvattsville,  Maryland  02782 

i 

t 

I 


I 


